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MULTIPLE OTMPONENT RNA CATALYSTS AND USES THEREOF 
This application is a continuation-in-part of U.S. Serial 
5 No. 08/378,235, filed January 25, 1995, the contents of 
which are hereby incorporated by reference into the present 
application. 

Background of the Invention 
10 The invention disclosed herein was made with Government 
support under NIH Grants No. GM 50313-01, ROl G50313, from 
the Department of Health and Human Services. Accordingly, 
the U.S. Government may have certain rights in this 
invention. 

15 

Throughout this application various references are referred 
to within parentheses. Disclosures of these publications in 
their entireties are hereby incorporated by reference into 
this application to more fully describe the state of the art 
20 to which this invention pertains. Full bibliographic 
citation for these references may be found at the end of 
this application, preceding the sequence listing and the 
claims. 

25 In 1977 a discovery was made which radically modified our 
understanding of gene expression. Two laboratories observed 
that the protein coding sequences of an adenovirus gene are 
interrupted by regions of encoding DNA (Berget, 1977; Chow, 
1977) . The sequences of these interrupting regions, or 

30 "introns", are transcribed into precursor mRNA, and excised 
from the parent RNA in a process called RNA splicing. The 
coding ends of mRNA are then joined together, resulting in 
"processed" mRNA which is ready to code for protein. 

Comparison of genomic and RNA sequences revealed that 
35 introns are actually quite common in eukaryotic genes. Many 
genes contain multiple introns, and by splicing them out and 
religating the coding regions together in different 
combinations, a single gene can produce multiple protein 
products. This imparts an additional degree of control and 
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versatility to the expression of eukaryotic genes. 
It was subsequently learned that a host of proteins and five 
nuclear RNAs (snRNAs) assemble as a complex upon nascent 
transcripts of pre-mRNA. This protein/RNA assembly was 
5 found to be responsible for catalyzing the splicing reaction 
and was appropriately dubbed the w spliceosome ,, (Brody, 
1985) . Analysis of reaction products showed that the 
spliceosome facilitates nucleophilic attack of an adenosine 
within the intron upon a phosphodiester bond at the 5'- 
10 intron /exon boundary to generate a branched intermediate 
(Fig. 1). The liberated 5' -exon then attacks at the 3'- 
exon- intron boundary, resulting in ligation of the two exons 
(thus generating coding mRNA) and excision of a lariat - 
shaped intron (Guthrie, 1991) (Fig. 1) . 

15 

Not long after the discovery of RNA splicing, workers in the 
laboratory of Thomas Cech observed that precursor ribosomal 
RNA from Tetrahymena thermophila was capable of undergoing 
an in-vitro splicing reaction in the absence of protein 

20 enzymes (Kruger, 1982). Although it was assumed that a 
complement of endonucleases and ligases were involved in the 
two transesterifications required for splicing, the Cech lab 
showed that a structure within the folded intron itself 
catalyzed the splicing reaction in the presence of Mg2+ and 

25 a guanosine cof actor. This was the first example of a 
catalytic RNA, or ribozyme, in which RNA itself was capable 
of behaving like an enzyme. It appeared that "self- 
splicing" by Tetrahymena intron (phylogenetically classified 
as a group 1 intron) was fundamentally different from pre- 

30 mRNA splicing because of structural differences in their 
reaction products (Cech, 1986) . 

Several years later, members of the group II family of 
introns were observed to undergo a self -splicing reaction 
35 in-vitro (Peebles, 1986; van der Veen, 1986) . Group II 
introns are unrelated to the group I introns in sequence or 
secondary structure and, although group II introns undergo 
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a self -splicing reaction, their reaction products are most 
similar to those of pre-mRNA splicing (Fig. 1) . Like pre- 
mRNA splicing, the first step of group II self -splicing 
involves attack of an internal adenosine residue rather than 
5 a bound guanosine, and after splicing is complete, group II 
introns are excised as lariats. Additionally, the sequences 
surrounding the splice sites near group II introns and pre- 
mRNAs are very similar. For these reasons, group II introns 
have often been referred to as the "missing link" between 
10 group I introns, which are believed to represent a sort of 
molecular fossil of primitive RNA processing, and the 
splicing of pre-mRNA which is a relatively modern 
development (Cech, 1968; Jacquier, 1990) . It has been 
proposed that the subdomains of the group II intron are 
15 analogous to snRNA molecules within the spliceosome and, 
since characterization of the unimolecular group II ribozyme 
would be far easier than dissection structure and mechanism 
of the spliceosome, study of group II intron ribozymes 
provides an important basis for our understanding of RNA 
20 splicing in general. Additionally, a functional analogy 
between group II intron subdomains and snRNAs would support 
the theory that eukaryotic pre-mRNA processing is catalzyed 
by RNA components of the spliceosome. 

25 Instead of proteins and snRNAs, the active site of the group 
II intron is composed of a single strand of folded RNA. 
This RNA is divided into six domains (Fig. 2A-2B, Fig. 3), 
which may, like snRNAs and proteins, function separately in 
composing the active site and catalyzing the splicing 

30 reaction (Michel, 1989) . It has therefore been of interest 
to define the functions of the individual domains. 
Techniques of molecular genetics have shown that Domains 2, 
4 and most of 3 can be replaced by short hairpin loop 
regions without detectable effects on self -splicing (Bachl, 

35 1990; Kwakman, 1989) . Domain 1, which is the largest, 
recognizes the 5'-exon by base-pairing and participates in 
several other important stabilizing interactions {Jacquier, 
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1991) . The most conserved region of the group II intron and 
the one with the greatest effects on Bplicing is Domain 5. 
This short hairpin region of RNA (only 34 nucleotides in the 
ai5g intron) appears to direct the folding or activity of 
5 the group II intron in a very fundamental way (Jarrell, 
1988; Pyle, 1994) . 

In 1988, the laboratory of Philip Perlman found that a group 
II construct containing only Domains 1-3 was unable to 
10 undergo the first step of splicing. However, addition of 
the tiny Domain 5 in-trans rescued the splicing reaction 
(Jarrell, 1988) . Direct proof of the strong D5 tertiary 
interaction was provided by kinetic and thermodynamic 
studies showing D5 bonds in trans with a Kp of ~300nM. (Pyle, 
15 1994) . Given that Domain 5 is almost completely paired and 
there is no phylogenetic covariation between sequences in 
Domains 1 and 5 (Michel, 1989), it is highly unclear how it 
would associate with Domain 1 to generate the catalytically 
active structure. 
20 With the finding that group II intron structures are formed 
during the trans -splicing of mitochondrial energetics genes 
from higher plants (Wissinger, 1992), it is now believed 
that group II splicing reactions are of general importance. 
They are central to the metabolism of plants, yeast and 
25 possible other eukaryotes. Some group II introns contain 
open reading frames which code for "maturases', or proteins 
which assist in the splicing and reverse-splicing of group 
II introns in-vivo (Perlman, 1989) . Some of these maturases 
contain reverse transcriptase and RNase H domains 
30 (Lambowitz, 1990) . These findings have precipitated 
interest in the introns as mobile genetic elements, or 
"infectious" RNAs, readily able to integrate into new 
genomes (Lambowitz, 1989). Interestingly, open reading 
frames and base-pairing recognition domains for trans- 
35 splicing variants of the group II intron are usually found 
in Domains 2 or 4 of the group II structure, regions of high 
variability and low sequence conservation (Michel, 1989) . 
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Although group II intron reactions are important in their 
own right, the greatest interest in group II introns stems 
from the hypothesis that they are related chemically and 
structurally to the spliceosomal apparatus responsible for 
5 eukaryotic pre-mRNA processing (Sharp, 1985) . This theory 
implies that the group II intron may be used as a model for 
reactions which take place in the spliceosome, a reaction 
center which is extremely difficult to characterize because 
it contains hundreds of components and must be assembled in 
10 cell extracts. It also suggests that group II intron 
chemistry will describe "remnants" of RNA catalysis still 
present in higher biological systems and suggests that 
spliceosomal processing is RNA- catalyzed despite the 
presence of proteins. There is a great deal of chemical 
15 theory which now supports these ideas. Work by Moore, Sharp 
and coworkers on the stereochemistry at pre-mRNA splice 
sites has provided evidence for the involvement of two 
active sites in spliceosomal catalysis (Moore, 1993) . This 
suggests that group II introns may also have two active 
20 sites; one for the first step and one for the second step of 
splicing. Group I introns are generally believed to have 
only one active site, and that is why the stereochemistry of 
nucleophilic attack in the first step is opposite to that of 
the second step (McSwiggen, 1989; Rajogopal, 1989) . 
25 Spliceosomal processing and group I introns are most similar 
in the second step of splicing, where both prefer the pro-S 
isomer of phosphate. Spliceosomal processing and group II 
introns appear most similar in the first step which is 
initiated by adenosine 2' -OH attack. This implies that the 
30 first step in group II intron splicing is most likely to 
represent a new form of catalytic active site and one which 
we must characterize before we can understand the 
relationship to pre-mRNA processing. To this end, 
engineering multiple- turnover ribozymes from the group II 
35 intron will enable us to examine the chemical mechanism for 
reactions catalyzed by the active site(s) of the group II 
intron. This same approach was a prerequisite to 
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characterization of the transition state, the active-site 
structure and the detailed mechanism descriptive of the 
group I intron (Cech, 1992) . 

5 it was originally believed that in-vitro group II self- 
splicing was inefficient before appropriate reaction 
conditions were defined in 1988 (Jarrell, 1988) . Despite 
splicing reactions which are complete within minutes and the 
fact that splicing is a complicated multi-step process, the 

10 belief that group II introns have lumbering reactivity still 
persists. It is also commonly believed that the structure 
of the group II intron is metastable, with few molecules in 
a pool assuming an active conformation. Unfortunately, 
there is little quantitative data to substantiate either of 

15 these theories. Many generalized ideas about the nature of 
the group II intron come from observations of a single 
molecule, particularly the ai5g intron from yeast 
mitochondria. There have been no systematic comparisons of 
group II intron reactivity of stability, and virtually no 

20 quantitative analysis on the kinetic behavior of any single 
intron. We know nothing about the architecture of the 
folded structure. Clearly, before the group II intron can 
be discussed coherently as a model for pre-mRNA processing, 
we must characterize much of its fundamental behavior. As 

25 indicated above, the most rapid way to do this is to 
transform the uniraolecular group II self -splicing reaction 
into multiple-turnover ribozyme that function under simple 
conditions. This will impart significant control over the 
reaction and facilitate careful mechanistic analysis of the 

30 intron. 

Several research groups have been pursuing the use of 
catalytic RNA, also known as ribozymes, for uses in 
diagnostics, as restriction enzymes, see for example Cech et 
3 5 al. U.S. Patent No. 4,987,071. Ribozymes which are highly 
specific in the recognition sequences offer potential in 
therapy for plants or mammals, see for example Haseloff et 
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al. U.S. Patent No. 5,254,678 for hammerhead ribozymes, 
Jennings et al. U.S. Patent No. 5,298,612 for minizymes; 
Hampel et al . European Appn. No. 

89 117 424.5, filed September 20, 1989 for hairpin 
5 ribozymes; and Altman et al. U.S. Patent No. 5,168,053, for 
prokaryotic RNAase P. However, this is the first report of 
the use of a modified group II intron composition to act in 
trans, i.e. as a catalyst, to cleave a nucleic acid target. 
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SBfflHiaPC of ^ he Invention 

This invention is directed to a composition for catalyzed 
oligonucleotide cleavage comprising a synthetic non- 
naturally occurring oligonucleotide compound. The compound 
5 comprises nucleotides whose sequence defines a conserved 
group II intron catalytic region and nucleotides whose 
sequence is capable of hybridizing with a predetermined 
oligonucleotide target sequence to be cleaved, such target 
sequence not being present within the compound. The 

10 composition also includes an appropriate oligonucleotide co- 
factor. Preferably, the conserved group II intron catalytic 
region is a group II intron domain I catalytic region. In 
one embodiment the conserved group II intron domain I 
catalytic region may further comprise a conserved portion of 

15 a group II intron domain II, a group II intron domain III, 
a group II intron domain IV, a group II intron domain V, or 
a group II intron domain VI. The invention is also directed 
to methods of treatment and methods of use of such 
compounds . 
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BRIEF DESCRIPTI ON OP THE PT^ ttQ 

Figure 1: Shows a comparison of RNA processing by self- 
splicing introns and the spliceosomal apparatus. 

5 Figure 2A and 2B: Shows the secondary structure of the two 
group II introns from yeast mitochondria (seq ID No. 4 and 

5) . 

Figure 3: Shows a schematic of a group II intron. The 
10 structure is composed of six subdomains organized like 
spokes of a wheel. Domain 1 contains the EBS1 and EBS2 
regions (Exon-Binding Sites 1 and 2) which base pair to the 
IBSi and IBS2 segments of the intron. The first step of 
splicing is initiated by attack at the 5' -splice site from 
15 an adenosine 2' -OH group on Domain 6. Domains 2-4 are 
dotted because they are dispensable for this first step. 

Figure 4A and 4B: Cleavage of a short RNA substrate (S) by 
the group II ribozyme. (Fig. 4) Domain 1 (Dl) combines with 
20 Domain 5 <D5) and substrate (S) to catalyze specific 
cleavage at a site analogous to the 5 '-splice site. Shaded 
regions and dashes between Dl and S designate pairing 
interactions IBS1:EBS1 (closest to the cleavage site) and 
IBS2:EBS2. Dl RNA (425 nts) was transcribed from plasmid 
25 pT7Dl cut with Eco RI . This plasmid was constructed by PCR 
amplification of Dl from the full-length plasmid pJD20 
(Jarrell, 1988). Primers specifically amplified a 425 nt 
intronic fragment including all of Dl beginning with the +1 
nucleotide of the intron (G) . The upstream primer included 
30 a T7 promoter sequence and a BamHI restriction site. The 
downstream primer contained an EcoRl site. The PCR fragment 
was ligated into a pUC19 vector, expressed and sequenced. 
D5 was transcribed as previously described from plasmid 
pJDI5'-75 (kindly provided, together with pJD20 by Dr. 
35 Philip S. Perlman) (Pyle r 1994). Transcribed substrates 
were made using T7 RNA polymerase using synthetic DNA 
oligonucleotide templates according to published procedures . 
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All RNAs were purified as described (Pyle, 1994). (Pig. 4) 
Ribozyme cleavage of substrate <S) 5'-"P- 

(SEQ ID NO 1) to a 5' -end- 
labelled 19 nt. product (p) and an 8 nt. unlabelled product. 
5 The cleavage site is marked with a " Vff and IBS sequences are 
underlined. Lanes 1-8 correspond to a reaction of S(l nM) 
in the presence of saturating Dl (100 nM) and saturating D5 
(3 /iM) . The reaction was performed under standard reaction 
conditions: at 42°C in 80 nM MOPS (pH 7.5), 1M KCl, and 100 

10 nM MgCl 2 . RNAs were mixed at 1.67x final cone, in 80 nM MOPS 
pH7.5 and then heated to 95°C for 1 min. to denature 
potentially mis-folded RNAs formed during storage. After 
cooling to 42°C, reaction was initiated by addition of 100 
nM MgCl 2 , IM KCl and 80 nM MOPS (final concentrations). 

15 Aliquots of the reaction were removed at the time intervals 
specified in the figure, combined with denaturing dye on 
ice, loaded on a 20% denaturing polyacrylamide gel (PAGE) 
and quant itated as previously described. 

20 Figure 5A and SB: Single -turnover kinetics of S cleavage by 
the three-part group II ribozyme. Rate constants and 
kinetic parameters (Table 1) were obtained by plotting 
pseudo- first order rates (k^) , determined from the slope of 
inset plots against a react ant varied in concentration. It 

25 was determined that k ob , values at saturation did not vary 
from 0.01 to 3 Nm [S] , indicating that pseudo-first order 
kinetics apply (data not shown) . Pseudo-first order semi- 
log plots of the cleavage timecourse were linear for >3 
reaction half-times (see also Fig. 9B) . Early time points 

30 (generally <50% completion) were used in the calculation of 
k ob . (see also Fig. 9B) . Individual values q£ a k were 
determined over a 6 week time period with the same RNA 
stocks and varied by 18% (with 95% confidence and 5 trials) . 
The 24-nucleotide substrate, CGUGGUGGGACAUUUUC^GAGCGGU (SEQ 

35 ID NO 2), was made synthetically (Scaringe, 1990). The 
empirically determined endpoint for reaction of this 
substrate was 80% after >7 reaction half-times. This value 
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was used in correction of all pseudo-first order plots used 
in this figure. (Fig. 5A) Inset: Representative pseudo- 
first order plots for cleavage of S (0.05 nM) at saturating 
D5 (3 mm) and varying concentrations of Dl (0.5, 1.5, 3.25, 
5 5, 10, 20, 100 nM shown). Hyperbolic nlnt- • The k obfl values 
determined from inset plots were graphed against [Dl] to 
generate the apparent binding curve for Dl . In both A and 
B hyperbolic plots, the k^ value for the reaction is the 
horizontal asymptote of plots <~ 0.02 min' 1 ) , and Km is the 

10 concentration of enzymatic subunit at which k ob . is one-half 
k roM . (Fig. 5B) Inset : Representative pseudo-first order 
plots for cleavage of S (0.1 nM) at saturating Dl (25 nM) 
and varying amounts of D5 (0.2, 0.33, 0.67, 1, 3, 9 fi.M 
shown) • flyperbplip plpt; The k ob8 values determined from 

15 inset plots were graphed against [D5] to generate the 
apparent binding curve for D5 and the kinetic parameters 
shown in Table 1 (line 3) . 

Figure 6: Ribozyme cleavage of a single phosphorothioate 
20 linkage: Comparison of ribozyme cleavage of a transcribed 

all -phosphate RNA substrate (squares: connected by ) and 

a synthetic substrate containing a single phosphorothioate 
at the cleavage site (circles; connected by - - -). a 
standard ribozyme digestion was performed as in Figures 4 
25 and 5, except that, unlike these previous kinetic figures, 
final [KC1] was 2M in order to approach endpoint more 
rapidly. Doubling the KCl concentration increases the rate 
of reaction approximately ten- fold, while pseudo- first order 
plots remain linear (data not shown). The all -phosphate 
3 0 substrate was cleaved to 96% ±0.8% (with 95% confidence and 
4 trials), empirically determined after >7 half-times. This 
substrate was cleaved 52% ±4% (with 95% confidence and 5 
trials) at >7 half-times (see Fig* 7, lane M M tt ). Time- 
courses using newly synthesized RNA reproducibly resulted in 
3 5 endpoints slightly less than 50%. Phosphorothioate linkages 
slowly convert to phosphates during storage and this factor 
may explain why slightly more than 50% of the substrate is 
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cleaved (Herschlag, 1991) . Time-courses using this 
oligonucleotide after it was first synthesized were cleaved 
to <48%. Including 10 mM DTT in the reaction had no effect 
on the extent of cleavage. 

5 

Figure 7: Determining the stereochemical identity of the 
less reactive substrate. Three different synthetic 
substrates (0.4 nM) were digested with 3 U//il of snake venom 
phosphodiesterase (SVPD, Boehringer- Mannheim) , in a reaction 

10 containing O.lMg/^1 carrier tRNA, 50 mM Tris-HCl pH 8.5 and 
10 mM MgCl 2 . Incubations were at 37°C for 0, 1, 2, 5, 10 and 
30 minutes, respectively. The lanes marked "phosphate" 
represent digestion of a normal all -phosphate substrate 
(described in Pig. 5) . The lanes marked "Rp + Sp" are 

15 digestions of a synthetic substrate containing a single 
phosphorothioate at the cleavage site (described in Fig. 6) . 
The lanes marked "Sp" are digestions of the synthetic 
phosphorothioate substrate that was less reactive in the 
presence of the ribozyme. To isolate this substrate, the 

20 racemate (Rp+Sp) was cleaved to apparent completion by the 
ribozyme {as in Fig. 9) , resulting in the 50% cleavage shown 
in lane 'M n . The less reactive band was re-isolated by PAGE 
and purified in the same manner as the other two substrates 
before treatment with SVPD. Radioactive counts were 

25 quantitated and the amount of stalling at +1G was expressed 
as a fraction of the whole lane. This stall product 
fraction was integrated over the timecourse of the reaction 
to give relative amounts of exonucleolytic stalling between 
the three time courses. The two reactions (Rp+Sp, and Sp) 

3 0 that had not reached full digestion by 30 min. were 
extrapolated to theoretical endpoints using the slope of the 
last two time points. Both substrates gave extrapolated 
endpoints of 42 min. The final ratio of stall product at 
+1G was {0.05 : 0.59 : l) for {phosphate : Rp+Sp : Sp}, 

35 respectively. 



Figure 8A and 8B: Nuclease mapping of transcribed 
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substrates and their cleavage products. m Fig. 8 A: 
nuclease digestions of S'-»P end-labelled all-phosphate 
substrate described in the legend to Fig. 4. in Fig. 8B, 
the 5'-"P end-labelled substrate was transcribed from the 
5 same DNA template, but or-thio (Sp) GTP (Amersham) was used 
as the sole source of GTP in the reaction. This substrate 
contains an Rp phosphorothioate residue at every G, 
including the cleavage site. Enzymatic digestions were 
carried out under standard conditions (Donis-Keller, 1977) . 
10 For both A and B, Lane 1 is the full-length substrate. Lane 
2 and 3 are digestions with hydroxide. Lane 4: Digestion 
with endonuclease U2 (BRL, cuts ApN) . Lane 5: enzyme B. 
Cereus (Pharmacia, cuts UpN or CpK) , Lane 6: enzyme CL3 
(BRL, cuts CpN), Lane 7: enzyme T2 (Pharmacia, cuts GpN) , 
15 Lane 8: Ribozyme reaction product (19-mer RNAs) , Lane 9: 
PI nuclease (Boehringer Mannheim, no sequence specificity) , 
Lane 10: Hydroxide ladder on substrate, Lane lis Pi 
digestion of ribozyme reaction product (all-phosphate 
substrate only) . Sequence specific nucleases cleave RNA to 
20 leave terminal 3' phosphate, while alkaline hydrolysis 
leaves predominantly 2' -3' -cyclic phosphates. Pi nuclease 
cleaves RNA to leave 3 ' -OH . 

Figure 9A and 9B: (Fig. 9A) . Quantitative determination of 
25 (Rp) phosphorothioate cleavage by the ribozyme. 
Phosphorothioate substrate was made by transcription of DNA 
template as in figure 8B, except labelling was accomplished 
by including a 3J P CTP in the transcription reaction. 
Nuclease digestion of precursor substrates: RNA (3000 cpm) 
30 was incubated with 0.01 U/fil Pi nuclease in 25 mM acetate 
buffer pH 5.3, 10 mM DTT, and 0.05 /ig/pl tRNA competitor for 
30 min. at 37°c. Lane l. Digestion of all-phosphate 
substrate with PI nuclease. Lane 2= All-phosphate 
substrate digested with ribozyme, followed by PI. Lane 3: 
35 Digestion of thio-substrate with PI. Lane 4: Digestion of 
thio- substrate with ribozyme, followed by Pi. Lanes 5 and 
6: Digestion of thio-substrate with Dl and D5 alone. 
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respectively, followed by PI. All ribozyme incubations were 
performed to completion under the conditions described in 
Figure 6. following ribozyme digestion, RNAs were ethanol 
precipitated, resuspended and digested to completion with PI 
5 nuclease (Lanes 1+3: 0.01 U/pl PI nuclease in 25 mM acetate 
buffer pH 5.3, 10 mM DTT and 0.05 ng/fil tRNA competitor for 
30 min. at 37 Q C. Lanes 2, 4, 5 and 6: as in Lanes 1 and 3 
except tRNA was 0.025 pg/iil) . Any radioactivity in the 
middle (*pCp 9 G) band in lane 4 corresponds to counts from the 

10 small fraction of thio- substrate (~ 18%) that was 
uncleavable by the ribozyme. Note that there cannot be a C 
3' to the cleavage site, that is G+l cannot be a C. (Fig. 
9B) . Kinetics of transcribed ribozyme substrates : Cleavage 
of the (Rp) -phosphorothioate substrate (circles, described 

15 in Fig. 8 and 9A) and the transcribed all -phosphate 
substrate (squares) follow pseudo- first order kinetics for 
more than three half -lives. Reaction conditions and 
saturating RNA concentrations were identical to Figure 4, 
except [S] concentration was 0.1 nM. Values for k obs are 

20 0.00094 and 0.032 min* 1 , respectively (average of two 
determinations) . For the thio-substrate, data is corrected 
for an empirically determined end-point (at 7 half-times) of 
83%. The all-phosphate substrate data is corrected for an 
empirically determined endpoint of 96% (at >7 half-times) . 

25 

Figure 10A-10B: Secondary structures of domains used in 
investigating the first step of group II intron splicing. 
(A) Schematic secondary structure of a group II intron and 
the two transesterif ications involved in self -splicing. 

30 During the first step of splicing (labeled step 1), a 2'- 
hydroxyl group on the indicated adenosine of D6 attacks the 
5 '-splice site (shown as an open circle) . During the second 
step of splicing (labeled step 2), the free 3'-hydroxyl of 
the 5'-exon attacks the 3 '-splice site (shown as a closed 

35 circle) . Regions of Domain 1 (BBS1 and 2) that base-pair 
with the 5'-exon (IBS1 and 2) are shown as wide lines 
(Jacquier and Michel, 1987; Jacquier and Michel, 1990). 
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Parts of the RNA that were not present in this study are 
shown as dotted lines, while RNA that was present is shown 
as a solid line, including 5' -exon/Domains 1-3 (exD123, from 
the 5' -end of the RNA to the double bars between Domains 3 
5 and 4), and Domain 5-6 (D56, the region between double bars 
flanking Domains S and 6). This schematic does not attempt 
to illustrate the exact secondary structure of individual 
domains, which are much more complex than indicated here 
(Michel et al., 1989). (B) The nucleotide sequence of D56 
10 (91 nts). The line at the 5' -end represents 15 nucleotides 
of vector sequence left over from cloning. The 3 '-end of 
D56 terminates at the last nucleotide of the intron. 3'- 
exon nucleotides have been removed so that the second 
transesterification of splicing cannot occur. 

15 

Figure 11A, 11B and 11C: Competition of hydrolysis and 
transesterification during the first step of splicing. (A) . 
Time courses are shown for reaction of internally- labeled 
"P-exD123 (l nM) with excess D5S or D56A- (3 ftM) at 45 f C 
20 using the standard reaction buffer at two different pH 
values. The indicated bands are precursor (exD123, 1003 
nt), D56/D123 <D56 that is connected to D123 through a 2' -5' 
branch, 801 nt) , D123 (linear intron resulting from 
hydrolytic 5'SS cleavage, 710 nt) , D12 (a secondary product 
25 of hydrolysis, 560 nt) , and 5'-exon (293 nt) . The D56/D123 
band occurs in the presence of D56 but is not found in the 
presence of D56A- (see +D56 lanes) . However, hydrolysis 
products are observed in reactions of both D56 and D56A- 
(all lanes except t = 0) . Hydrolysis reactions dominate 
30 over the branching reaction at long reaction times, 
especially at higher pH (+D56 lanes, pH 7.5). (B and C) 
Detailed time-courses of the reaction shown in (A) , between 
"P-exD123 and D56 at pH 7.5. The evolution of hydrolysis 
product (B) simultaneous with branched product (C) shows 
35 different kinetic behavior. The branch bursts quickly and 
then declines while hydrolysis products arise in pseudo- 
firBt order fashion, fitting an exponential curve. 
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Figure 12A, 12B and 12C. Timecourse of branching 

independent of hydrolysis and characterization of branched 
product. (A) Autoradiogram of "P-D56 (1 nM) converting to 
5 32 P-D56/D123 in the presence of exD123 (2 /iM) under standard 
reaction conditions at pH 7.0. At longer times, bands 
develop below 32 P-D56 and 32 P-D56/D123 . These are degradation 
products of both species that accumulate over time and do 
not affect kinetic analysis of the reaction. (B) Plot of 
10 fraction reacted versus time, quant itated from the data 
shown in (A) . The reaction shows exponential behavior that 
levels off at 48% total reaction after several hours. <C) 
Digestion of 32 P-D56 (left) and "P-D56/D123 (right) using PI 
nuclease, which cannot cleave branched trinucleotides. The 
15 "P-D56 is internally labeled with a- 32 P-UTP, so PI digestion 
will result in a radioactive trinucleotide only if branching 
occurs specifically to produce pA(2' -pG> 3 ' -pU (Reilly et 
al., 1989). The trinucleotide was resolved from 

mononucleotide on a one-dimensional 20% polyacrylamide gel 
20 (Reilly et al . , 1989; Moore and Sharp, 1993). As expected, 
only Pi treatment of the branched species (right) resulted 
in 32 P- labeled pA(2 ' -pG) 3 ' -pU, while PI digestion of "P-D56 
results only in mononucleotide (left) . RNAs were reacted 
with PI for 12.5 minutes in the experiment shown, although 
25 the same results were obtained from 5-15 minutes. The ratio 
of branched trinucleotide to mononucleotide was 1:54, 
whereas the expected ratio was 1:21 given the uridine 
composition of D56. This higher ratio may be due to 
hydrolytic instability of the branch during workup or 
30 occasional branching from the wrong nucleotide. The latter 
explanation is unlikely given that D56A- does not branch 
(Figure 11A) . 

Figure 13A,13B: Single -turnover kinetic analysis of the 
35 forward branching reaction results in a K« curve for D56. 
Figure 13 A. Reactions such as those shown in Figure 12A were 
performed with increasing concentrations of exDl23. The 
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resultant k ob . values (determined as shown in inset) were 
plotted as a function of [exD123] and an apparent binding 
curve is generated, reflecting the affinity of D56 for 
exD123 . Using a graphing program, the data were fit to a 
5 standard equation describing 1:1 bimolecular association, 
resulting in a K,, of 510 ± 160 nM and J^ t of 0.041 ± 0.004 
min' 1 (Pyle and Green, 1994) . Standard errors shown were 
obtained from these fits as described previously. All 
semilog plots used for determination of k oba values were 
10 corrected for a 51% endpoint determined empirically after 7 
half-times. Variance in individual k^. values, the points 
on this curve, is described below. (Fig 13B) : Semilog plot 
of 32 P-D56 (1 nM) reaction with exD123 (2 jiM) , having 
corrected the reaction for a 48% endpoint (shown in Figure 
15 12B) . The observed rate of reaction is obtained by solving 
the equation of the line for the reaction half-time, as 
described in the experimental section. The k ob8 for the inset 
plot shown is 0.033 min* 1 . A small burst of initial rate (at 
4% reaction) can be seen by the fact that a line 
20 reproducibly crosses the y-axis below zero. Individual k ob . 
values obtained on different days, with different RNA 
stocks, were found to have a variance of ±37% with 95% 
confidence based on six measurements. This is a very 
conservative estimate of maximum error given that the 
25 variance in timecourses performed . on the same day is 
observed to be much lower. For example, the D56 and D56A- 
hydrolysis reactions, studied side by side on the same day, 
both yielded k obs values of 0.018 min' 1 (see text). 

30 Figure 14A, 14B, 14C : Kinetic analysis and 

characterization of the reverse -branching reaction: (A) 
Timecourse for the reverse branching reaction: Upon 
addition of saturating 5'-exon, the "P-D56/D123 branched 
species (1 nM) converts to "P-D56 under standard reaction 

35 conditions at pH 7.0 (lanes 1-10, reaction timepoints of 
0.75, 1.0, 1.5, 3.5, 7.5, 10, 15, 30, €0 and 120 minutes). 
In the absence of 5'-exon, 32P-D56 is liberated very slowly. 



r 
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as a product of hydrolytic attack on the 2' -5' link (lanes 
11-18, reaction timepoints of 0, 2.5, 5, 10, 15, 30, 60, and 
120 minutes) . By quant itating the amount of random 

breakage in the branched molecule and comparing it to 
5 hydrolysis of the 2' -5' linkage, we determined that the 2'- 
5' -linkage is approximately 130-fold more labile than 
surrounding 3' -5' linkages in this context (lane 18). (B) . 
Covalent attachment of 5'-exon to D123 in the reverse 
branching reaction: The 5'-exon (2 nM) was 32 P-end labeled 

10 and reacted with D56/D123 (0.1 nM, contained a trace of 
radioactivity for purposes of isolation, lanes 3-5) and 32 P- 
D56/D123 (0.4 nM, lanes 6-8) under standard reaction 
conditions at pH 7.0. Timepoints were 30 (lanes 3 and 6), 
60 (lanes 4 and 7) and 120 minutes (lanes 5 and 8) . Lane 1 

15 is "P-D56; lane 2 is 32 P-5'-exon (-17/0); and lane 9 is 32 P- 
D56/D123. In lanes 3-5, 32 P-5'-exon reacts and shifts up to 
form band "P-ex/D123 . In lanes 6-8, the same shift is 
observed, concomitant with generation of equivalent molar 
amounts of free "P-D56 (there are two- fold more counts in 

20 the 32 P-ex/D123 band because it is end-labeled while D56 is 
body- labeled) . (C) Kinetic analysis of the reverse 
reaction: The data shown in (A) were plotted, including 
additional early reaction timepoints. The fraction reacted 
vs. time fit a biexponential equation representing a 

25 combination of two rates (experimental section) . The first 
phase has a rate of 0.95 ± 0.17 rain" 1 . A total of 19 ± 1.3% 
of the precursor reacts through this phase. The second 
phase proceeds with a rate of 0.011 ± 0.0010 min* 1 . The 
variances shown are the standard errors calculated from fit 

30 to the equation described in Experimental Details below. 

Figure 15: Overall reaction mechanism and kinetic 
parameters for the first step of group II intron self- 
splicing. D56 binds to exD123 with a Kj of 0.51 ± 0.16 jiM 
35 (Figure 13) . There are three reactions that begin to 
compete after complex formation occurs: Whether or not it 
is in the proper conformation, the D56/exD123 complex can 
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catalyze hydrolysiB at the 5' -splice site with a chemically 
limited rate of 0.018 ± 0.0070 miir 1 (k^, Figure 11 and 
text) . At the same time, the D56/exD123 complex samples 
conformational states, eventually reaching one that is 
5 capable of branching (schematized as a bend in D56) . This 
process occurs with a rate of 0.041 ± 0.0040 min' 1 (k f , confl , 
Figure 13). At this point, the 5' -exon can dissociate (K d 
~ 10 nM) or it can attack the branch site. The branched 
species is more conformationally constrained, and 19% of it 
10 reacts rapidly with 5' -exon at a rate of 0.95± 0.17 min 1 
(kriburso' Figure 14C) . The remaining molecules sample an 
active conformation and react with a rate of 0.011 ± 0.001 
min" 1 (k r(conf)f Figure 14C) . 

15 Figure 16A-16B: Catalytic mechanism of RNA cleavage by the 
minimal group II ribozyme and schematic representation of 
substrate recognition and cleavage. Figure 16A. The 
catalytic mechanism of RNA cleavage is shown; water (or 
hydroxide) is the nucleophile. Figure 16B. The minimal group 

20 II ribozyme comprises domain l (Dl f nts«=nucleotides) of the 
six domains present in the full-length intron and contains 
two exon-binding sites (EBS1 and 2) . The EBSs hybridize to 
intron binding sites (IBS1 and 2) on the oligonucleotide 
substrate, separated by a three nucleotide link. In the 

25 presence of a second RNA moleucle containing domain 5 <D5) 
of the group II intron, the substrate is cleaved at the site 
indicated by the arrow. 

Figure 17: Determination of pseudo-first order rate 
30 constants for the cleavage of chimeric substrates at 
ribozyme saturation. Rates were measured simultaneously, 
under standard reaction conditions where the concentration 
of Dl was 100 nM, the concentration of D5 was 3 pM 
(experimental section) . Rates of reaction (k^,) were 
35 determined by correcting each data set for the amount of 
total inactive substrate and plotting the data as shown, 
where p = fraction product at each timepoint. By solving 
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the equation for each line, the half-time of reaction and 
subsequently the k^ were determined as described previously 
[Michels and Pyle, 1995] . For the the rapidly reacting 
substrates shown, the rates were 0.021 rain* 1 for rS, 0*018 
5 min' 1 for -3dU, 0.026 rain* 1 for -4dU and -5dU, 0.025 rain* 1 for 
-6dA and 0.020 min* 1 for -7dC. The variance in these rates 
was ±16% determined using eight separate trials. For -ldC 
the reaction rate was 0.0013 rain" 1 , with a variance of ±23% 
determined using eight different trials. For -2dU, the 

10 rate was 0.0042 rain" 1 , with a variance of ± 29V based on 
three trials. All errors were calculated using a 95% 
confidence limit and trials were performed on different 
days. For experiments performed side-by-side as shown here, 
relative error is expected to be much lower than the 

15 variances conservatively stated above. 

Figure 18: Mapping the position and termini of substrate 
cleavage products: High-resolution mapping with 
endonucleases and alkali: Lanes 1 to 7 represent mapping 

20 of 32 P 5' -end labeled rS, lanes 8 to 14 are maps of 32 P 5' -end 
labeled -ldC. For each substrate, lanes 1 and 14 are full- 
length substrate; lanes 2 and 13 are partial digestions of 
substrates by Tl nuclease {Pharmacia, cuts GpN) ; lanes 3 and 
12, are partial digestions with enzyme S. Cereus 

25 (Pharmacia, cuts UpN or CpN) ; lanes 4 and 11 are partial 
hydroxide digestions of the full length substrates; lanes 5 
and 10 are partial digests on the 50/50 mix of precursor and 
product from a ribozyme reaction; lanes 6 and 9 are the 
50/50 mix of product and precursor from a partial ribozyme 

30 reaction; lanes 7 and 8 are partial digestions with nuclease 
PI (Boehringer Mannheim, cuts DNA or RNA) . The mobilities 
for substrate (S) and product (P) are marked on the side of 
the gel. Bands in lanes 6, 7 r 8 and 9 have a 3'-hydroxyl 
terminus, causing them to migrate more slowly than their 

35 corresponding bands in other lanes, which possess a 3'- 
phosphate or 2 '-3 '-cyclic phosphate terminus. There are 
also mobility differences due to the presence of a single 
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deoxynucleotide, which tends to make chimeric 
oligonucleotides migrate more rapidly (lanes 6-9; lanes 7 
and 8) . Verification of deoxynucleotide cleavage by 
characterizing the 3' termini of reaction products. a 
5 50/50 mixture of oligonucleotide substrate and product from 
a partial D1/D5 ribozyme reaction was subjected to oxidation 
(+ NalOJ and then 3 '-end labeled with a hydrazine dye 
derivative (+ Lucifer Yellow). Lanes 1 to 4 are rS 
substrate, and lanes 5 to 8 are -ldC substrate. The 
10 identity of each band is indicated to the right. After 
oxidation and covalent dye-labeling, both the precursor and 
product band of the rS substrate undergo a large mobility 
shift (lane 4). For -ldC, only the substrate band shifts to 
slower mobility (lane 5), indicating that the cleavage 
15 product of -ldC has a 3' deoxy terminus, which cannot be 
oxidized with NalO,. Doublets in lanes 4 and 5 are most 
likely due to the products of Lucifer Yellow hydrazide 
attack on one or both of the aldehydes produced by oxidative 
attack of the ribose sugar (Hannske et al . , 1974]. 
20 Following oxidation with NalO,, there is a small mobility 
shift for species with a 3 '-ribose because the sugar is 
oxidized to a larger dialdehyde moiety (compare lanes l and 
2). Additional bands in the +NaIO« lanes are attributed to 
periodate attack on residues other than the ribose sugar 
25 [Rizzino and Freundlich, 1975] . Verification of 
deoxynucleotide cleavage by the D123 ribozyme. This 
experiment and figure is the same as in (B) , except that 
D123 was the ribozyme instead of Dl. As in (B) , the -ldC 
cleavage product is shown to have a 3' -deoxynucleotide 
30 terminus by inspection of lane 5. 

Figure 19A, 19B and 19C: Determination of kinetic parameters 
for rS, -ldC and -2dU substrates: By plotting k oh , as a 
function of [Dl] at saturating [D5] , one can fit the data to 
35 a binding curve and obtain kinetic parameters for a reaction 
that follows Michaelis-Menten kinetics. Data were fit as 
previously described [Michels and Pyle, 1995; and Pyle and 
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Green, 1994] . (A) For the rS substrate, « 4.9 nM and k cat 
« 0.025 rain" 1 . (B) For -ldC, K„ ■ 14 nM and )^ at - 0.0019 
rain" 1 . <C) For -2dU, K„ = 10 nM and k c . c « 0.0051 min' 1 . All 
k cat values determined from these fits are in good agreement 
5 with the saturating k ob8 values reported in Fig. 1. Errors 
and k c . t /K B values determined from this experiment are 
provided in Table 2. 

Figure 20: Determination of saturating k^ values for 

10 oligonucleotide cleavage by the fast D123 ribozyme. 
Cleavage of rS and -ldC substrates by D123 (100 nM) and D5 
(3 /iM) . The apparent k^ values from this data are 0.88 ± 
0.089 min* 1 for cleavage of rS and 0.17 ± 0.021 min -1 for - 
ldC. Errors were calculated with 95V confidence from 3 

15 trials, using average k oba values of 0.84 min" 1 and 0.16 min" 1 
for rS and -ldC, respectively. Because D123 cleaves 
oligonucleotides so efficiently, rate-limiting 
conformational processes become evident in the seconds 
timescale. A small initial lag exists at timepoints earlier 

20 than 15 seconds for both substrates. By conducting 
fluorescence binding experiments and altering the order of 
component addition, the lag was found to derive from a 
folding event related to the association of D5 (Z. F. Qin, 
data not shown) . After this initial lag, the data are almost 

25 perfectly linear such that, at longer reaction times, the 
folding event is not rate -limiting. Rates were calculated 
from the fit to the linear portion of this data (after 15 
seconds) . 

30 Figure 21: Summary of deoxynucleotide effects on the 
chemical rate of D1/D5 ribozyme cleavage: The magnitude of 
k cbctB diminution is plotted for the single and block 
deoxynucleotide substitutions along the oligonucleotide 
substrate strand. The k^ effect for IBS2 is labelled as WT 

35 for block because the k cat for cleavage of the d(IBS2) 
oligonucleotide was the same as the rS substrate, provided 
that additional Dl was added to overcome slight weakness in 
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the chimeric EBS2/IBS2 helix. 



Figure 22A, 22B and 22C: The three reactions observed during 
in-vitro self -splicing of group II intron ai5g. The 

5 reaction products shown at right are the species observed on 
high-resolution polyacrylamide gels of self -splicing 
reaction products. The second step of splicing is extremely 
rapid and the first step of splicing is rate limiting. 
Therefore, although splicing intermediates are shown in (A) 

0 and <B) , these are barely detectable, if they are observed 
at all. SER(C) can take place immediately after splicing, 
so its products (free exons} are also evident on separation 
gels. 



15 Figure 23A, 23B and 23C: High- resolution sequencing gel of 
reaction products from self -splicing. The product band are 
labeled as follows: A is lariat intron; B is precursor RNA; 
C is linear intron (which comigrates with broken lariat) ; D 
is spliced exons; E is free 3'-exon; F is free 5'exon. In 
20 all cases, linger timecourses than those shown were required 
for kinetic analysis (A) Splicing in 100 mM MgCl 2 and 500 mM 
(NH t ),S0«: Lanes 1-15 correspond to timepoints of 
0,1,2,3,4,8,10,15,20,30,40,60,120,240, and 900 minutes,, 
respectively. (B) Splicing in 100 mM MgCl 2 and 500 mM KC1: 
25 Lanes 1-6 correspond to timepoints of 0,10,15,20,30, and 40 
minutes, respectively. (C) Splicing in 100 mM MgCl 2 and 500 
mM NH 4 C1: Lanes 1-5 correspond to timepoints of 12,20,30,45. 
and 60 minutes, respectively. Splicing products D-F are 
also observed under the latter two sets of conditions. 
30 although the lower half of the gels are omitted. These 
shorter products were identical to those shown in full lanes 
of Figure 8. 



Figure 24A, 24B, 24C and 24D: Timecoureses for self -splicing 
under low-salt conditions. Experimental data for the 
decrease in precursor RNA (□) , and the evolution of linear 
(•) and lariat (▲) intron products are plotted and fit 
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using equations 1-3 in the text. Both reactions proceed to 
completion at very long times, but the significant amount of 
RNA degradation that occurs after five hours at pH 7.5 
precludes quantitative analysis after that time. (A) 
5 Reaction profile in 100 mM MgCl a . The k, value determined 
from decrease in precursor was 0 . 0049±0 . 00022 min* 1 . The 
coefficients f hy and f br , indicating the relative amounts of 
linear and lariat intron were 0.75±0.11 and 0.18±0.0060, 
respectively. These represent the horizontal asymptotes of 

10 the product curves shown on the plot. The k T estimations 
determined from the fit of linear and lariat evolution were 
0.0047 min 1 and 0.0096 min 1 , respectively. Because they are 
based on a relatively smaller number or radioactive counts, 
these values are less accurate than k T determined from 

15 decrease in precursor. (B) Reaction profile in 10 mM MgCl 2 
and 5 nM protamine. The k T value determined from decrease 
in precursor was 0 . 0029±0. 00014 min* 1 and the coefficients f hy 
and f br were 0.77±0.088 and 0.14±0.020, respectively. The k T 
estimations determined from the fit of linear and lariat 

20 evolution were 0.0030 min 1 and 0.0056 min" 1 , respectively. (C 
and D) Kinetic parameters of ribozymes at low salt 
conditions. 

Figure 25A-1, 25A-2, 25B, 25C; Semi-log plots of reaction 
25 showing biphasic kinetic profiles. One of the conventional 
ways to plot first-order kinetic data is a semi-log plot. 
The reaction timecourses for self -splicing show biphasic 
behavior for which each phase is linear. The inflection 
point between the two phases of reaction occurs at 
30 approximately 50% (-0.3 on the y-axis) , regardless of what 
type of monovalent in is added. All of the reactions shown 
were conducted in 100 mM MgCl 2 and 500 mM of the monovalent 
slats indicated. The inset to the pot of (NH 4 ) 2 S0 4 shows a 
considerable burst of activity during the first 13% of 
35 reaction (-0.06 on the y-axis). Note that the data for the 
next phase of splicing (after the burst) is completely 
linear when shown in detail. 
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Figure 26A, 26B f 26C; Kinetic analysis of self -splicing in 
the presence or monovalent cations: The experimental data 
for decrease in precursor (□) , and the evolution of linear 
(•) and lariat <▲) intron products are plotted and fit 
5 using equations 1-10 of the text. The kinetic parameters 
cited below and in Table 5 were determined from these fits. 
(A) Kinetics of splicing in 500 mM Kcl and 100 mM MgCl 2 : From 
the fit of the data, the k\ and k B T values for the reaction 
of precursor population A and B were 0.028±0.0031 min' 1 and 
10 0.0 042±0. 00026 min" 1 , respectively. The fractional 

coefficients for A(f\ y and A £ r ) were 0.32±0.040 and 
0.093±G.022, respectively. The coefficients for B(f B hy and 
f B br ) were 0.48±0.068 and 0.16±0.073 respectively. These 
were used to compute the individual rate constants for 
15 branching and hydrolysis that are reported in Table 5. 
Values of k A x determined from linear and lariat intron 
evolution were 0.019 min 1 and 0.037 min'\ respectively. The 
k 8 T determined from lariat evolution was 0.0030 frtin 
Based on other estimations, the k* T value for linear 
20 evolution was fixed at 0.0035 min' 1 . (B) Kinetics of 
splicing in 500 mM NH 4 C1 and 100 mM MgC^ : From the fit of 
the precursor data, the k\ and k D T values for reaction of 
precursor were 0.038±0.0044 min' 1 and 0 . 0039±0 . 00027 min' 1 , 
respectively. From the fit of product evolution data, the 
25 fractional coefficients for A(f A hy and f\ r ) were 0.38±0.044 
and zero, respectively. The coefficients for B(f B hy and f 8 te ) 
were 0.10±0.018 and 0.51, respectively. The latter value 
was fixed based on the linear population at long timepoints. 
Under this reaction condition, early timepoints for linear 
3 0 evolution could not be accurately differentiated from 
background counts so they were not included. These values 
were used to compute the individual rates of branching and 
hydrolysis shown in Table 5. The value of k\ determined 
from lariat formation was 0.037 min - 1 . Values of B k T 
35 determined from branching and hydrolysis were 0.0057 min' 1 
and 0.0032 min' 1 , respectively. (C) Kinetics of splicing in 
500 mM <NH«) 2 S0 4 and 100 mM MgCl 2 : From the fit of precursor 
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reaction using a fixed 13% burst population, the rates of 
the burst, k\ and k* T were 0.22±0.037 min* 1 , 0.047±0.0033 min 
- 1 and 0.0028±0. 00014 min * l , respectively. The fractional 
coefficients for A(£\ y and f\ r ) were 0.39±0.0044 and zero, 
5 respectively. The coefficients for B(f B hy and ^ ) were 
0.22±0.050 and 0,39, respectively. The latter value was 
fixed based on the linear population observed at long 
timepoints. These values were used to compute the 
individual rates of branching and hydrolysis shown above and 
10 in Table 5. 

Figure 27: Comparison of ai5g self -splicing products and 
reactions: In order to qualitatively illustrate the relative 
amounts of reaction products and to show the differential 

15 amounts of nonspecific RNA degradation under all the 
conditions described here and in the literature, short 
timecourses were performed and subjected to electrophoresis 
on a 4% denaturing high-resolution polyacryl amide gel. The 
wavy mobility of RNA in lanes 14-18 is due to stretching of 

2 0 the low -percent age gel during work-up and is not a 
characteristic of the indicated reactions. As in Figure 
23A, the bands are assigned as follows: A, lariat intron; B, 
precursor; C, linear intron; D, spliced exons; E, free 3'- 
exon; F, free 5'-exon. Timepoints were chosen in order to 

25 show similar amounts of reaction under different conditions. 
Lane 1 is the precursor RNA at t=0. For lanes 2-20, the 
reaction time are 5, 20, 65, 900, 1.5, 20, 45, 6, 14, 50, 
30, 70, 160, 35, 75, 170, 30, 70 and 160 minutes, 
respectively, under the reaction conditions shown. Lane M 

30 contains RNA markers of the molecular weights shown. For 
reactions shown in lanes 2-11, the reaction conditions 
included 100 mM MgCl 2 , together with 500 mM of the indicated 
salt. Lanes 12-14 are self-explanatory. Lanes 15-17 
include 10 mM MgCl a , together with 2 mM spermidine. Lanes 

35 18-20 include 10 mM MgCl 2 , together with 5 nM protamine. 
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DBTAILED DBSCRI PTTON OF TOE Ttiy^Tn^ 

This invention is directed to a composition for catalyzed 
oligonucleotide cleavage comprising a synthetic non- 
naturally occurring oligonucleotide compound. The compound 
5 comprises nucleotides whose sequence defines a conserved 
group II intron catalytic region and nucleotides whose 
sequence is capable of hybridizing with a predetermined 
oligonucleotide target sequence to be cleaved, such target 
sequence not being present within the compound. The 

10 composition also includes an appropriate oligonucleotide co- 
. factor which may be provided in cis or trans. Preferably, 
the conserved group II intron catalytic region is a group II 
intron domain I catalytic region. In one embodiment the 
conserved group II intron domain I catalytic region may 

15 further comprise a conserved portion of a group II intron 
domain II, a group II intron domain III, a group II intron 
domain IV, a group II intron domain V, or a group II intron 
domain VI. 



20 In another embodiment the conserved group II intron domain 

I catalytic region may further comprise a portion of a group 

II intron sufficient to maintain catalytic activity. The 
conserved portions of the domain I region may comprise 
structural features or conserved secondary structural 

25 features sufficient for such activity. As used herein, 
"conserved structural features" include the EBS1 domain, 
EBS2 domain, al domain, a' domain, € domain, e' domain, and 
the bulge between the D domain and the D' domain. The 
sequenes of the EBS1 and EBS2 domains may be modified to 

3 0 hybridize with the desired target sequence, so long as 
sufficient secondary structure is preserved. 

In the composition above, the nucleotides whose sequence is 
capable of hybridizing with a predetermined oligonucleotide 
35 target sequence to be cleaved may be two hybridizing 
regions, each region having 2 to 12 nucleotides capable of 
hybridizing with the oligonucleotide target to be cleaved. 
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In one embodiment, each of the regions has € or 7 
nucleotides. Preferably, the synthetic non-naturally 
occurring oligonucleotide compound comprises between 80 and 
1000 nucleotides. In one embodiment the compound comprises 
5 300-600 nucleotides. 

In one embodiment of the invention, the appropriate co- 
factor is an oligonucleotide compound from a conserved 
portion of a group II intron domain V having the formula: 

10 3' (X) n (X) n . 5' 

X * X 
X t X 
X * G 
G * C 

15 (X) c * (X) c , 

G ♦ U (X) d . X 

(X) a * A 
(X) d A X 

20 wherein each X represents a ribonucleotide which may be the 
same or different; wherein each of (X) n , (X) n ., (X) c , (X) c . , 
(X) dw (X) a and <X) d represents an oligonucleotide; wherein 
n, n' c, c' , d, d' , and a each represents an integer which 
defines the number of nucleotides in the oligonucleotide 

2 5 with the provisos that n and n' are greater than or equal to 
1; a represents an integer which is greater than or equal to 
1; d and d' represent an integer which is greater than or 
equal to 5; c and C represents an integer which is greater 
than or equal to 4; wherein each * represents base pairing 

30 between the nucleotide located on either side thereof; 
wherein each t may or may not represent base pairing between 
the nucleotide located on either side thereof; and wherein 
each solid line represents a chemical linkage providing 
covalent bonds between the nucleotide located on either side 

35 thereof. In one embodiment, 5' (X) n .XX may be 5' (X) n .XA. 

The compositions above may further comprise a divalent 
cation such as Mg 2+ or Mn 2 *. Preferably, the predetermined 
oligonucleotide target sequence to be cleaved and/or ligated 
4 0 is mRNA e.g. a mammalian mRNA, a yeast mRNA, a bacterial 
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mRNA or a viral mRNA. However, the predetermined 
oligonucleotide target sequence to be cleaved may be DNA. 
The mRNA may encode a growth factor such as an angiogenic 
factor, a basic fibroblast growth factor, a colony- 
5 stimulating factor 1, cystic fibrosis transmembrane 
conductance regulator, an epidermal growth factor, an 
erythropoietin, a fibroblast growth factor, a G-protein, a 
granulocyte -macrophage colony stimulating factor, a growth 
hormone, IL-1, IL-1R, IL-2, IL-2R, IL-4, IL-6, an insulin- 

10 like growth factor, an insulin-like growth factor 1, an 
interferon, an interleukin, a keratinocyte growth factor, 
luteinizing hormone receptor, MDR-1, a nerve growth factor, 
a platelet derived growth factor, a scatter factor, a 
transforming growth factor a, a transforming growth factor 

15 p, a transforming growth factor, or a tumor necrosis factor. 

The mRNA may encode an oncogene or a tumor suppressor gene 
such as bcl-2, Jbcr-aJbl, Jbe*, BPV, c-aJbi, c-fes, c-fms, c- 
fos, c-H-ras, c-^it, c-myb, c-myc, c-mos, c-oea, cerbB, DCC, 

20 erbA, erbB-2, ets, fig, FSFV gp55, Ha-ras, HIV tat, HTLV-l 
tat, JCV early, jun, L-myc, Ick, LPV early, met, N-jnyc, NF- 
1, N-ras f neu, p53, Py mTag, pim-1, ras, RB, rel, 
retinoblastoma-1, SV-40 Tag, TGF-a, TGF-^, trk, trkB, v-ajbi, 
v-H-ras, v-jun, or WT-1. The mRNA may be associated with 

25 a chromosomal translocation or have one or more point 
mutations. In addition, the mRNA may be an mRNA whose 
overproduction is associated with a disease or condition. 

The mRNA may be a viral mRNA associated with a 
3 0 cytomegalovirus, an Epstein-Barr virus, a hepatitis B virus, 
a hepatitis C virus, a herpes simplex virus, a herpesvirus, 
a HIV-l virus, an immunodeficiency virus, an influenza 
virus, a papillomavirus, a picomavirus, a polio virus or a 
T-cell leukemia virus. Alternatively, the mRNA may be a 
35 plant mRNA associated with alfalfa, apples, asparagus, 
bananas, broccoli, carrots, celery, chicory, coffee, 
cabbage, mustard, corn, cottonseed, squash, cucumber, 



WO 96/22689 



PCT/US96/01337 



-30- 

cantaloupe, grapes, lettuce, palm, potato, rapeseed, 
raspberry, soybean, sunflower, strawberry, tomato, or wheat. 

The invention is also directed to an oligonucleotide 
5 transfer vector containing a nucleotide sequence which on 
transcription gives rise to the synthetic non-naturally 
occurring compounds above. The transfer vector may be a 
bacterial plasmid, a bacteriophage DNA, a cosmid, or an 
eukaryotic viral DNA. 

10 

The invention is also directed to host cells transformed by 
the transfer vector above. The host cells may be 
prokaryotic host cells or an eukaryotic host cells such as 
E. coli, a monkey COS host cell, a Chinese hamster ovary 
15 host cell, a mammalian host cell or a plant host cell. 

The invention is also directed to a method for cleaving an 
oligonucleotide comprising the steps of: a) providing a 
synthetic non-naturally occurring oligonucleotide compound 
20 or compounds as defined above and b) contacting the 
synthetic non-naturally occurring oligonucleotide 
compound (s) with a target oligonucleotide to cause cleavage 
of the target oligonucleotide. 

25 A transgenic, non-human vertebrate animal, having one or 
more cells bearing a DNA sequence encoding the composition 
above. The transgenic animal may be a mammal such as a cow, 
goat, sheep, pig, horse, dog, cat, or rodent. In addition, 
the animal may be a fish, a chicken, or a turkey. 

30 

The invention is also directed to a method of treatment 
which comprises: a) recovering a subject's appropriate 
cells; b) transforming the cells with a DNA sequence 
encoding the composition of above; and c) introducing the 
35 resulting transformed cells into the subject so as to treat 
the subject. 
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Synthetic preparations of mRNA are well known (see Sambrook 
et al., 1989, Molecular Cloning: ft Laborat-n^y M rTri1 , 2 nd 
Ed. ) . Mixed DNA-RNA oligomers with modified base pairs for 
the ribozyme or minizyme can be prepared by commercially 
5 available DNA synthesizers such as those produced by Applied 
Biosystems. For derivatives see (Uhlmann, E. and Peyman, 
A., 1990, Chemical Reviews 90:543-584.), for H-phosphonate 
monomers (Agrawal et al U.S. Patent No. 5,149,798). 

10 The compounds of this invention may be covalently or non- 
covalently associated with affinity agents such as proteins, 
antibodies, steroids, hormones, lipids, specific nucleic 
acid sequences, intercalating molecules (such as acridine 
derivatives, for example 9 -amino acridine) or the like to 
15 modify binding affinity for a substrate nucleotide sequence 
or increase affinity for target cells, or localization in 
cellular compartments or the like. For example, the 
compounds of the present invention may be associated with 
RNA binding peptides or proteins which may assist in 
20 bringing the endonuclease into juxtaposition with a target 
nucleic acid such that hybridization and cleavage of the 
target sequence may take place. Nucleotide sequences may be 
incorporated into the 5' and 3' ends of the group II intron 
increase affinity for substrates. Such additional 
25 nucleotide sequences may form triple helices with target 
sequences (Strobel, S.A., et al., 1991, Nature 350:172-174.) 
which may enable interaction with intramolecularly folded 
substrate. Alternatively, modified bases within the 
additional nucleotide sequences may be used that will 
30 associate with either single stranded or duplex DNA 
generating base pair, triplet, or quadruplet, interactions 
with nucleotides in the substrate. 



The compounds of the claimed invention may be further 
35 stabilized using methods in the literature for example the 
use transcription terminators on the 3' end such as the T7 
terminator, p- independent terminator, cry element (Gelfand 
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et al. U.S. Patent No. 4,666,848) or the TrpE terminator. 
Furthermore, sequences such as the poly (A) addition signal 
AATAAA may be added and strategies involving changing the 
length of the 3' non-coding region (see Gillies, U.S. Patent 
5 No. 5,149,635). These techniques can be used to stabilize 
RNA in the compound. 

The invention also embodies methods of production of the 
nucleic acid based compounds described above comprising the 

10 steps of: (a) ligating into a transfer vector comprised of 
DNA, RNA or a combination thereof a nucleotide sequence 
corresponding to said compound; (b) transcribing the 
nucleotide sequence of step (a) with RNA polymerase; and (c) 
recovering the compound. The invention also includes 

15 transfer vectors, bacterial or phage, comprised of RNA or 
DNA or a combination thereof containing a nucleotide 
sequence which on transcription gives rise to the compounds 
described above. 

20 The invention described herein also provides a method of 
cleavage of a specific RNA target sequence which comprises 
reacting the compound with the target sequence so as to 
thereby cleave the specific target sequence. Such target 
sequences may be indigenous to mammals or plants . In one 

25 embodiment, the target sequence is in a viral gene. The 
invention also provides a method for the treatment of viral- 
diseases in plants and animals'. 

Further, many methods have been developed for introducing 
30 cloned eukaryotic DNAs into cultured mammalian cells 
{Sambrook et al., 1989): calcium phosphate- or DEAE-dextran- 
mediated transf ection; polybrene; protoplast fusion; 
elect roporat ion; and direct microinjection into nuclei. The 
present invention, however, extends to other means of 
35 transfer such as genetic bullets {e.g. DNA-coated tungsten 
particles, high-velocity micro projectile bombardment) and 
elect roporat ion amongst others (Maliga, P. (1993) Towards 
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plastid transformation in flowering plants. Tibtech 11:101- 
106; Bryant, J. (1992) Transgenic wheat plants: the end of 
the beginning. Tibtech 10:342-343.; or Shimamoto, K., 
Terada, R., Izawa, T. f and Fujimoto, H. (1989) Fertile 
5 transgenic rice plants regenerated from transformed 
protoplasts. Nature 338:274-276.). 

Further, the compound described herein may be used in plants 
to cleave undesirable mRNA. The appropriate cleavage would 

10 lead to phenotypic changes. Phenotypic changes in plant 
cells or plants may include drought resistance, salinity 
resistance, resistance to fungal, viral or bacterial 
infection; modifications of growth characteristics; 
sterility; fruit production; flowering; senescence; altering 

15 oil seed metabolic pathways to increase production; and the 
like (see Shewmaker et al . U.S. Patent No. 5,107,065). It 
is evident that one or more RNA involved in determining 
phenotype are identified, such RNAs may be inactivated by 
• cleavage utilizing the endonuclease of this invention and 

20 thus the phenotype of the plant or plant cell altered. 
Diseases or infections which may be treated in plants with 
endonucleases of this invention include fungal infection, 
bacterial infections (such as Crown-Gall disease) and 
disease associated with plant viral infection. 

25 

Phenotypic modifications within animals (including in some 
applications man) which may be effected by cleaving and thus 
inactivating target RNAs associated with phenotype would 
include growth characteristics of animals, fertility, 

30 skin/cosmetic modifications, reproductive characteristics, 
disease resistance and the like. A myriad of applications 
arise for phenotypic modifications in animals, and plants as 
previously mentioned. One or more RNAs associated with a 
given endonucleases may be targeted against such RNAs for 

35 their inactivation with consequential phenotypic 
modification. 
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Prokaryotic and eukaryotic cell cultures may be 
phenotypically modified by treatment with endonucleases of 
this invention. For example, bacterial cultures or yeast 
cultures involved in production of food components (such as 
5 cheese, bread and dairy products) and alcoholic beverage 
production may be treated so as to modify enzyme content, 
flavor production, cell growth rate, culture conditions and 
the like. Eukaryotic and prokaryotic cells in culture may, 
for example be protected from infection or disease 
10 associated with mycoplasma infection, phage infection, 
fungal infection and the like. 

The compounds of this invention may also be used to treat 
diseases or infection in humans, animals, plants, or 
15 prokaryotic or eukaryotic cells. The ability to treat 
disease or infection is based on the fact that the compounds 
of this invention are capable of cleaving any RNA which 
contains a suitable cleavage site. Most RNAs will contain 
one or more suitable cleavages sites. 

20 

The period of treatment would depend on the particular 
disease being treated and could be readily determined by a 
physician. Generally treatment would continue until the 
disease being treated was ameliorated. 

25 

Examples of human and animal targets may be bacterial and 
prokaryotic infection, viral infection and neoplastic 
conditions associated with the production of aberrant RNAs 
such as occurs in chronic myeloid leukemia, viruses, growth 
30 factors as described above. 

The virus may be a plant virus such as a tobamovirus, a 
tobravirus, a hordeivirus, a potexvirus, a carlavirus, a 
potyvirus, a closterovirus, a tymovirus, a tombusvirus, a 
35 sobemovirus, or a luteovirus. The plant virus may be a 
potato yellow dwarf virus, a cucumber mosaic virus, a tomato 
spotted wilt virus, a tomato mosaic virus, a potato virus X 



WO 96/22689 



PCTOJS96/01337 



-35- 

(PVX) , a potato virus Y (PVY) , a carnation latent virus, a 
tomato rattle virus, a pea early browning virus, a barley 
stripe mosaic virus, a turnip yellow mosaic virus, a barley 
yellow dwarf virus, a beet yellows virus, a potato leaf roll 
5 virus, a tomato bushy stunt virus, a southern bean mosaic 
virus, a maize chlorotic virus, beet necrotic yellow vein 
virus, or a tobacco necrosis virus. 

The transgenic organisms described herein may be 

10 characterized in that they contain in their genome a 
sequence which gives rise, on transcription, to the 
compounds mentioned above. The transgenic organism may be 
a plant or a non-human mammal. Examples of transgenic 
mammals include Leder et al. U.S. Patent Nos. 4,736,866 and 

15 5,175,383; Krimpenfort et al. U.S. Patent No. 5,175,384; 
Wagner et al. U.S. Patent No. 5,175,385; and U.S. Patent 
Nos. 5,183,949 and 5,347,075. The transgenic plant, 

including fruits, and seeds thereof, may be from alfalfa, 
apple, bean, canola (oilseed rape), cantaloupe, corn, 

20 cotton, courgette, cucumber, melon, papaya, pepper, potato, 
rice, soybean, squash, strawberry, sunflower, sweet pepper, 
tobacco, tomato, or walnut. Also included are the plant 
cells transformed by the above-mentioned transfer vector, as 
well as a prokaryotic or eukaryotic cell, plant or animal, 

25 comprising a nucleotide sequence which is, or on 
transcription gives rise to, the compounds described herein. 

An effective amount of a compound of the present invention 
would generally comprise from about 1 nM to about 1 mM 

3 0 concentration in a dosage form, such as a cream for topical 
application, a sterile injectable composition, or other 
composition for parenteral administration. In respect of 
topical formulations, it is generally preferred that between 
about 50 fiM to about 500 /xM endonuclease be employed. 

35 Compounds comprising nucleotide derivatives, which 
derivatives may involve chemically modified groups, such as 
phosphorothioate or methyl phosphonate derivatives may be 
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active in nanomolar concentrations. Such concentrations may 
also be employed to avoid toxicity. 

Therapeutic strategies involving treatment of disease 
5 employing compounds of this invention are generally the same 
as those involved with antisense approaches, such as 
described in the antisense literature (Chrisley, L.A. (1991) 
Antisense Research and Development, 1:65-113; Stein, C.A. et 
al. Science (1993) 261:1004-1012). Particularly, 
10 concentrations of compounds utilized, methods and modes of 
administration, and formulations involved may be the same as 
those employed for antisense applications. 

An "effective amount" as used herein refers to that amount 
15 which provides a desired effect in a mammal having a given 
condition and administration regimen. Compositions 
comprising effective amounts together with suitable 
diluents, preservatives, solubilizers, emulsifiers, 
adjuvants and/or carriers useful for therapy. Such 
20 compositions are liquids or lyophilized or otherwise dried 
formulations and include diluents of various buffer content 
(e.g., Tris-HCL, acetate phosphate), pH and ionic strength, 
additives such as albumin or gelatin to prevent absorption 
to surfaces, detergents (e.g., Tween 20, Tween 80, Pluronic 
25 F68, bile acid salts), solubilizing agents (e.g., 
Thimerosal, benzyl alcohol), bulking substances or tonicity 
modifiers (e.g., lactose, mannitol) , covalent attachment of 
polymers such as polyethylene glycol to the oligonucleotide, 
complexation with metal ions, or incorporation of the 
30 material into or onto particulate preparations of polymeric 
compounds such as polylactic acid, polyglycolic acid, 
polyvinyl pyrrolidone. etc. or into liposomes, 
microemulsions, micelles, unilamellar or multilamellar 
vesicles, erythrocyte ghosts, or spheroplasts . Such 
35 compositions will influence the physical state, solubility, 
stability, rate of in vivo release, and rate of in Yiyc: 
clearance of the oligonucleotide. Other ingredients 
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optionally may be added such as antioxidants, e.g., ascorbic 
acid; low molecular weight (less than about ten residues) 
polypeptides, i.e., polyarginine or tripeptides; proteins, 
such as serum albumin, gelatin, or immunoglobulins; amino 
5 acids; such as glycine, glutamine acid, aspartic acid, or 
arginine; chelating agents such as EDTA; and sugar alcohols 
such as mannitol or sorbitol. Possible sustained release 
compositions include formulation of lipophilic depots (e.g., 
fatty acids, waxes, oils) . Also comprehended by the 

10 invention are particulate compositions coated with polymers 
(e.g., polyoxamers or polyoxamines) and oligonucleotides 
coupled to antibodies directed against tissue- specific 
receptors, ligands or antigens or coupled to ligands of 
tissue-specific receptors. Further, specific nucleotide 

15 sequences may be added to target the oligonucleotides of 
this invention to the nucleus, cytoplasm or to specific 
types of cells. Other embodiments of the compositions of 
the invention incorporate particulate forms protective 
coatings, protease inhibitors or permeation enhancers for 

20 various routes of administration, including parenteral, 
pulmonary, nasal and oral. 

Suitable topical formulations include gels, creams, 
solutions, emulsions, carbohydrate polymers, biodegradable 

25 matrices thereof; vapors, mists, aerosols, or other 
inhalants. The oligonucleotides may be encapsulated in a 
wafer, wax, film or solid carrier, including chewing gums. 
Permeation enhancers to aid in transport to movement across 
the epithelial layer are also known in the art and include, 

30 but are not limited to, dimethyl sulfoxide and glycols. 

For An-vjtro use, the compounds of this invention are 
generally reacted with a target RNA which contains one or 
more suitable cleavage sites. Optionally, the target RNA 
35 may be purified or substantially purified. The nucleotide 
sequences rhe compounds of this invention are selected so as 
to specifically hybridize or form a double -stranded duplex 
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with a target RNA whereafter cleavage takes place. 
Accordingly, target RNA may be specifically cleaved toitffl 
in the presence of other RNAs which themselves would not be 
cleaved. 

5 

The compounds may be utilized in a manner similar to 
restriction endonucleases, that is for the specific cleavage 
of RNA to facilitate RNA manipulation. All that is required 
for such manipulations is that the target RNA to be cleaved 
10 contains a uracil base and thus a suitable cleavage site. 

The compounds of this invention may be utilized in 
diagnostic procedures, such as the mapping or finger- 
printing of RNA. Specifically, the compounds of this 
15 invention would enable mapping of RNA and may be used to 
detect mutations in RNA sequence. Such procedures may be 
used in research and may also have forensic and other 
diagnostic applications. 

20 RNA cleavage products in-vitro may be readily detected, for 
example, by visualization on acrylamide or agarose gels 
where the amounts of RNA cleaved are sufficiently large for 
direct visualization after separation and reaction with 
nucleotide visualization agents, such as ethidium bromide. 

25 Alternatively, where the target RNA cleaved is present in 
small amounts, such as in a sample containing many RNAs, 
cleavage products may, for example, be detected by using 
radiolabeled probes for sequence complementary to the 
target sequence, or amplification techniques such as PCR 

30 (Sambrook et al., 1989). 

A target RNA for cleavage in-vitro may be derived from any 
source, and may be of animal, viral, bacterial, plant, 
synthetic, or other origin. As RNA is common to all known 
35 living organisms, this invention may be utilized to cleave 
any RNA species having a suitable cleavage site as mentioned 
previously. 
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This invention is illustrated in the Experimental Detail 
sections which follow. These sections are set forth to aid 
in an understanding of the invention but are not intended 
to, and should not be construed to, limit in any way the 
5 invention as set forth in the claims which follow 
thereafter . 

EXPERIMENTAL DPTMfrS 
10 Example 1: 

A self-splicing group II intron was transformed into a 
three-part ribozyme that site-specifically cleaves small 
oligonucleotide substrates with multiple turnover. The 
multi-component structure of the ribozyme facilitated 

15 quantitation of individual rate constants, interactions 
between reaction components (5'-exon, Domain l and Domain 
5) , and reaction stereospecif icity. Domain l and Domain 5 
may or may not be located on the same molecule. The 
efficient reaction catalyzed by this ribozyme appears to be 

20 limited by the chemical rate of hydrolysis. An Rp 
phosphorothioate substrate may be the preferred diastereomer 
of the ribozyme. 

Self -splicing group II introns are essential for the 
25 processing of many organella genes in plants and yeast 
(Umesono, 1989; Kuck, 1990). Group II introns are 
characterized by a conserved secondary structure that can be 
organized into six domains. The splicing reaction is 
initiated by attack of a 2'-hydroxyl group projecting from 
3 0 the bulged adenosine of Domain 6(D6), which attacks the 5'- 
splice site and generates a "lariat- intermediate (Peebles, 
1986; van der Veen, 1986). 

In group II introns, both steps of splicing can proceed if 
35 the initial nucleophile is a water molecule rather than the 
2' -OH of a bulged adenosine in D6 (Jarrell, 1988(a); 
Jarrell, 1988(b); van der Veen, 1987). In addition, 
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hydrolysis promotes the apparent first step of group II 
intron trans- splicing by 5'-exon-G oligonucleotides 
(Jacquier, 1986; Altura, 1989). Domain 5 (D5) is the most 
phylogenetically conserved region of the group II intron and 
5 it has been observed to act in- trans, promoting specific 
hydrolytic cleavage tat the 5' -splice site (Jarrell. 1988). 
It is clear that D5 and Domain 1 (Dl) are absolutely 
required for the first step of splicing (Koch, 1992). 
Recent kinetic analysis have examined the hydrolysis of 5'- 

10 exon from domains 1-3 when D5 is provided in-trans (Franzen, 
1993; Pyle, 1994). These studies showed that D5 associates 
with' high affinity (-300 nM) to the other active site 
components. The strength of D5 binding is especially 
surprising considering that this small RNA (-40 nts) is 

15 largely duplex in form and, having no apparent phylogenetic 
covariation with the rest of the intron. seems to bind 
without the aid of base-pairing interactions. This D5- 
catalyzed reaction obeys a straightforward Michaelis-Menten 
mechanism in which the chemical step of hydrolysis appears 

20 to be rate -limiting (Pyle, 1994) . 

M^-i mal Group TT introns that r e tain substrate binding and 
^t-Mytic activity 

Further deletion studies have demonstrated that group II 
25 introns will retain substrate binding capability and 
catalytic activity without certain subdomains. Please refer 
to the secondary structure of the yeast mitochondrial group 
II intron depicted in Figure 2A.2B. These experiments have 
shown that the group II intron will remain fully active for 
30 splicing with the removal of the D2b domain, the C2 domain, 
and the lower stem of the CI domain. Additionally, domain 
II and domain IV can be completely deleted and splicing 
activity of the molecule will remain intact. We have 
demonstrated that the combination of domain III and domain 
35 I results in a very efficient ribozyme. Domain III appears 
to speed up the kinetics of the chemical steps involved in 
splicing. 
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In order to independently evaluate the role of Dl, D5 and 
5'-exon sequences on active-site folding and catalytic 
activity, we constructed a modular group II intron ribozyme 
5 consisting of a 425 nt Domain 1 RNA (Dl) that binds small 
oligonucleotide substrates (s, 24-27 nts) in trans. 
Efficient cleavage at a single site is promoted by addition 
of a reaction cof actor (D5, 58 nts) that is also provided in 
trans (Fig. 4a) . This ribozyme is analogous to the 
10 Tetrahymena ribozyme except that the cof actor is a short RNA 
which is unchanged in the reaction instead of a nucleophilic 
guanosine that is consumed during reaction (Zaug, 1986 (a) ; 
Zaug, 1986(b); Pyle, 1994). With the reaction substrate and 
cof actors cleanly separated, it is now possible to hold any 
15 one component rate limiting in order to determine its 
precise role in binding or catalysis. The modular structure 
permits facile incorporation of mutations and single-atom 
changes for probing mechanism {Cech, 1992) . We have used 
this true ribozyme to begin construction of a more complete 
20 kinetic framework for group II intron catalysis and to 
determine the stereospecif icity of a reaction analogous to 
the first step of group II intron self -splicing. 

Kinetics of Substrate Cleavage by the Ribozyme: The 
25 reaction substrate (S) is a small oligonucleotide with a 
sequence analogous to the 5' -exon/intron boundary in the 
ai5g group II intron (Fig. 4a). In the presence of Dl and 
D5, this oligonucleotide is cleaved specifically at the CpG 
sequence that corresponds to the normal 5' -splice site (Fig. 
30 4b) . If either Dl or D5 is absent from the reaction, no 
cleavage of the substrate occurs. Previous studies showed 
that a 52-nt 5'-exon terminating in the CpG cleavage site 
was the minimal substrate necessary for trans -splicing or 
trans -cleavage by the ai5g group II intron (Jacquier, 1986; 
35 Jacquier, 1987) . m our study, only 17 nts of exon, 
encompassing both IBS1 and IBS2 (Fig. 4) were required for 
efficient cleavage. This substrate also contains 7-8 
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nucleotides of intronic sequence, including nucleotides (+3 
and +4) involved in the €-e' interaction believed to 
stabilize 5' -splice-site association (Jacquier, 1990). 

5 K„ 9 and K**/*« values for S were obtained using single- 
turnover kinetics where [S] « K„ s and [D5] was saturating 
relative to [Dl] (Table 1, line 1) according to standard 
kinetic procedures (Fig. 5. Table 1, line 1) (Pyle, 1994; 
Herschlag, 1990; Fedor, 1992). K* 5 determined from this plot 

10 is 6.3 Nm, implicating strong association between S and Dl 
(see IBS/EBS pairings, Fig. 4a). This is the first 
quantitative analysis of the strength of binding between the 
5'-exon and Dl. and the Kj approximates a dissociation 
constant for the IBS/EBS pairings plus any other 

15 interactions mediating intron/exon association. At very 
high concentrations of Dl (the plateau of the plot in Fig. 
5A) k obfl reflects the conversion of E-S -> E-P, where E 
represents Dl that is fully saturated with D5 cofactor and 
P represents reaction products. Under these saturating, 

20 single -turnover conditions, the apparent rate (defined here 
as k^J represents the rate of the chemical step and any 
conformational changes that take place within the E-S 
complex. The k^, value at saturating [Dl] was 0.021 min' 1 . 



25 
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Table 1 - Kinetic Parameters for Substrate Cleavage by the 
Group II Introns 



Substrate 


rate- 
lim 

compo 
nent 


(min- 1 ) 1 


(nM) 
s 


^max/^m 

(M- 1 min-») 1 ' 


^nax 

(minV 


thio 


S i synthetic) 


Dl 


0.21± 
0.001 


6.3± 
0.88 


(3.3±0.49) 
X 

10 6 


0.012 


1.8 


^(aynchecicj 


Dl 


0.032± 
0.0048 






0.009 
4 


3.4 


S| synthetic) 


D5 


0.019± 
0.0010 


870± 
110 


(2.2±0.30) 
X 

10* 







Table 1 legend: 

1 Described in legend to Fig. 5. 5 represents the maximum 
10 rate at saturation and represents the horizontal asymptote 
derived from the binding curves shown in Fig. 5. 
Theoretical k^ values at «D1 and »D5 were calculated from 
the data shown in Fig. 5 and found to be in agreement with 
a value of 0.026 min* 1 determined experimentally at 100 nM Dl 
15 and 9 /xM D5. K„. represents the apparent binding constant of 
the limiting component. These values were determined from 
the fit of the curves shown in Fig. 5 to an equation 
describing 1:1 bimolecular association". The standard error 
for the fit is shown. The appropriateness of the standard 
20 error calculation was confirmed using a jackknife 
approximation of standard error 12 . * k^/X {analogous to 
Kat/K values determined from the slopes of )^ bs vs. [Dl] or 
[D5] at very low concentrations. Reported error was 
propagated from the k^ and JC standard errors. ' The 
25 reported rates represent k ob , values for phosphorothioate 
substrates at saturation, where [Dl] * 100 nM and [D5] * 3 
/xM (analogous to the plateau of the curve in Fig. 5A) . Each 
value is the average of two trials. Like all reported k^,, 
variance is 18% (see legend to Fig. 5A) . The empirically 
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determined endpoint (after 7 reaction half-times) used in 
calculation of k^ was 82% for the transcribed substrate and 
48% for the synthetic substrate. .. The ratio of k M)t to Jc^ 
for phosphorothioate substrates at saturation. " This value 
5 is the k eb . at Dl and D5 saturation for cleavage of a 
transcribed all -phosphate substrate. 

For conditions under which [D5] was rate-limiting, values 
for K.**, KJ >i and k«,/K B , I>5 were determined by monitoring the 

10 cleavage rate of S in the presence of saturating Dl through 
a range of D5 concentrations (Fig. 5B, Table 1. line 3) . In 
this case, each molecule of S was bound to Dl and overall 
rate was limited by the concentration of D5. The k B ns is 870 
nM, k^/lC" is 2.2 x 10 min and Js is 0.019 "min . 

15 Values""of k™ from both Dl and D5-limited experiments are in 
close agreement (Table 1). The consistency of the data 
obtained under limiting Dl and D5 conditions suggests that 
under either condition, reaction proceeds through the same 
ternary complex for which the individual components have 

20 different affinities. 

Specific, efficient cleavage of small substrates by the 
group II ribozyme has now been observed in a kinetically 
characterized system. The results are consistent with 

25 kinetic behavior of larger, more complex constructs of the 
group II ribozyme has now been observed in a kinetically 
characterized system. The results are consistent with 
kinetic behavior of larger, more complex constructs of the 
group II intron (Jarrell, 1988; Pyle, 1994). The ribozyme 

30 construct can now be used to carry out detailed 
enzymological studies that examine reaction mechanism at an 
atomic level. Of particular interest is the 

stereospecificity of group II catalysis of the first and 
second steps of splicing (Lamond, 1993; Moore, 1993(a); 

35 Moore, 1993(b)). Group II introns are invoked as an 
autocatalytic model for the spliceosome and there are 
indications that group II introns and pre-mRNA splicing 
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share a common evolutionary ancestry (Cavalier -Smith, 1991) . 

Stereospecificity of Cleavage by a group II Ribozyme: 
Synthetic incorporation of a chiral phosphorothioate linkage 
5 results in an equal mix of Rp and Sp diastereomers. A 
synthetic substrate for the ribozyme was machine - 
synthesized, containing a single phosphorothioate linkage at 
the cleavage site. The ribozyme cleaves this substrate to 
an endpoint of -50%, suggesting that it has a preference for 
10 one of the two diastereomers (Fig. 6) . However, the 
fraction uncleaved by the ribozyme could represent RNA that 
is unreactive for other reasons, such as the presence of 
protecting groups left over from chemical synthesis. To 
determine the stereochemical identity of the fraction 
15 uncleavable by the ribozyme, the uncleaved fraction was gel- 
purified and examined by digestion with snake venom 
phosphodiesterase (SVPD) . SVPD is a 3S 5' exonuclease that 
readily cleaves phosphate and Rp phosphorothioate linkages. 
This enzyme is markedly inhibited by Sp phosphorothioates, 
20 cleaving them 2000X more slowly than their corresponding Rp 
isomers (Burgers, 1979). When a 5' -end labelled RNA 
contains an internal Sp phosphorothioate linkage, SVPD will 
readily digest the 3 ' -end of the strand and suddenly pause 
at the site immediately adjacent to the Sp linkage ( + 1G) . 
25 This apparent stall will accumulate and then diminish at 
long times as the SVPD slowly cleaves through it and then 
rapidly digests the remaining strand to mononucleotides. 
Three 5' -end labelled RNAs were subjected to cleavage by 
SVPD (Fig. 7) . An all -phosphate substrate was readily 
30 digested, with no apparent stalling at position 18, adjacent 
to the ribozyme cleavage site (Fig. 7, phosphate panel) . 
Digestion of the racemic substrate (never treated with 
ribozyme) resulted in prominent exonucleolytic stalling 
adjacent to the single phosphorothioate incorporated at the 
35 ribozyme cleavage site. After incubating the racemic 
substrate with ribozyme (for > 7 reaction half-times), the 
ribozyme -uncleavable fraction was digested with SVPD to 
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yield a predominant stall product that was much greater in 
intensity than observed for the racemic mix (Fig. 7, Sp 
panel) . For each lane, the number of counts in the stall 
band were normalized for the total number of counts in the 
5 entire lane (from 1-24 nts) . Using this value, we 
calculated the integrated population of the stall product 
over the duration of the timecourses shown. The integrated 
intensity of the stall product in the "Sp" panel (Fig. 7) 
was 1.0 relative to the integrated intensity of 0.59 in the 

10 racemic "Rp + Sp» panel. An integrated intensity 0.05 for 
the "phosphate" panel may be taken as a measure of 
background. Together with the fact that the racemic 
substrate is cleaved by the ribozyme to -50%, these results 
suggest that the group II ribozyme may preferentially cleave 

15 the Rp diastereomer. 

Taking into account the 50% endpoint for cleavage of the 
mixed substrate, rate constants for cleavage of the 
phosphorothioate substrate were readily obtained (Table 1) . 

20 A two-fold diminution in rate (the "thio-ef f ect" ) was 
observed for cleavage of a phosphorothioate linkage. The 
magnitude of this rate-decrease, upon substitution of 
phosphate with phosphorothioate, has been cited as evidence 
for a reaction rate limited by chemistry rather than 

25 conformational changes (Herschlag, 1991) . This would be 
consistent with results from previous studies of 5' -splice- 
site hydrolysis by a two-piece ribozyme. Correspondence 
between multiple -turnover and single -turnover kinetic 
parameters, together with the observations that K,, 1 * - K,, 05 and 

30 that rate was log-linear with pH led to the hypothesis that 
the reaction was limited by the chemical rate of hydrolysis 
(Pyle, 1994). 

in order to directly test the hypothesis that the ribozyme 
35 cleaves an Rp diastereomer, a new substrate was prepared by 
in-vitro transcription. T7 RNA polymerase incorporates only 
Rp phosphorothioates into RNA polymers (Griffiths, 1987) . 
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Using or-thio-GTP, we transcribed an RNA substrate containing 
an Rp phosphorothioate at every G residue, including the 5'- 
splice site (5'SS) . Just as the ribozyme specifically 
cleaves a phosphodiester linkage at the 5'SS sequence (Pig. 
5 8, lane A8), the ribozyme specifically cleaves a substrate 
containing an Rp phosphorothioate linkage at the cleavage 
site (Fig. 8, lane B8) . To confirm that the ribozyme 
cleaves both phosphodiester and phosphorothioate RNA 
substrates uniquely at the 5'SS sequence, the transcribed 
10 substrates and products were mapped with a battery of 
sequence- and stereo- specific nucleases (Fig. 8, lanes 4-7, 
9) . The substrate containing all -phosphate linkages was 
mapped next to the substrate containing phosphorothioate s at 
G resides (Fig. 8, Sections A and B) . The proper sequence 
15 at the cleavage sites was confirmed by digestion of 
substrates with C13, B-cereus and Tl endonucleases (Boguski, 
1979? Donis-Keller, 1977). Product oligonucleotides migrate 
at the expected position relative to specific nuclease cuts 
at -1C and +1G on the substrates. The endonucleo lytic 
20 digests described above, together with U2 digestion, 
confirmed the absolute sequence throughout both substrates 
(Fig. 8) . The migration of the ribose product with a 19-mer 
PI endonuclease digestion product suggests that the group II 
ribozyme leaves 3' -OH and 5 '-phosphate termini. The 
25 presence of a phosphorothioate at the cleavage site was 
confirmed by a gap in the PI nuclease ladder of the thio- 
substrate (lane 9), together with I 2 cleavage of the thio- 
substrate (i 2 data not shown) (Potter, 1983; Schatz, 1991) . 
Taken together, these data show that the ribozyme is 
3 0 selecting the proper cleavage site on both ribose and 
phosphorothioate -containing substrates. Thus, the ribozyme 
readily cleaves the Rp diastereomer of a transcribed 
substrate and does not select adjacent phosphate sites for 
attack. 

35 

Even though the ribozyme specifically cleaves transcribed 
substrate at an Rp phosphorothioate, the partial digests do 
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not provide quantitative evidence that the substrate 
sequence and linkage sulfur content are homogeneous at that 
site. Significant de-thioation at the splice-site linkage 
would result in conversion to phosphate and apparently 
5 efficient cleavage by the ribozyme. We sought a 
quantitative demonstration that the group II ribozytne 
cleaves a linkage that is uncleavable by PI endonuclease . 
To address these concerns, the thio-substrate was 
transcribed in the presence of or- 32 P-CTP. Because PI 
10 endonuclease cannot cleave RP phosphorothioates , a total PI 
digest of this substrate leads to three products in 
approximately equimolar ratios: *pC # *pCp.G and *pCpsGp 8 G 
where p a indicates a phosphorothioate linkage (Fig. 9A, lane 
3) . Because the intensity of the three bands are 
15 approximately equimolar, this indicates that the substrate 
has been faithfully transcribed by T7 RNA polymerase and 
that the cleavage site dinucleotide must be an Rp linkage 
with the same level of phosphorothioate incorporation as 
positions +5G and +6G (which end up in the trinucleotide) . 
20 When a sample of internally 32 P-labelled substrate is cleaved 
to completion by the ribozyme and then immediately digested 
with PI under similar conditions, two predominant bands are 
observed: *pC and *pCp,Gp.G (Fig. 9A, lane 4). The middle 
p*Cp a G band, corresponding to the ribozyme cleavage site, is 
25 considerably diminished in intensity while the *pC band is 
enhanced. Incubation of substrate with dl or D5 alone, 
followed by PI digestion does not result in a diminution of 
the middle *pCp 8 G band (Fig. 9A, lanes 4+5). These results 
quantitatively demonstrate that this group II ribozyme 
3 0 cleaves an Rp phosphorothioate in a transcribed substrate 
that is homogeneous in sequence and linkage at the cleavage 
site . 

For the transcribed phosphorothioate substrate, the cleavage 
35 kinetics were quantitated at Dl and D5 saturation, resulting 
in an apparent k,^ value approximately 3- fold slower than 
that of a phosphate at the ribozyme cleavage site (Fig. 9B) . 
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Like the "thio-ef f ect" observed using synthetic substrates, 
this value is consistent with an apparent rate- limiting 
chemical step {Table 1, column 7). The rate of cleavage in 
this homogeneously Rp case is the same (within error) as 
5 that of the 50% synthetic substrate that is cleavable by the 
ribozyme. This provides further evidence that the cleaved 
fraction in the synthetic experiments is predominantly Rp in 
configuration. 

10 The new group II ribozyme described here has provided a 
quantitative description of inter-domain interactions, 
individual rate constants and reaction stereospecif icity . 
This is the first example of a group II intron that has been 
transformed into a ribozyme construct that catalyzes an 

15 efficient, single reaction on small oligonucleotide 
substrates. The ribozyme acts as a sequence-specific RNA 
endonuclease and consists of two RNA subunits which are 
essential for catalysis. Kinetics for the reaction are 
unusually clean and suggest a stable, homogeneous 

20 combination of reaction components. Like the Tetrahymena 
ribozyme, the design of this group II construct facilitates 
biochemical and kinetic studies required for identification 
of active- site residues and tertiary interactions involved 
in catalysis. 

25 

The new three -part system shows that interactions between 
Domain 1 and the S'-exon are very strong and that the EBS- 
IBS pairings plus the first several intronic nucleotides are 
all that is required for 5' -splice site recognition in this 

30 reaction. Although compensatory base substitutions have 
been performed, this is the first study to quant itate the 
magnitude of 5-exon binding tot he intron (~ 6 nM) 
(Jacquier, 1987). The k^/Km 5 value for this ribozyme, (3.3 
x 10* M" 1 min 1 ) iB -20X slower than that of the Tetrahymena 

35 ribozyme, which is rate-limited by duplex association rather 
than chemical or conformational steps (Herschlag, 1990) . A 
slow chemical step (0.02-0.1 min" 1 ) appears to limit the rate 
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of RNA hydrolysis by ai5g group II ribozymes and as a 
result, Kmx/Xx* is bIowt than the rate of duplex formation. 
The results using this construct confirm that the 
interaction of D5 with the other components is very strong, 
5 especially given the ill-defined nature of its tertiary 
interactions. The kinetic parameters describing the kD5- 
dependent reaction are similar in this three -part system to 
those in a related two-part system that included D2 and D3 
and covalent attachment of 5'-exon to Dl (Pyle, 1994). 
10 However, the 2X higher [KCl] used here stimulates the rate 
in both systems by 2 -3 -fold. 

The ribozyme construct described here is relevant to 
splicing because a hydrolytic first step does not prevent 

15 the second step from occurring (Jarrell, 1988; and van der 
Veen, 1987) . Furthermore, the first step is believed to be 
rate- limiting in two-step splicing by the ai5g group II 
intron and the rate of reaction by our ribozyme is 
comparable to that of overall splicing under similar 

20 reaction conditions (Jarrell, 1988; Pyle, 1994; Jacquier, 
1987; Peebles, 1987; Chanfreau, 1994; and Peebles, 1993). 
The results observed here are not unique to the KCl 
conditions known to promote 5' -splice site hydrolysis. 
Efficient trans -cleavage of the small substrate occurs if 

25 KCl is replaced with any of other monovalent salts 
previously reported to promote the full splicing reaction. 

Our results demonstrate that a group II ribozyme 
preferentially cleaves an rp diastereomer at the 5' -splice 

30 site in a reaction that appears to be analogous to the first 
step of splicing. There are several possible 

interpretations of this result, given that there are two 
distinct D5 -catalyzed reactions that involve formation of 
the EBS-IBS pairings (Suchy, 1991). One possibility is that 

35 our system mimics the forward first step of splicing and 
that group II introns generally proceed with Rp 
stereospecificity during the first step. If this is so, it 
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is remarkable given that spliceosomal processing {for which 
group II introns are commonly invoked as models) proceeds 
with Sp stereospecif icity during the first step (Moore, 
1993; Maschhoff, 1993). It has been proposed that Rp 
5 stereospecificity, in either step of group II splicing, 
would indicate that group II introns are unrelated to the 
spliceosomal apparatus (Lamond, 1993) . However, it is 
possible that hydrolysis reactions or reactions where the 
substrate is presented in-trans may result in different 
10 stereospecificity than that of cis-splicing. This might 
occur if the e-f interaction (between the +3 ,+4 and 
+116, +117 nts of ai5g) cannot form properly in-trans or if 
branch formation requires participation of additional 
subunits not present in our construct. In any case, we have 
15 characterized a minimal first step reaction which proceeds 
with comparable rates and subunit binding constants to those 
of more complicated systems and is specific for cleavage of 
sequences immediately surrounding the 5' -splice site. Thus, 
if stereospecificity of our system turns out to be different 
20 than that in cis-splicing, it means that we can build on our 
system to identify additional structures solely responsible 
for Sp branch formation. It would remain relevant to pre- 
mRNA splicing because snRNAs are believed to have evolved 
from group II domains that became separated and started to 
25 function in-trans, much like our system (Sharp, 1985; Cech, 
1986) . It is therefore important to understand how the 
group II intron can be faithfully divided into trans- 
ribozyme constructs. 

30 A second possibility is that our system mimics a reversal of 
the second step of splicing and that the reaction may not 
proceed through the same active site as that utilized for 
branch formation. There are many indications that a 
reaction equivalent to the reverse of the second step of 

35 group II splicing can result in hydrolysis of linkages that 
follow 5'-exon sequences, particularly when certain 
structural features are lacking or under alternative 
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reaction conditions (Jarrell, 1988; Suchy, 1991). However, 
there are several indications that our reaction is not a 
reversal of the group II intron second step. D3 has 
previously been shown to be important for spliced exon- 
5 reopening, independent of its role in efficiency of the 
first step (Koch, 1992; Bachl, 1990). Our construct does 
not contain D3 . Perhaps most importantly, our ribozyme does 
not readily cleave an oligonucleotide composed of ligated 
exons (nts -19 to -1 of 5'exon followed by the first 9 nts 

10 of the 3'-exon, data not shown). When provided in excess, 
this oligonucleotide is significantly less inhibitory to our 
reaction than an equal amount of 5' -splice-site substrate, 
thus implicating effects on both binding and chemistry. 
Therefore, the intronic sequences of our normal substrate 

15 are important for catalysis (perhaps verifying the 
importance of the €-€' interaction) and our reaction mimics 
the forward first step in both sequence specificity and 
kinetic profile. 

20 In addition, like the first step of splicing this reaction 
is inhibited by the presence of a C residue immediately 3' 
to cleavage site. 

Relationship to group II intron self -splicing: 
25 The ribozyme described here behaves in a manner consistent 
with the first step of group II self -splicing and therefore 
serves as a representative model for first step processes 
(Table 2) . The fact that the reaction proceeds 
hydrolytically, rather than through branch-point attack, 
30 does not exclude the reaction from representing the first 
step for several reasons. The first step of splicing can 
proceed by hydrolysis without preventing the second step 
from occurring (van der Veen et al. f 1987; Jarrell et al . , 
1988(b)). Therefore, the hydrolytic first step is a 
3 5 significant catalytic strategy that group II introns can 
exploit in order to effect splicing. Spontaneous attack of 
water at specific internucleotide linkages is not a 
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favorable process in the absence of a catalyst. The 
uncatalyzed rate for hydrolytic cleavage of the 
phosphodiester linkage has been estimated to be 3 x 10 _9 min* 1 
(Herschlag and Cech. 1990) . Given the kinetic parameters 
5 reported here, the group II intron provides a 1 x 10 7 -fold 
rate -enhancement for hydrolytic cleavage at the 5 '-splice 
site. Like RNase P, the group II intron contains active 
site functionalities that specifically promote this 
reaction. Recent studies of in-vitro self -splicing show 

10 that, under all reaction conditions, the first step 
continuously proceeds through both competing hydrolytic and 
transesterif ication pathways. It is well documented that 
the first step of self -splicing is rate-limiting in ai5g 
group II intron (Jacquier and Michel, 1987; Peebles et al., 

15 1987; Chanfreau and Jacquier, 1993). Therefore, concrete 
kinetic parallels link the ribozyme reported here and the 
overall rate of group II intron self -splicing . 

The reaction reported in this study might be considered 

20 analogous to the "spliced-exon-reopening" or SER reaction 
promoted by group II introns in high concentrations of KC1 . 
This is unlikely given that ligated exons are a poor 
substrate for the ribozyme and that our results are not 
unique to the KC1 conditions exclusively thought to promote 

25 SER (Jarrell et al., 1988(b). Efficient trans -cleavage of 
the small substrate occurs if KC1 is replaced with any of 
other monovalent salts previously reported to promote the 
full splicing reaction. Another distinguishing feature is 
that the SER reaction is unaffected by C substitution at the 

30 +1 nucleotide (Peebles et al., 1993). The ribozyme reaction 
described here, like the first step of the full splicing 
reaction, is inhibited by the presence of a +1C, while any 
other nucleotide can function at that position. Although 
the reaction reported here is clearly distinct from SER, we 

35 consider SER to be a reaction related to the first step of 
group II intron splicing. If a group II intron does not 
rapidly release ligated exon products, they may be 
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susceptible to another nucleophilic attack at the splice 
junction. A common misconception is that the SER reaction 
is a reversal of the second step of splicing. Strictly 
speaking, this is impossible given that the second step of 
5 group II intron splicing occurs exclusively (and quite 
rapidly) through a transesterif ication in which the leaving 
group is the 3' -hydroxy 1 of the lariat or linear intron 
(Jarrell et al. f 1988(b). Water is not the leaving group in 
the forward reaction, so water cannot be the attacking 
10 nucleophile in a reaction that is truly the reverse. 



Table 2: Kinetic effects of base changes at the 
ribozyme cleavage site 



30 



SubstrateY 


k ob- (min-l)§ 


relative rate 


+1G (WT) 


0.036t 


1.0 


+ 1A 


0.053t 


1.5 


+ 1U 


0.028t 


0.79 


+ 1C 


0.0019* 


0.05 


ligated exons 


0.0068* 


0.19 



Table 2 Legend: 

Y The sequence of transcribed substrates used to monitor 
35 mutational effect at the 5 '-splice site is 5'- 
GGUGUGGUGGGACAUUUUC^NAGCGGUU , (Seq. ID No. 6) where N marks 
the +1 nucleotide at the cleavage site. The sequence of the 
ligated exons substrate is 5 ' -GGAGUGGUGGGACAUUUUC'ACUAqGUft, 
(Seq. ID No. 7> where the underlined nucleotides are 3'-exon 
40 sequences. Rates are not affected by identity of the first 
three 5' -nucleotides . 

§ All rates were obtained simultaneously so it is exptectd 
that their relative error is very low (< 5%) . The absolute 
error for k oba values reported in this study is + 19% with 
45 95% confidence (see legend, Table 1) . 

t Pseudo- first order plots were linear for at least four 
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analogous to that in which 32 P-exD123 was treated as 
substrate, while D5 was added in excess as the enzyme for 
hydrolytic cleavage of the 5' -splice site (Pyle and Green, 
1994) . 

5 

Time-courses of the reaction reveal that D56 forms a high- 
molecular weight branch, becoming attached to D123 (note 
bold band D56/D123, Figure 11A) . At the same time, the 
previously observed 5' -splice site hydrolysis reaction 

1C competes with branching (bands D123 and D12 , Figure 11A) 
(Jarrell et al . , 1988a). During both the branching and 
hydrolysis reactions, 5' -exon is liberated (Figure 11A) . 
The quantity of hydrolytic cleavage at the 5' -splice site is 
the same for D56 as it was for D56A- , but the latter 

15 molecule is unable to form a branch with exD123 (compare D56 
and D56A- panels, Figure 11A) . This establishes that a 
bulged adenosine in D6 is essential for 2'-hydroxyl group 
attack, as previously observed for full-splicing reactions 
(van der Veen et al., 1987). The lack of branched product 

20 for D56A- suggests that 2'-hydroxyl groups elsewhere on D6 
do not attack the 5 '-splice site and that the fidelity of 
branch-site selection is not compromised in our system. 

Once it was apparent that D56 and exD123 underwent 
25 transesterification, it was important to establish reaction 
conditions that were appropriate for careful analysis of 
this reaction. A survey of monovalent and divalent ion 
effects gave results consistent with other studies 
indicating that high concentrations of ammonium salts (0.5 
30 to 1.5 M) stimulate branching relative to hydrolytic 
pathways for 5' -splice site cleavage (data not shown) 
(Peebles et al., 1987; Jarrell et al., 1988b). Ammonium 
sulfate readily stimulated branching, but only after long 
and irreproducible lags in reactivity. Ammonium chloride 
35 also promoted branching, but without the concomitant lags 
and poor kinetic behavior associated with ammonium sulfate. 
Although different divalent ions, including Ca 2+ , Mg 2 * and Mn 2 * 
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(100 mM) were examined, only Mg 1 * supported activity, 
provided that at least 50% of the divalent ion was Mg J \ Mn" 
could substitute for some of the divalent ion requirement. 
Based on these findings, the buffer salts chosen for this 
5 study were 500 mM NH t Cl and 100 mM MgQl , which are the 
lowest salt concentrations capable of promoting efficient 
branching. Temperature was also varied and 45oC was found 
to be the optimal condition for this reaction. 

10 When a single -turnover ribozyme reaction is highly pH- 
dependent, it is likely to be limited by the rate of 
chemistry (Dahm et al. ( 1993; Pyle and Green, 1994; Zaug et 
al., 1994). Conversely, a lack of pH dependence under 
these conditions suggests that reaction rate is limited by 

15 conformational change. The D5-dependent hydrolysis of the 
5- -splice site was previously shown to be highly pH 
dependent and limited in rate by the chemical step of 
reaction (Pyle and Green, 1994) . Although it contains a 
relatively intact D6 duplex, DS6A- reacted in an almost 

20 identical manner to that previously reported for D5 (Figure 
11A) (Pyle and Green, 1994). D56 and D56A- hydrolysis 
reactions of the 5' -splice site are both highly pH-dependent 
and, by analogy to D5 reactions, they are probably limited 
by the rate of chemistry. The maximum saturating rate of 

25 hydrolysis by both D56 and D56A- at pH 7.0 under standard 
reaction conditions (k c « for hydrolysis) is 0.018 ± 0.0066 
min" 1 (determined as described for k^ values in the legend 
to Figure 13, inset). However, the branching reaction 
promoted by D56 is not pH dependent from pH 5.5 to 8.0 

30 (Figure 11A and data not shown) . This indicates that the 
rate- limiting step for branching is something other than 
chemistry . 

In addition to having different rate-limiting steps, the two 
competing reactions promoted by D56 have markedly different 
35 profiles of product evolution. The hydrolysis reaction 
proceeds slowly and exponentially, as expected for a pseudo- 
first order reaction similar to that catalyzed by D5 alone 



WO 96/22689 



PCT/US96/01337 



-63- 

( Figure 11B) . Conversely, the branching reaction is very- 
rapid in the first few minutes, resulting in a burst of 
activity that halts and then decreases before all the exD123 
precursor is consumed (Figure 11C) . This effect is 
5 particularly pronounced at pH values greater than or equal 
to 7.5, at which point hydrolysis of the 5' -splice site 
dominates over branching (Figure 11A) . Late in reaction, 
the observed decrease in branched molecules could be due to 
instability of the 2' -5' linkage, causing it to break and 
10 produce linear D123 . However, this is ruled out by 
identical kinetics of hydrolysis for the D56A- and D56 
reactions (Figure 11A) . If breakage accounted for the 
decline in branched species, D5 6 -mediated hydrolysis would 
appear to go faster, because it would possess an additional 
15 pathway for generating the D123 product. The plateau 
observed for D56 branching could be due to a limited 
population of conformations that are competent for 
transesterif i cat ion. Alternatively, the observed behavior 
could be due to a branching reaction that is reversible 
20 while hydrolysis is not. If transesterif ication were 
reversible, the population of branch would plateau early 
and, in the presence of an irreversible competing reaction, 
it would appear to eventually decline and lose ground to an 
irreversible process that proceeds to completion. By this 
25 model, the hydrolysis reaction behaves as a trap or sink, 
through which all of the material ultimately reacts. If the 
latter model is correct, there are three reactions going on 
at the 5' -splice site: branching, reverse branching and 
hydrolysis. One might be able to set up a complicated 
3 0 mathematical model describing this phenomenon and fit the 
data to it. This would be difficult given that one cannot 
obtain an individual reaction endpoint for each of the many 
competing reactions observed here. Instead, we decided it 
was better to simplify the experiment so that more 
35 straightforward kinetic analysis could be applied. Our 
approach involved redesigning the experiment so that 
branching could be examined without interference from the 
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competing hydrolysis reaction. 

Branching is studied in the absence of hydrolysis: One way 
to study branching independent of any other reaction is to 
5 follow the conversion of 32 P-D56 into branched species <"P- 
D56/D123). Because the nucleophile during branching is the 
bulged adenosine 2' -hydroxy 1 on DS6, D56 is changed in the 
course of reaction, becoming incorporated into a higher- 
molecular weight molecule. Although the hydrolysis reaction 

10 is occuring, there are several reasons why it should not 
influence the kinetic parameters of branching determined in 
this experiment. Provided that exD123 is unlabeled and in 
excess, the hydrolysis reaction will be undetectable because 
D56 catalyzes that reaction in a multiple-turnover manner, 

15 without becoming incorporated into the products. The 
nucleophile in the hydrolysis reaction is water or metal - 
coordinated hydroxide, neither of which is covalently 
connected to D56. The noncovalent binding properties of D56 
are similar to those of D5 (see following sections) and, by 

20 analogy to D5, the on-and off-rates of D56 from exD123 are 
expected to be fast (Pyle and Green, 1994; Michels and Pyle, 
1995) . Therefore, trace quantities of "P-D56 (InM) and 
excess exD123 (3 jiM) should react in a manner purely 
descriptive of the branching reaction. As anticipated, 

25 under these reaction conditions, the only process observed 
is covalent attachment of "P-D56 to D123 (Figure 12A) . In 
a denaturing gel, the product of that reaction comigrates 
with the branched species observed when "P-exD123 is labeled 
and D56 is in excess. 

30 

The branched D56/D123 molecule was physically characterized 
in two ways: PI nuclease was used to confirm the presence 
of a pA(2' -pG)3' -pU branched trinucleotide and the position 
of the branch was confirmed by sequencing. To characterize 
35 the branched trinucleotide, D56 was transcribed with 32 P-UTP 
and then reacted with exD123. The resulting "P-D56/D123 
branched product was isolated and digested to completion 



WO 96/22689 



PCT/US96/01337 



-65- 

with PI nuclease, which cannot cleave a branched 
trinucleotide (Reilly et al. f 1989). Products of this 
reaction were then subjected to standard electrophoretic 
techniques which revealed the presence of an abundant 
5 branched trinucleotide (Figure 12C) (Moore and Sharp, 1993), 
formed through the proper combination of nucleotides. The 
sequencing experiments, which were performed using 
phosphorothioate interference {Gish and Eckstein, 1988; 
Schatz et al., 1991), resulted in ladders that terminated at 
10 the bulged A branch-point, rather than at the 3' -end of D6 . 

Branching cannot proceed to completion: Having simplified 
the transesterification (branching) reaction to a single 
conversion of unbranched to branched species, and having 
15 characterized the products of this transformation, it was 
then possible to derive kinetic parameters describing the 
reaction. The conversion of "P-D56 to" P-D56/D123 was 
monitored as a function of time (Figure 12A) . Quantitation 
of the data revealed that, although the initial reaction is 
20 rapid and follows a single exponential, it proceeds to only 
-50% completion (Figure 12B) . Again, there are at least two 
explanations for this observation: Either only 50% of the 
D56 molecules are catalytically active, or the reaction is 
reversible. A hallmark of reversible reactions is that they 
25 will proceed to an apparent endpoint of less than 100% 
completion because a steady-state equilibrium is established 
between starting materials and products (Jencks, 1987) . 
Experiments were therefore designed to test for this 
possibility and to measure the rate of any reverse reaction. 
3 0 Approaches for examining the reverse reaction are described. 
However, before proceeding with a discussion of the reverse 
reaction, it is important to point out mechanistic 
information that was obtained from timecourses of the 
forward reaction (Figure 12A, B) . 

35 

Branching kinetics reflect D56 tertiary interactions: Rates 
of the forward branching reaction were obtained by 
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correcting the data for an observed 50* reaction endpoint 
(Figure 12B) and then plotting the evolution of product in 
a semilog format (Figure 13B) . The data fell on a straight 
line and observed rates of reaction (k„J at a given 

5 [exD123] were taken from the slope of plots such as the one 
shown. The fact that the data are linear through four half- 
lives of reaction indicates that the reaction is "pseudo- 
first order", meaning that the "P-D56 and exD123 react as if 
they were one molecule, in a complex. Linearity in this 

10 type of plot also indicates that the observed forward 
reaction is kinetically "well-behaved", representing a 
single process that is uniform from early to very late in 
the reaction. As shown in Figure 11. the rate of this 
reaction is not pH dependent, indicating that it is not 

15 limited in rate by the chemical step of reaction. In 
addition, the linear kinetic behavior takes place after a 
tiny burst of initial reaction, representing 4% of the total 
population of molecules (y-intercept. Figure 13B) . Taken 
together, these data indicate that the rate-limiting step in 

20 the forward reaction for most of the molecules is not 
chemistry and is likely to represent the rate of 
conformational change of the D56/exD123 complex into a 
catalytically active state. 

25 Although the forward rate of branching <k ob .) does not 
represent chemistry, the single-turnover design of this 
experiment and pseudo-first order kinetic behavior allowed 
us to use k ob , values for estimating the relative affinity 
between D56 and exD123. This was accomplished by 

30 determining the rate of branching (k^.) at a variety of 
[exD123]. The k^, values obtained from a collection of 
semilog plots (such as that shown in Figure 13B) were 
plotted as a function of [exD123] . As previously observed 
for D5-mediated reactions (Pyle and Green, 1994), these 

35 plots show saturable reaction behavior and fit a standard 
bimolecular binding curve (Figure 13A) . Kinetic parameters 
for the branching reaction were obtained from the fit of 
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this binding curve to the data. The k^, or maximum rate of 
the forward reaction at saturation is 0.041 ± 0.0040 min" 1 
and is 510 ± 160 nM. This K,,, for D56, as revealed through 
single-turnover branching reactions is remarkably similar to 
5 that of D5 (270 ± 25 nM) , determined through single turnover 
kinetics, multiple turnover kinetics and direct binding 
assays (Pyle and Green, 1994) . D56A- shows similar binding 
parameters, revealed through kinetic analysis of the 
hydrolysis reaction. These data suggest that D6 does not 
10 contribute additional tertiary interactions for 
stabilization of the group II active site and that most of 
the inter- domain binding interactions are contributed by D5 
for both branching and hydrolysis reactions. 

15 Reverse-branching requires a 5' -exon nucleophile: The 
reverse reaction is expected to occur when the 3 ' -hydroxy 1 
group of the 5 '-exon acts as a nucleophile, attacking the 
phosphate at the branch-site and releasing the 2'-hydroxyl 
of D56 as the leaving group. The first experiments 
20 conducted using internally labeled exD123 suggested that, 
because branched material decreased after an initial 
plateau, the branching reaction waB fully reversible {Figure 
11C) . This was supported by the 50% endpoint observed in 
reactions of * 2 P-D56 with exD123 (Figure 12B) . It has been 
25 proposed that the 5' -exon is not released immediately 
following attack at the 5' -splice site, but rather that it 
remains noncovalently attached to the intron (Jacquier and 
Rosbash, 1986) . In the bound configuration, the associated 
exon can either perform the second step or, as observed 
30 here, it can attack the branch in a reverse of the first 
step. However, if the exon were removed and purified, it 
might be expected to have difficulty re-entering the active 
site. 

35 In the previous experiments (Figure 11C and 12B) , the 5'- 
exon was the long, 293 -nucleotide exon normally attached to 
D123. To study reverse -branching, we generated this exon 
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from a large-scale reaction. It was gel-purified and then 
added to a reaction containing the branched "P-D56/D123 
material. However, after isolation and gel -purification, 
the 293- nucleotide 5'-exon was unable to react. Without 
5 the context of the adjacent intron, or perhaps as a 
consequence of purification, the long exon could not 
function as a nucleophile, possibly having adopted an 
unfavorable conformation that occluded its 3 '-end. 
However, a short 5 '-exon that contains all sequences 
10 necessary for interaction with the intron readily stimulated 
reverse branching after addition to 32 P-D56/D123 . This RNA, 
which was used in all subsequent kinetic analyses of reverse 
branching was 17 -nucleotides in length, consisting of intact 
IBS1 and IBS2 sequences and terminating with the 3' -hydroxy 1 
15 group of -1C (the last nucleotide of the 5' -exon). The 
sequence of this oligonucleotide is 5 ' -CGUGGUGGGACAUUUUC-3 ' 

(SEQ I.D. No. ). Reactions of ai5g and other introns 

have been observed in the presence of short 5'-exons 
containing only the IBS sequences (Figure 10A) (Suchy and 
20 Schmelzer, 1991; Michels and Pyle, 1995). This is made 
possible by the fact that EBS-IBS interactions between exon 
and intron are very strong, with K d determined to be ~ 10 nM 
in the case of the ai5g intron (Michels and Pyle, 1995) . 
For studies of reverse -branching, "P-D56/D123 was isolated 
25 from a denaturing polyacrylamide gel and purified according 
to standard techniques. Using the same reaction buffers and 
incubation temperature as the forward reaction, trace 
quantities of branched material (1 nM) were reacted with an 
excess of 5' -exon ( Figure 14A) . During the ensuing 
3 0 reaction, the branched material rapidly gave rise to free 
"P-D56 in a manner that reflected the strength of 5' -exon 
association. Timecourses of the reverse reaction were 
conducted in the presence of varying concentrations of 5'- 
exon. Consistent with strong interactions IBS-EBS 
35 reactions, the efficiency of reverse -branching (see below) 
decreased only at concentrations of 5' -exon lower than 10 
nM, which serves as a rough estimate of K<,. 
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In the absence of 5'-exon, the branched species slowly 
hydrolyzed at the 2' -5' linkage (Figure 14A) , but this 
process was extremely inefficient. Although addition of 5'- 
exon to 32 P-D56/D123 stimulates efficient release of "P-D56, 
5 this same apparent behavior would be observed if 5'-exon 
served only to enhance folding of D123, facilitating 
hydrolytic cleavage of the 2' -5' -linkage. Before reverse- 
branching could be conclusively established, it was 
necessary to observe covalent attachment of 5'-exon onto 

1C D123. This would prove that a transesterif ication, rather 
than a hydrolysis reaction, had occurred. To this end, 
5'-exon was labeled with 32 P and incubated with the branched 
species. As expected, upon reaction with "P-D56-D123 , the 
32 P-5'-exon shifts to the appropriate higher molecular weight 

15 upon reaction (727 nts) , with concomitant appearance of 
proportional quantities of "P-D56 (Figure 14B) . Together, 
these results establish that the D56/D123 branch covalently 
reacts with 5-exon, releasing the 2'-hydroxyl of the bulged 
adenosine as the leaving group. 

20 

The reverse -branching reaction is rapid: Timecourses of 
reverse -branching reveal a rapid burst behavior indicative 
of a large population of highly reactive molecules. The 
kinetic data fall in two distinct phases, and rates of 

25 reverse branching were calculated by fitting the raw data to 
a biexponential equation (Figure 14C) . This type of 
analysis allows one to obtain the rates of two different, 
simultaneous reactions. The burst phase of reaction, 
representing 19% of the molecules, proceeded with a rate of 

30 1 min* 1 . The remaining population of molecules reacted with 
a rate of 0.01 min" 1 . For both phases of reaction, 32 P-exon 
shift experiments analogous to those shown in Figure 14B 
confirmed that free D56 evolution was the result of 5'-exon 
attack, and not slow hydrolysis of the branched linkage. 

35 The fast phase of reaction (k^ = 1 min' 1 ) represents a lower 
limit for our estimate of the chemical rate for reverse 
branching. By contrast, the second slow phase of reaction. 
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like the forward reaction, is likely to represent the rate 
of conformational change from an inactive state to a state 
that is capable of reverse branching. It is highly 
significant that a higher fraction of molecules is 
5 conformational^ predisposed to burst and react quickly in 
the reverse reaction (19*) in comparison to the forward 
reaction (4%) . This suggests that the branched molecules 
are more preorganized and that the 2' -5' linkage, which 
brings the 5' -splice site close to catalytic functionalities 

10 on Domain 5, helps more of the population to assume a 
conformation capable of efficient reaction. These data 
confirm that, as indicated by behavior of the forward 
reaction, the reverse reaction is fast and can compete 
effectively during the first step of splicing {Figure 15) . 

15 were it not for a highly efficient second step of splicing, 
which engages every free 5 r -exon as soon as it forms, the 
forward reaction would be unable to proceed to completion 
through branching. 
DISCISSION 

20 One of the most distinctive features of self -splicing by 
group II introns is that the first step of splicing takes 
place through branch-point attack, resulting in production 
of a lariat intron (Peebles et al., 1986; Schmelzer and 
Schweyen, 1986; van der Veen et al., 1986). Here we show 

25 that this reaction occurs simultaneously with two other 
reactions, all of which have important implications for the 
mechanism of productive group II self -splicing and intron 
mobility. We observe that the forward branching reaction 
cannot proceed to completion on its own. The reason for 

30 this is that a rapid reverse reaction becomes dominant once 
a population of branched molecules is established. The 5' - 
exon that is released during the forward branching reaction 
is a highly reactive nucleophile in the reverse reaction, 
readily debranching the 2 '-5' linkage. Our data suggests 

35 that the reverse reaction is actually more facile than the 
forward reaction due to greater proximity of active-site 
components. Another reaction that competes during forward 
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branching is hydrolytic cleavage of the 5' -splice site. 
Although this route is slower than branching, it is 
irreversible and is observed to drive the first step to 
completion in the absence of other reactions (such as the 
5 second step of splicing) . 

An alternative pathway for the first step of splicing: The 
hydrolytic pathway makes a significant contribution to the 
forward reaction and to splicing in general. Although it is 
10 slow, it is limited by the chemical step and not by 
conformational heterogeneity of the intronic components 
(Pyle and Green, 1994). Catalyzed by D5, this reaction 
appears to serve as a default reaction that promotes proper 
5' -splice site cleavage regardless of intronic 
15 conformational problems. The reaction appears to take place 
because of an inherently large strain on the scissile bond 
of the 5' -splice site, rendering that linkage susceptible to 
attack by nucleophiles such as water or Mg^-bound hydroxide. 
The hydrolysis reaction is irreversible, and might therefore 
2 0 be an important means by which genes ensure that forward 
splicing takes place without risking intron mobility into 
other parts of their genome. In- vivo, there may be proteins 
that facilitate the hydrolytic splicing pathway and others 
that suppress it in favor of branching. Proteins that may 
25 belong to the latter class have already been implicated in 
intron mobility, a process that would benefit from splicing 
through branching (Lambowitz, 1989; Lambowitz and Belfort, 
1993) . 

30 Conformational limitations on 5 f -splice site 
transesterif ications: When monitored in the absence of 
hydrolysis, the kinetics of the forward branching reaction 
are remarkably linear for up to four half-times of reaction 
and they do not display lags when plotted in a semi -log 

35 format. This indicates that there is no set population that 
is particularly slow to fold and that each molecule has the 
capability to react through branching. The fact that the 
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observed rates (slopeB of semilog reaction plots) do not 
represent the chemical step and that the reaction shows 
saturable behavior at 510 nM exD123 indicates that D56 
readily binds exD123, but that in the bound configuration it 
5 samples a set of conformational states, only a subset of 
which are competent for branching. For the forward and 
reverse reactions, the population of D56/exDl23 complexes 
that react immediately (4 and 19%, respectively) represent 
the fraction of molecules that started in the proper 

10 conformational state for reaction. Once this population is 
consumed, the remaining D56/exD123 molecules react more 
slowly, as if limited by the rate of a single conformational 
change. It is important to point out that this is probably 
not a directed, single conformational change from one state 

15 to another. Macromolecules typically sample many 
conformations in a random fashion. The slow rates of 
transesterification following bursts of activity may 
therefore represent the average time of stochastic sampling 
before the complex (either D56,exD123 in the forward reaction 

20 or 5' -exon,D56/D123 in the reverse) reaches a conformation 
that is kinetically competent to react. 

It is interesting that the X. for D56 (510 nM) , as reflected 
through the kinetics of the branching reaction (Figure 13), 

25 is similar to that of D5 (270 nM) , as reflected through 
direct binding assays and catalysis of the hydrolysis 
reaction (Pyle and Green, 1994) . This indicates that D6 
does not contribute additional tertiary interactions beyond 
those already provided by D5. Our data support the notion 

3 0 that D5 functions to dock into the core and position D6 
correctly such that the bulged adenosine is in position to 
attack the 5' -splice site (Dib-Hajj et al., 1993). Further 
evidence for the idea that D5 behaves similarly in the 
presence and absence of D6 comes from rates of hydrolytic 

35 cleavage at the 5' -splice site. Just as D5 functions as a 
true enzyme for exD123 hydrolysis (Franzen et al., 1993; 
Pyle and Green, 1994), so does D56. The rates of hydrolysis 
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catalyzed by D56 and D56A- (0.018 min' 1 ) are only slightly 
slower than that of D5 alone (0.041 min* 1 ) (Figure 15) . The 
small difference in their hydrolysis rates may suggest that 
D6 shields the 5' -splice site from attack of water, or that 
5 it sterically occludes the approach of solvent -based 
nucleophiles, but such effects are probably small given the 
minor differences in hydrolytic rate between D5 and DBS 
analogs. 

10 Rapid transesterification by a group II intron: It is 
important to consider why the branching reaction proceeds so 
readily in reverse. After all, a 2'-hydroxyl group is a 
better nucleophile than a 3'-hydroxyl group (Sulston et al., 
1968) so if nucleophilicity were the only consideration, the 
15 reaction would probably proceed only in the forward 
direction. However, the forward reaction creates a higher- 
energy species than the reverse reaction. The Y-shaped, or 
lariat species created after the first step places 
catalytically essential Domains 5 and 6, together with the 
20 3 '-splice site, in immediate proximity to intronic portions 
of the 5 '-splice site. The high-energy nature of the 
branched species is also supported by our data, which 
indicate that the 2' -5' bond is weaker than 3' -5' linkages 
that might replace it (Figure 14A) . It has been proposed 
25 that the constrained lariat structure formed after the first 
step of splicing drives the second step (Jacquier and 
Jacquesson-Breuleux, 1991) . Here we show that this high- 
energy intermediate can also drive the first step backwards, 
so it is not necessarily a determinant in second-step 
30 efficiency. Conformational preorganization of the branched 
species is expected to be especially pronounced with intact 
group II introns because, in addition to covalent 
connection, a lariat intron circularly constrains active- 
site residues so that they are held in spatial proximity. 
35 The efficiency of reverse branching in this two-part system 
should not be attributed to the fact that catalytic domains 
are located on one molecule in the reverse reaction, while 
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they are divided on two molecules in the forward direction. 
This entropic effect is countered by the fact that all of 
the reactions reported here were studied at saturation. 
That is. even though D5 and D6 were on a separate molecule 

5 in the forward reaction, concentrations were manipulated 
such that they were fully bound to exD123 . This was 
confirmed by the pseudo- first order (acting as one molecule 
instead of two) kinetic behavior of the complex. It is 
potentially harder to ensure docking of D56 when these 

0 domains are physically connected to D123 in-cis. This is 
because intervening RNA conformations can prevent an RNA 
molecule from assuming the docked state and it would be 
difficult to kinetically detect this effect. 



5 The reverse reaction appears to be so rapid not because it 
is necessarily faster in the chemical step (although the 
lower limit for that step is a very rapid 1 min'M . but 
because more of the molecules are conformationally competent 
to chemically react. Note that the burst in the forward 

20 reaction represents approximately 4% of the total 
population, while the burst in the reverse reaction is 19% 
of the population. In both forward and reverse reactions, 
bursts are followed by slower reaction phases (k conf(lwl , = 
0.041 min 1 ; k cmtitev > = 0.010 min' 1 ) 

25 that represent the rate of conformational change into 
catalytically competent forms (Figure 15) . Thus, there is 
a direct catalytic advantage conferred by constraining the 
branched molecules. In-vivo there are probably proteins 
that facilitate folding into fully-active conformations for 

3 0 reverse or forward reactions, depending on whether the 
protein is designed to stimulate splicing or intron mobility 
(Lambowitz and Belfort, 1993). 

In addition to providing a mechanistic picture of competing 
35 reactions during the first step of splicing, this study 
reveals for the first time that group II introns have the 
innate ability to catalyze fast reactions. For the most 
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active population of molecules, the reverse of the first 
step proceeds with rate of 1 rain* 1 , and judging from the 
small burst observed in the forward direction, its rate is 
at least that fast. A rate constant of this magnitude is 
5 observed for the action of many protein enzymes during other 
biologically significant reactions, such as DNA cleavage by 
restriction endonuclease Eco Rl (Walsh, 1979) . it is also 
similar to the h M values observed for other ribozymes (Pan 
et al., 1993). Previous studies indicate that the rates of 

10 group II intron- catalyzed reactions are relatively slow, 
with best half-times of - 10 minutes under ionic conditions 
similar to those used here (Michel et al., 1989). It is now 
clear that these apparently slower rates are due to 
constructs that react through hydrolytic pathways or that 

15 are limited in overall rate by unfavorable conformations. 
Without examining early sets of timepoints during 
transesterification, it would have been difficult to discern 
if these reactions were fast. In the absence of kinetic 
data demonstrating that a chemical step of RNA catalysis by 
20 group II introns could be efficient, there remained the 
possibility that intron-associated proteins were required 
for composition of the active site in- vivo. This study 
demonstrates that, given the proper conformational state, 
RNA structures within group II introns can catalyze 
25 reactions on a biologically feasible timescale. 

Relevance of a facile back-reaction during the first step of 
splicing: While two-step splicing would obscure the 
apparent effects of branch reversion, the reaction is 

30 nonetheless important. In fact, a reversible branching 
reaction confers significant advantages. First, it 
provides a check for the fidelity of 5' -splice site 
cleavage. If an improper 5' -splice site were chosen during 
branch-point attack, the intron would lose interactions 

35 required for the second step of splicing, such as the 
interaction between first and penultimate nucleotides of the 
intron (Chanfreau and Jacquier, 1993? Peebles et al., 1993). 
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Facile branch reversion, particularly when mis-branching 
slowed the second step, would give the intron another 
opportunity to start over again and find the correct 5' 
splice site. This proof-reading mechanism would need to be 

5 verified by experiments showing that mis-cleaved 5'-exons 
can attack and reverse formation of a cryptic branched 
species. Reversibility might also prevent mis-splicing 
during transcription. Any branching reaction that takes 
place before important residues in the 3'-exon are free of 

10 template and polymerase would readily reverse. Thus, until 
the RNA is fully competent for an accurate second step, 
branching reversal will effectively prevent splicing. A 
similar mechanism may help control the fidelity of splice- 
site selection during spliceosotnal processing. Even in that 

15 case, the formation of the 2' -5' link after the first step 
of splicing might create a higher-energy intermediate in 
which a reverse reaction is favorable. However, in that 
case the control of this reaction is probably mediated by 
snRNPs and associated factors rather than by intronic 

20 substructures. The spliceosomal protein Prp 16 has recently 
been implicated in the controlling the fidelity of branch- 
site selection (Burgess and Guthrie, 1993). Similarly, a 
protein may control fidelity of 5' -splice-site selection by 
exploiting the potential reversibility of spliceosomal 

25 branching reactions. 

The second step of group II intron self -splicing is known to 
be a very rapid reaction that proceeds almost entirely to 
completion, while the first step of splicing is known to be 

30 rate-limiting (Jacquier and Jacquesson-Breuleux, 1991) . Our 
results provide a clear reason for the rapidity of the 
second step: it serves as a kinetic trap for those 
molecules that succeed in branching. Before the first step 
of splicing can reverse, the second step rapidly ensues and 

35 the full splicing reaction is driven to completion. Thus, 
there is necessarily kinetic coordination of the first and 
second steps of splicing in group II introns. While 
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precision timing of snRNP binding and release helps to 
coordinate the two steps of spliceosomal processing, RNA 
substructures in the group II intron have adopted the 
thermodynamic balance and kinetic control necessary for 
5 driving a series of competing reactions to a common 
conclusion, resulting in fully spliced RNA. 

One of the most intriguing consequences of efficient 
branching reversal is its potential role in intron mobility. 
10 Group II introns are infectious elements that have scattered 
widely through the genomes of many different organisms 
(Lambowitz, 1989). There are several models for the 
mobility of self -splicing introns, some of which do not even 
involve the splicing of intronic RNA (Lambowitz and Belfort, 
15 1993) . However, one model for intron mobility proposes that 
liberated lariat introns can bind other RNA molecules, 
recognizing sequences with a high degree of complementary to 
their own EBS regions. Through a series of two successive 
reverse transesterif ications, the intron is believed to 
20 insert itself into the host RNA. Group II introns often 
contain open reading frames for proteins with reverse - 
transcriptase domains (Saldhana et al . , 1993), so a DNA copy 
of the modified transcript might be made and incorporated 
during recombination into the host genome. This model for 
25 intron mobility rests on the observation that excised 
introns can reverse -splice into new RNA molecules, targeting 
sequences complementary to endogenous EBS sequences. 
Reverse splicing has been observed for both group I (Woodson 
and Cech, 1989) and group II introns (Augustin et al., 1990; 
30 MSrl and Schmelzer, 1990; Muller et al., 1991). This 
observation does not necessarily reflect the relative 
efficiency of forward and reverse reactions. After all, the 
first step of group I intron self -splicing does not proceed 
with sufficient efficiency for its products to affect the 
35 apparent endpoint of forward reactions. However, an 
important implication of reverse- splicing reactions is that 
the branching reaction of group II introns is capable of 
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reversal. The attack of 5'-exon on mutant lariats of intron 
bll was examined previously, and these studies suggested 
that the most difficult step in reverse-splicing was 
reversal of the branching reaction (Augustin et al., 1990). 
5 This may have been due to a lack of kinetic control on the 
system or the use of mutants rather than wild-type 
sequences. It is now well-established that many 

modifications that affect the second step of group II self- 
splicing also affect the first step (Chanfreau and Jacquier, 

10 1994) . Therefore, an important ramification of our analysis 
is that the final step necessary for reverse splicing, the 
reverse -branching reaction, is actually facile and that the 
intron can readily surmount this particular kinetic barrier 
to reintegration into an RNA site. Proteins that promote 

15 intron mobility, such as the maturases (Lambowitz and 
Bel fort, 1993), may act by simply inhibiting the forward 
second step of splicing, thereby exploiting the natural 
propensity of the branching reaction to proceed in reverse. 

20 Conclusions 

Together, these data show that the first step of splicing 
involves the interplay of three reactions, all of which have 
the capability to serve important biological functions. The 
hydrolysis reaction results in proper splice-site selection, 

25 and although it is a slow chemically limited process, it is 
able to drive the first step of splicing to completion. The 
two transesterification reactions, branching and reverse 
branching, are chemically rapid and partially limited in 
rate by slow conformational change. The forward branching 

3 0 reaction provides kinetic parameters indicating that 
tertiary interactions with D5 mediate association of D56 in 
the active site. The forward reaction cannot go to 
completion because it is a highly reversible transformation. 
The reverse reaction is very fast, proceeding through a 

35 higher initial population of conformational ly competent 
molecules. The efficiency of reverse branching would 
permit it to aid in proof-reading of 5' -splice-site 
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selection and may explain the facile mobility of group II 
introns. The presence of an efficient back-reaction during 
the first step of splicing necessitates a rapid second step 
in order to drive group II intron self -splicing to 
5 completion. 

ESB&IMEBSSL PBSCRJFTXON 

Plasmids, templates, and transcription of RNA: The group II 
intron examined in this study is intron By of the COX1 gene 
10 encoding the mitochondrial cytochrome oxidase subunit 1 of 
Saccharomyces cerevisiae (ai5g) . Plasmids pJDI3'-673 and 
pJDI5'-75 were kindly provided by Dr. Philip S. Perlman 
(Jarreli et al., 1988a). Plasmid pT7-D56 was constructed by 
Kunkel mutagenesis of pJDI5'-75 with the following 
15 modification to published protocol (Kunkel et al., 1991): 
Annealing of the mutagenic oligonucleotide was done by 
heating to 65 ;C for 5 minutes and then placing samples 
directly on ice. Mutagenic oligonucleotides were 

synthesized on an ABI 392 DNA/RNA synthesizer and purified 
20 by PAGE. Plasmid pT7-D56 contains domains 5 and 6 with a 
unique EcoRV site at the 3' end of D6 . Plasmid pT7-D56A- 
was similarly constructed except that the branch point 
adenosine has been deleted. Plasmid templates were 
linearized by restriction with EcoRV before transcription 
25 with T7 RNA polymerase (Davanloo et al . , 1984). 
Transcription of pJDI3'-673 and pJDI5'-75 to make exD123 and 
D5 RNA, respectively, has already been described (Pyle and 
Green, 1994). Transcriptions of pT7-D56 and pT7-D56A- to 
make D56 and D56A- RNA, respectively, were done under 
30 identical conditions. Synthetic 5'-exon analog 5'- 

CGUGGUGGGACAUUUUC-3' (SEQ I.D. No.: ) was synthesized 

from 2'-silyl protected RNA phosphoramidites on an ABI 392 
DNA/RNA synthesizer. This RNA was purified using a 
reverse -phase OPC (ABI) column followed by PAGE. 

35 

Forward reaction conditions: Reactions were performed in 
0.65 ml Slickseal tubes (National Scientific) in a final 
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volume of 10 or 20 ptl . RNA provided in excess was 
unlabeled, while RNA provided in trace was labeled with 32 P, 
as indicated. To begin a reaction, aliquots of exD123 and 
D56 RNA (stored in lOmM MOPS and lmM EDTA pH 6.5) were 
5 individually diluted and combined with 1M MOPS pH 7.0 to a 
final concentration of 55 mM MOPS pH 7.0. The RNA 
components were then heated to 95 iC for 1 minute to unfold 
alternative structures formed by storage at -20 jC (Groebe 
and Uhlenbeck, 1988; Walstrum and Uhlenbeck, 1990). RNA 

10 samples were cooled at room temperature for 15 seconds and 
then briefly spun in a microcentrifuge. Reaction was 
initiated by simultaneously combining the exDl23 sample, the 
D56 sample and a concentrated salt solution such that the 
final reaction contained both RNA components at indicated 

15 concentrations and a reaction buffer of 40 mM MOPS pH 7.0, 
100 mM MgCl 2 and 500 mM NJCl. Results were identical if 
reaction was initiated by combining D56 and exD123 RNAs that 
had been previously incubated individually with reaction 
buffer. After either initiation protocol, reactions were 

20 incubated at 45°C in a water bath, aliquots were removed at 
selected timepoints and combined with equal volume of gel 
loading buffer {10 M urea, 0.1% xylene cyanol and bromphenol 
blue dyes, 8.3% sucrose, 40 mM tris pH 7.5 and 0.83 mM EDTA, 
final concentrations) before placing on ice. Samples were 

25 then subjected to electrophoresis on denaturing 4% 
polyacrylamide gels, which were then dried and quantitated 
on a Packard Instant Imager using methods previously 
described {Pyle and Green, 1994). Experiments were 
performed to determine optimal reaction conditions 

3 0 (described above) for quantitation of branching kinetics: 
a) the pH was varied from 5.5 to 8.0, b) the temperature was 
varied from 37° to 50°C, and c) the salts tested were 100 mM 
MgCl 2 (99.99% pure, Aldrich) , MgS9 or MnC£ combined with 
500 mM NH<C1, <NH 4 ) 2 S0 4 , KC1 or NaCl. 

35 

Determination of Icinetic parameters for the forward 
reaction: The forward reactions fit a single exponential 
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with a plateau at 50% reaction completion. To quantitate 
the rates of these reactions, the data was corrected for a 
50% reaction endpoint and the natural logarithm of the 
quantity (1- fraction product) was plotted vs. time in a 
5 semilog format as shown in the inset to Figure 13. The 
equation of the line was solved for the reaction half-time, 
which was converted to k^ (the = 0.693 / t 1/2 ) . The same 
rate constant would also be obtained by fitting the raw data 
to a single exponential, but it would be harder to see 
10 deviations from reaction homogeneity at long reaction times, 
initial bursts of reaction were estimated from the y- 
intercept of semilog plots. The rate of the forward 
hydrolysis reaction was determined by guantitating D56A- 
cleavage of exDl23, using published procedures for studying 
15 D5-catalyzed reaction (Pyle and Green, 1994). The for 
the D56-exD123 complex and the overall k cat for the forward 
reaction were determined by plotting k ob , values as a 
function of [exD123] (Figure 13). The data were fit using 
a published equation describing 1:1 bimolecular association 
20 (Pyle and Green, 1994) . and k c . t values, with accompanying 
standard errors were determined by best fit to the equation. 

Branch isolation: Trace amounts of 32 P-labeled D56 (10 nM, 
labeled internally or at the 5 '-end) was incubated with 

25 excess exD123 (l^M) in 40 mM MOPS pH 7.0, lOOmM MgCl 2 , 500mM 
NH 4 C1 at 45 |C for 2-3 hours. After incubation, the sample 
was ethanol precipitated and then resuspended in 15^1 MOPS 
storage buffer (10 mM MOPS, 1 mM EDTA, pH 6.5) and 15 jil gel 
loading buffer (described above) before electrophoresis on 

30 a 4% denaturing polyacrylamide gel. The gel was transferred 
to a piece of X-ray film that had been developed and then 
washed thoroughly. The gel was analyzed for 5-10 minutes 
on a Packard Instant imager and the actual -size image was 
laser-printed on a transparency sheet. By aligning the 

35 printout with the scanner coordinates of the gel, bands were 
easily identified and were excised with precision. Excised 
bands were eluted into MOPS storage buffer overnight before 
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ethanol precipitation and resuspension. 

Characterization of branched product: In order to 

sequence the position of the branch site, phosphorothioate- 
5 doped D56 was transcribed using small amounts of specific Qf- 
thio A, C, U or 6 nucleotide triphosphates. The transcripts 
were 32 P-end-labeled, reacted with exD123 to form the 
branched species, purified and then isolated. Iodine 
cleavage of phosphorothioates in the branched molecules, 

10 followed by electrophoresis on a high resolution 
polyacrylamide gel resulted in base- specif ic sequencing 
ladders terminating at the bulged A branch-point in D6. A 
large gap separated this terminal cleavage from the branched 
material at the top of the gel. For analysis of the 

15 branched trinucleotide, PI nuclease was used to digest 32 P- 
D56/D123 (typically 10,000 to 20,000 cpm branched RNA} in a 
buffer containing 25 mM sodium acetate -pH 5.3, 0.025 pq/fil 
tRNA, and 0.005 units/>l PI nuclease (Boehringer) . 
Reactions were incubated for 5-15 min at 37 |C before adding 

20 to equal volumes of quench buffer (36% formamide, 0.1X TBE, 
1.8% sucrose, 0.02* xylene cyanol dye and 25 mM EDTA final 
concentrations) and loading on a denaturing 20% 
polyacrylamide gel. Following the separation of 

mononucleotide (pU) from branched trinucleotide pA (2' -pG) 3' - 

25 pU by electrophoresis, gels were dried and quantitated on a 
Packard Instant imager . 

Conditions for reverse reaction and kinetic analysis: 
Stocks of branched D56/D123 species (labeled either at the 
30 5' -end of D56 or internally in D56) and 5'-exon were diluted 
and reacted using the same procedures, buffers and reaction 
conditions described above for the forward reaction. The 
k^, values for the reverse reaction were obtained by fitting 
the fraction reacted vs. time to a biexponential equation: 

35 

y = Ae" k l l + (l-Aje'^ 
where A = the fraction of molecules that participates in the 
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fast phase of reaction, k x = rate of the fast phase (burst, 
possibly chemistry) , k 2 * rate of the slow phase 
(conformational change) and t = time along the x-axis. 

5 Example 3: Group XI intron ribozvmes that cleave DMA and y m 
linkages with Pimilar efficiency, and lack co ntacts with 
substrate 2'-hvdroxvl arming 

Summary; 

10 Group II introns are self -splicing RNAs that catalyze a 
variety of important biochemical transformations. These 
autocatalytic molecules can be reconfigured into highly 
specific, multiple -turnover ribozymes that cleave small 
oligonucleotides in- trans. Unlike other catalytic RNA 

15 molecules, the group II ribozymes cleave DNA linkages almost 
as readily as RNA linkages. A ribozyme variant containing 
intron Domain 3 cleaves DNA linkages with an efficiency 
comparable to that of restriction endonuclease EcoRI. By 
probing the effects of single deoxynucleotides along the 

20 substrate strand, the ribozymes were found to bind substrate 
without engaging 2'-hydroxyl groups for ground-state 
stabilization. The ribose 2'-hydroxyl group at the cleavage 
site plays a minor role in transition- state stabilization by 
these ribozymes, suggesting a catalytic strategy that 

25 broadens the scope of potential substrates. Because they 
also do not engage 2' -hydroxy 1 groups in substrate binding, 
group II intron ribozymes may rely exclusively on base- 
pairing for molecular recognition. 

30 Introduction: 

The ribose 2'-hydroxyl group is widely regarded as an 
important functionality for molecular recognition and 
chemical catalysis by RNA [Musier-Forsyth and Shimel, 1992; 
Pyle et al. 1992; Herschlag et al., 1993; Smith and Pace, 

35 1993; and Stobel and Cech, 1993]. Evidence for its 
participation as either a hydrogen bond donor or acceptor 
demonstrate that the 2 / -hydroxyl group is a multifunctional 
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recognition determinant [Pyle et al. 1992? Herschlag et al., 
1993; Smith and Pace, 1993; Yap and Musier-Forsyth, 1995; 
Pley, H.M. et al., 1994; and Pley, H. M. et al., 1994]. For 
ribozymes in the hammerhead mechanistic class (class A) , the 
5 2'-hydroxyl group at the substrate cleavage site behaves as 
a nucleophile in the phosphodi ester cleavage reaction [Pan 
et al., 1993] . For ribozymes in another mechanistic class 
(class B) , which includes RNase P, group I introns and group 
II introns, the 2'-hydroxyl group adjacent to the scissile 

10 bond is not the attacking nucleophile (Scheme 1) [Pyle, 
1993] . Nonetheless, in two examples (group I intron 
ribozymes and RNase P) , this 2' -hydroxy 1 at the cleavage 
site has been observed to play an important role in 
transition- state stabilization [Smith and Pace, 1993, 

15 Herschlag et al., 1993; and Perreault and Altman, 1992]. 
Neighboring 2' -hydroxy 1 groups have also been observed to 
contribute to substrate recognition and binding [Pyle et al. 
1992; Herschlag et al. 1993; Strobel and Cech, 1993 Pyle and 
Cech 1991; and Bevilacqua and Turner, 1991]. Although 2'- 

20 hydroxyls within a core region (Domain 5) of group II 
introns are essential for binding and catalysis, the 
specific role of 2'-hydroxyl groups on substrates for group 
II intron ribozymes has not been investigated. 

25 Through the process of directed molecular evolution, a 
group I intron ribozyme was evolved to cleave DNA linkages 
with relatively high efficiency [Robertson and Joyce, 1990; 
and Tsang and Joyce, 1994]. But to date, there have been 
no examples of "natural ribozymes" that cleave DNA with an 

30 efficiency comparable to their cleavage of RNA. The term 
•♦natural" is used rather loosely here since most ribozymes 
are man-made pieces, comprising the actives site of self- 
cleaving RNA molecules. However, the term is used here to 
differentiate wild-type active-site sequences from those 

3 5 that have been optimized by the process of directed 
molecular evolution or selection. Regardless of their 
derivation, ribozymes that make little use of substrate 2'- 
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hydroxyl groups are an important source of information about 
alternative strategies in the catalytic repertoire of 
nucleic acid enzymes. In addition, this behavior by natural 
ribozymes may point to a biological impact of ribozymes on 
5 genomic DNA. The data presented here indicate that certain 
ribozymes derived from a group II intron stimulate efficient 
catalysis without engaging the 2'-hydroxyl groups on their 
target substrates. 



10 Self -splicing group II introns have been found in the 
organellar genes of plants, lower eukaryotes and also in 
prokaryotes [Pyle, 1995; and Michel and Perat, 1995] . The 
excision of ubiquitous group II introns is essential for 
metabolism in many of these of organisms. Group II 

15 introns can be arranged into six domains, of which domain 1 
(Di) and domain 5 <D5) have been found to be essential for 
catalytic activity [Koch et al., 1992; and Michels and Pyle, 
1995] . Other domains are dispensable, although certain 
deletions result in reduction of particular activities 
20 [Pyle, 1995] . Self -splicing can proceed via either one of 
two pathways. One pathway involves initial 5' splice site 
attack by the 2'-hydroxyl group from the bulged adenosine of 
domain 6 (D6) , which generates a "lariat" intron that is 
ultimately released after the second step of splicing 
25 [Peebles et al., 1986; and Schmelzer and Schweyen, 1986; and 
van der Veen et al., 1986]. Alternatively, water can act as 
the nucleophile in the first step, resulting in ligated 
exons and release of linear intron [Jacquier and Rosbash, 
1986; van der Veen et al., 1987; Jarrell et al., 1988b; and 
30 Daniels et al., 1995]. A recent study has demonstrated that 
hydrolytic attack is a legitimate and common pathway for in- 
vitro self -splicing of group II introns [Daniels et al., 
1995] . 



in spite of their biological importance, there have been 
relatively few investigations of the structural features and 
mechanistic enzymology of group II introns [Jarrell et al., 
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1988; Franzen et al., 1993; Pyle and Green, 1994; Padgett et 
al., 1994; and Podar et al., 1995]. Recently, a detailed 
kinetic study on a ribozyme derived from the group II intron 
ai5g was reported [Michels and Pyle, 1995]. That work 
5 established a minimal group II ribozyme system in which Dl, 
acting as an enzyme along with D5 as cofactor, cleaves short 
RNA oligonucleotides analogous in sequence to the 5' splice 
junction (Scheme 2) . Similar to hydrolytic splicing 
reactions (in either forward or reverse splicing), the 

10 nucleophile is water or hydroxide (Scheme 1) . As expected 
for base pairing between the intron binding sites (IBS1 and 
2) and exon binding sites (EBS1 and 2), Dl binds to its 
substrate strongly with a X* of 6.3 nM. The cleavage reaction 
is rate limited by chemistry, with k CBt = k chem = - 0.03 min" 1 

15 at saturation. Using this kinetic framework, it was 
possible to incorporate single -atom changes into the 
substrate and probe the reaction mechanism. 

In this investigation, the role of 2' hydroxyl groups on the 

20 oligonucleotide substrate for a group II ribozyme was 
examined by incorporating deoxynucleotides into specific 
positions along the strand. By monitoring the effects of 
deoxynucleotide substituents on catalysis and binding to the 
ribozyme, it was possible to evaluate the role of these 2' 

25 hydroxyl groups. The data indicate that, for two different 
ribozyme constructs, 2' hydroxyl groups along the substrate 
backbone are minimally involved in tertiary contacts that 
stabilize binding to the ribozyme or stimulate catalysis 
during the chemical step. In one case, deoxynucleotide 

30 linkages are cleaved with the highest efficiency (k cat /KJ 
observed for any ribozyme to date, having a k cat comparable 
to that of certain protein restriction enzymes. Although 
significant effects are observed for deoxynucleotides at the 
cleavage site and immediately adjacent, these effects are 

35 small compared to those observed for other ribozymes in 
mechanistic class B, which show ~l0 3 -fold effects on k^ 
[Herschlag et al., 1993; Smith and Pace, 1993; and Herschlag 
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et al., 1993]. Trivial explanations for the effects 
reported here are ruled out by kinetic experiments that 
reveal, for both ribose and deoxyribose substrates, that the 
reactions proceed with the same rate-limiting step. These 
5 data suggest that, for certain types of reactions, group II 
intron ribozymes may utilize a novel catalytic strategy that 
bypasses nearby 2' hydroxyl groups and exploits other 
factors for transition- state stabilization. 

10 Experimental Pr oe edu T?ff? 

Preparation of RNA: RNA transcriptions of Dl (425 nt) and D5 
{58 nt) , from plasmids pT7Dl and pJDI5'-75 respectively, 
were carried out as previously described [Michels and Pyle, 
1995) . Ribozyme D123 (710 nts) was transcribed from 

15 plasmid pT7D123, which was cut with Hind III. This plasmid 
was constructed by Kunkel mutagenesis [Kunkel et al . , 1991] 
from plasmid pT7I3'-675 [Jarrell et al., 1988; and Pyle and 
Greene, 1994]. Mutagenesis removed the 5'-exon from the 
parent plasmid, such that the first nucleotide of the 
20 resulting transcript corresponds to the first nucleotide of 
the ai5g group II intron from yeast mitochondria. Plasmids 
pJDI5'-75 and pJDI3'-675 were used. Oligonucleotide 
substrates, all ribose (rS) or with deoxyribonucleotides 
(chimeric substrates) selectively placed at specific 
25 positions, were synthesized on an ABI 392 DNA/RNA 
synthesizer, and worked up according to standard procedures 
[Scaringe et al . , 1990]. Oligonucleotides were composed of 
the same base sequence which, for the all ribose substrate 

(rS) was S'CgjMSM^QmEm^SZSCGGU (SEQ I.D. No. ). 

30 Underlined positions represent IBS2 and IBS1 respectively, 
dashed underline represents the intronic sequences and w 
indicates the cleavage site. The d(IBSl) block substrate was 

5 ' CGUGGUGGGAdCdAdUdUdUdUC-GAarnan f fiEn i.d. No. ) in 

which all IBS1 residues other than the nucleotide at the 
35 cleavage site were deoxynucleotides (bold). The d(IBS2) 
substrate contained 5 ' CdGdUdQdQdUdGGGACAUUUUC w GAGCGGU (SEQ. 
10 No - >• The d(link) substrate contained 
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deoxynucleotide substitutions at the three positions between 
IBS1 and IBS2 (as above). The d(intron) substrate contained 
the residues S f CGilS^O^cmM^ dq^AdOd^dQdGdXf . (SEQ I.D. 
No. _J For substrates containing single deoxynucleotide 
5 substituents, they are named according to the specific 
position relative to former exon and intron sequences. For 
example, -ldC is 5 ' rgmguGGGACAIJUUUdC v GAGCGGU, (SEQ I.D. No. 

) and -2dU represents 5 ' CGUGGUGG GACAUUUdUC v GAGCGGU , (SEQ 

I.D. No. ) etc. 

10 

Reaction conditions for kinetic analyses: Characterization 
of preliminary reaction rates for the d(IBS2), d(link), 
d(IBSl) and d(intron) substrates was carried out under 
single turnover conditions with saturating ribozyme. The 

15 reaction buffer contained 1 M KC1, 100 mM MgCl 2 , and 80 mM 
3- [N-Morpholino]propanesulfonic acid (MOPS) pH 7.5. 
Reactions were usually carried out in 20 ul volumes at 42QC, 
with excess ribozyme (100 nM - 250 nM Dl or D123, 3 uM - 6 
uM D5, as indicated) and trace 5' 52 P labeled substrates 

20 (0.25 - 1 nM) . Initiation of the reaction, measurement of 
v values and single turnover kinetic analysis were 
conducted as described previously (Michels and Pyle, 1995] 
with one exception: Preliminary results from fluorescence 
experiments showed that Dl can slowly fold and bind to the 

2 5 substrate without D5 . To ensure that both rS and chimeric 
substrates reached binding equilibrium with ribozyme before 
reaction, substrates containing single deoxynucleotides were 
preincubated for 20 minutes with Dl or D123 prior to the 
initiation of the cleavage reaction by addition of D5. 

30 Appropriate amounts of substrate and Dl were mixed into 4 ul 
of B0 mM MOPS pH 7.5, and heated at 95 °C for 1 minute to 
denature potentially misfolded structure. After cooling to 
42°C, 8 ul concentrated salt in 80 mM MOPS pH 7.5 was added 
to this Dl/S mix to give a final concentration of 1M KC1 and 

35 100 mM MgCl 2 . The solution was then incubated at 42°C for 20 
minutes. The reaction was initiated by the addition of 8 ul 
D5 mixture, preincubated in parallel under the same salt 
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concentration. This additional procedure was merely a 
precaution, however, since control experiments showed that 
for rS and mutant substrates, preincubation of Dl or D123 
with substrate results in the same kinetic parameters as 
5 analyses without preincubation, performed as described 
previously [Michels and Pyle, 1995] . 

Some of the ribozymes described herin are capable of 
catalytic activity at physiological conditions. Kinetic 

10 parameters of a ribozyme which contains domain 1 and domain 
3 (D13) is compared with the kinetic parameters of a 
ribozyme that contains domain 1, domain 2 and domain 3 
(D123) (Figures 25C and 25D) . Salt concentrations used are 
the following: lOmM MgCl 2 , 25mM MgCl, and 50 mM MgCl 2 . These 

15 low salt conditions are used in the absence of any other 
monovalent cations such as ammonium or potassium. It is 
common to include such monovalent cations at molar 
concentrations in reaction mixtures in order to effect 
catalytic activity. However, these compounds do not require 

20 such salt concentrations. The rate for the D13 ribozyme at 
lOmM MgCl 2 is 9.5 min, at 25 mM MgC^ is 22.9 min and at 50 
mM MgCl 2 is 518 min (Fig. 25C) . The rate for D123 ribozyme 
at 25 mM MgCl 2 is 12 min and at 10 mM MgCl 2 is 54 min. Thus, 
these compounds are capable of functioning catalytically at 

25 physiological conditions. 

Determination of individual kinetic parameters: In order to 
determine ^ and k^,. for reaction of substrates with the Dl 
ribozyme, k^, was monitored at varying concentrations of Dl, 

30 while holding [D5] at saturation. The data fit an equation 
describing bimolecular association, from which the kinetic 
parameters and £ t were extracted. The data were 
quantitated and analyzed as described previously [Michels 
and Pyle, 1995; and Pyle and Green, 1994]. In order to 

35 determine if, in each case, k Cit represented / it was 
necessary to conduct parallel pH/rate profiles of the 
reaction. Values of k cat for rS, -IdC and -2dU substrates 
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were compared at pH 6.5, 7.0 and 7.5, in 80 mM MOPS, 100 mM 
MgCl 2 and 2 M KC1. MOPS was chosen because it maintains 
high buffering capacity throughout this pH range. In 
addition, MOPS maintains a stable pH over a wide temperature 
5 range. For these experiments, the KC1 concentration was 
raised to 2 M in order to accurately quantitate the slow 
rate of -ldC cleavage at pH 6.5. Doubling the KCl 
concentration speeds up the reaction by approximately 10 
fold for all substrates, at all pH values. With the 
10 exception of -ldC cleavage at pH 6.5, which is too slow to 
measure, the same pH/rate behavior is seen for all three 
substrates (with lower precision) at 1 M KCl. 

Mapping of Cleavage Products and Chemical analysis of 3'- 

15 ends: Partial alkaline degradation and enzymatic 

digestions with nucleases Tl, B. Cereus (Pharmacia), and PI 
(Boehringer) were carried out as described [Michels and 
Pyle, 1995] . In order to examine whether cleavage of the - 
ldC substrate was indeed occurring at the deoxyribose 

20 cleavage site, 5' "P labeled substrates, rS and -ldC, were 
cleaved to 50% completion, under conditions described above. 
The excess salt was removed with oligo purification 
cartridges (OPC, from Applied Biosystem, Inc.). The 
reaction products were then oxidized with 0.1 M NaIO< in 0.1 

25 M sodium acetate pH 5.0 for 90 minutes at room temperature 
in the dark (Stephenson and Zamecnik, 1967] . Excess NaI0« 
was removed using an OPC cartridge. The oxidized mixture was 
then reacted with 0.1 M Lucifer Yellow hydrazide (Lucifer 
Yellow CH, Aldrich) in 0.1 M sodium acetate pH 5.0 at room 

30 temperature in dark for approximately four hours. If the 
3' -terminus of an oligonucleotide is a ribose, reaction with 
NaI0 4 will oxidize the terminal sugar to a dialdehyde that 
reacts readily with hydrazine dyes [Stephenson and Zamecnik, 
1967; and Hansske et al., 1974], such as Lucifer Yellow. 

35 Any 3 '-terminal deoxyriboses are resistant to attack of 
NaI0 4 . After passing the dye reaction through an OPC 
cartridge to remove excess Lucifer Yellow, the mixture was 
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combined with denaturing dye, and loaded onto a 20% 
denaturing polyacrylamide gel. Attachment of Lucifer Yellow 
to the 3' -terminus of an oligonucleotide substantially 
reduces its mobility in a denaturing gel; the identity of 
5 3 '-terminal sugars was then analyzed by autoradiography of 
the products after electrophoresis. 

RESULTS: 

Preliminary Characterization oj Deoxvnucle otide Rat* 

10 Effectg; 

Effects of block substitution on reaction rate: In order to 
locate regions of important 2' hydroxyl groups, chimeric 
substrates were synthesized in which the cleavage site (-1C) 
was always a ribose and blocks of surrounding residues were 

15 substituted with deoxynucleotides. The blocks of 

deoxynucleotides were clustered into four functional 
regions: IBS2 (-16 to -11), link {-10 to -8), IBS1 (-7 to - 
2) and intron ( + 1 to +7) (Scheme 2). The IBS2 and IBS1 are 
stretches of nucleotides that pair with complementary 

20 sequences on the Dl ribozyme [Michel and Ferat, 1995; and 
Michel et al., 1989]. The "link" region is a three 
nucleotide spacer between these paired sequences and the 
"intron" region is composed of seven nucleotides 3 '-to the 
ribozyme cleavage site (see experimental section) . The 

25 chimeric substrates were tested using single turnover 
cleavage reactions under saturating ribozyme conditions (100 
nM Dl, 3 uM D5) , and cleavage rates were compared to those 
of an all-ribose substrate (rS) (Table 3, column 1) . 

30 
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Table 3. Effects of block deoxyribose substitutions on 
ribozyme rate. 



10 





l a 


2 


3 


Modified 


Effect of dN 


250 nM Dl b 


SfM D5 b 


Substrate 


{rel rate) 1 


(rel rate) 1 


(rel rate)' 1 










dUBS2) 


3.2±0.51 


1.4±0.22 


2.9±0.46 


D(lBSl) 


inactive 


inactive 


inactive 


d(link) 


1.4±0.22 


ND 


ND 


d(intron) 


3.3±0.53 


3.3±0.53 


3 .1±0.50 
1 — rrns — -J — lJ 



limit for saturation of the unmodified rS substrate. See 
Table 4 for of rS. Relative rates are for the 

chimeric substrates divided by that of rS. Variance 
15 calculated to 95% confidence, from multiple trials on a 
single derivative. ND, not determined. 

b Rate relative to rS cleavage under equivalent conditions. 
Except for d(IBSl), the block-deoxynucleotide substrates 
showed small rate decreases relative to rS, indicating that 

20 the 2' hydroxyl groups within IBS2, the linker and intronic 
nucleotide sections of the substrate are not essential for 
catalytic activity. Within experimental error, the d(link) 
substrate was as active as rS, so it was excluded from 
further investigation. The rate of d{IBS2) cleavage was - 3 

25 fold slower than rS. That this effect was due to well- 
documented decreases in hybrid stability [Bevilacqua and 
Turner, 1991; Martin and Tinoco, 1980; and Hall and 
McLaughlin, 1991] was supported by restoration of d(IBS2) 
cleavage rates at higher concentrations of enzyme components 

30 (Table 3, column 2). Conversely, the three-fold difference 
in cleavage rates between rS and d(intron) substrates 
remained the same at higher concentrations of both Dl and D5 
(Table 3, columns 2 and 3), suggesting that the decrease of 
cleavage rate in d(intron) is due to effects on the chemical 

35 rate or some effect other than decrease in binding affinity 
to either component of the ribozyme. 



WO 96/22689 



PCT/US96/01337 



-93- 

Under all conditions, d(IBSl) was inactive within the 
experimental window. Given that this oligonucleotide has a 
ribose at the cleavage site, there are at least two 
explanations for its lack of activity: Either important 2'- 
5 hydroxyl groups lie within the IBS1 region, or a DNA/RNA 
hybrid duplex adjacent to the cleavage site is destabilizing 
for other reasons. The latter possibility is reasonable 
given that DNA/RNA duplexes are thermodynamically 
destabilized [Hall and McLaughlin, 19911 , their 

1C microstructural features can differ slightly from RNA/RNA 
duplexes [Egli et al . , 1993] and they can deviate 
significantly from A-form geometry and pitch, an extreme 
case being the recent crystal structure of a B-form DNA/RNA 
helix [Chen et al., 1995]. But without a signal for 

15 activity, it was impossible to differentiate these models 
and determine which, if any 2' hydroxyl groups of IBS1 were 
important in ground state or transition- state stabilization. 

Effects of single deoxynucleotide substitutions within IBS1: 
2 0 To further investigate the importance of 2' hydroxyl groups 
within IBS1 and at the cleavage site, rS was substituted 
individually with single deoxynucleotides, and the mono- 
substituted chimeras were tested using single turnover 
reactions at saturating ribozyme conditions (Figure 16) . 
25 The slowest cleavage rate is observed for a substrate 
containing a deoxyribose at the cleavage site <-ldC) . An 
observed rate of 0.0013 min" 1 indicates that the 2 f hydroxyl 
group at the cleavage site is the most important of all 2' 
hydroxyl groups along the substrate. However, this is only 
30 a 16 fold decrease relative to the rate of rS cleavage. 

Cleavage of the -2dU substrate proceeded with a k^ of 
0.0042 min 1 , which represents a 5- fold decrease in rate 
relative to rS. Within experimental error, single deoxy 
35 substitutions at all other positions within IBS1 (-3 to -7) 
had only minor effects on rate (Fig. 16). However, as 
observed for a block deoxynucleotide substitutions of IBS2, 
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when the -2 to -7 IBS1 nucleotides are simultaneously 
converted to deoxy ribonucleotides, there are likely to be 
additional effects due to duplex weakening that raise the 
Dl-S K„. Taken together, the single -deoxynucleotide 

5 investigations reveal that the most important 2'-hydroxyl 
groups along the substrate are the 2' -hydroxy 1 at the 
cleavage site (-1C) and one that is adjacent to it (-2U) in 
a region that is normally considered the 5'-exon sequence. 
Under the single -turnover, saturating ribozyme conditions 

10 used in these experiments, the rates were likely to 
represent effects on the chemical step of the reaction. 
This was not definitively concluded until further analysis 
was performed, establishing the rate-limiting step for k obI 
in each case. Before that could be accomplished, however, 

15 it was important to establish that cleavage was taking place 
at the proper site on the chimeric substrates. 

Ph*rarferization of Reaction Products: 

Mapping the cleavage site of the -ldC substrate: One reason 
2C that cleavage of a deoxynucleotide linkage might appear to 
be efficient is that the ribozyme might "skip" the modified 
linkage and choose an adjacent ribose phosphodiester moiety 
for attack [Kahle et al. r 1993]. As a first approach to 
investigating that possibility, the cleavage products from 
2 5 the D1/D5 ribozyme reactions were mapped using nucleases and 
alkaline hydrolysis, followed by high-resolution sequencing 
(Fig. 17A) [Kuchino and Nishimura, 1989] . Reaction of -ldC 
with alkali (lane 11) results in a blank space at the 
position corresponding to mobility of reaction product {a 
30 17-mer, lane 9) . 

To more precisely evaluate the mobility of the cleavage 
products, the substrates and products were mapped by partial 
PI nuclease digestion, which cleaves both DNA and RNA 
35 linkages, resulting in fragments containing the same 3' 
hydroxyl and 5' phosphate termini observed upon reaction 
with the ribozyme [Kuchino and Nishimura, 1989] . The 
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products of ribozyme cleavage of both the rS (lane 6) and 
the -ldC (lane 9) substrates were found to migrate at the 
identical position as the 17-mer generated by partial PI 
nuclease digestion of each substrate (lanes 7 and 8, 
5 respectively) . Taken together, the mapping studies suggest 
that cleavage of the -ldC substrate occurs at the proper 
sequence and is not the result of attack at an adjacent 
ribose linkage. 

10 High -resolution mapping of reaction products can be 
complicated by many factors. It is important to directly 
probe the terminus of the cleavage product, since correct 
cleavage of the -ldC substrate will generate a product with 
a 3' -deoxy terminus. NaI0 4 oxidation of an oligonucleotide 
15 with a 3' ribose terminus will convert the ribose sugar ring 
to a dialdehyde, which can be labeled with a hydrazine 
reagent [Hansske et al., 1974], such as Lucifer Yellow CH 
hydrazide. The labeled oligonucleotide can be distinguished 
from the unlabeled one because it migrates more slowly on 
20 denaturing gels. An oligonucleotide product with a 3' deoxy 
terminus would be resistant to oxidation/dye labeling and 
would not show a mobility shift. In this experiment, 
oxidation and labeling was done on a 50:50 mixture of 
substrate and product. The precursor oligonucleotide served 
25 as an internal control, since this component always 
possesses a 3' ribose terminus and should show a mobility 
shift if the above scheme is correct. The results of this 
analysis are shown in Figure 17B. For the rS substrate, 
after oxidation and labeling, both substrate and product 
30 shifted up (lanes 3, 4), indicating that both of them have 
3' ribose termini. For the -ldC substrate, the substrate 
mobility shifted after oxidation and labeling, while the 
product mobility was unchanged (lanes 6,5] . This indicates 
that the cleavage product of -ldC has a 3' deoxynucleotide 
35 terminus and that the -ldC substrate is cleaved at the 
correct site. 
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».^™ t *i~ analysis of Bsaaamaslsatids Bate Effect*: 

pH/rate analysis of the chimeric substrates: The kinetic 
effects of deoxynucleotide substitution could result from a 
number of factors: including effects on ribozyme 
5 conformation, binding or on the chemical step of catalysis. 
Furthermore, there have been several reports that 
deoxynucleotide substitution can change the rate-limiting 
step of reaction [Smith and Pace, 1993,- and Herschlag and 
Cech, 1990] , making comparisons of uncharacterized k,^ or k cit 

10 values relatively difficult to interpret. Therefore, it was 
essential to investigate whether k*. for cleavage of rS and 
-ldC under saturating Dl/DS conditions was governed by the 
same rate- limiting step. That cleavage of rS is limited by 
chemistry under these conditions was established previously 

15 using a comprehensive kinetic analysis, together with 
observation of a predicted "thio-ef f ect » on the rate of 
reaction. As in other hydrolytic cleavages of group II 
intron substrates, this k chTO is expected to vary 
logarithmically with pH [Pyle and Green, 1994] . 

20 

values of k^, for the rS, -ldC and -2dU substrates were 
obtained at several pH values using the same buffer (Fig. 
18) . In the pH range tested (pH 6.5 to pH 7.5), both rS and 
chimeric substrates have a remarkably linear pH/rate 

25 profile, in which the rate changes logarithmically with each 
pH unit. The slope of each pH/rate plot is very close to l, 
being 0.87 for rS, 0.94 for -ldC, and 0.85 for -2dU. A 
linear pH/rate profile with slope of one is regarded as 
evidence that chemistry is rate-limiting [Herschlag et al., 

30 1993; Pyle and Green, 1994; Dahm et al., 1993; and Khosla 
and Herschlag, 1994] and that some form of general base 
catalysis is occuring [Fersht, 1985] . The fact that all 
three substrates elicit this type of behavior suggests that 
k ch „„ is limiting in each case and that their cleavage rates 

35 can be directly compared. 



Effects of deoxy nucleotide substitution on the binding 



WO 96/22689 



PCT/US96/0I337 



-97- 

between substrate and ribozyme: The previous pH/rate 
experiment indicates that 2' -hydroxy 1 groups at -1C and -2U 
directly affect the chemical rate of rS cleavage. Although 
transition-state effects are evident, this does not rule out 
5 the possibility that ground state binding might also be 
af-fected by the deoxynucleotide substituents, since the 
large excess in ribozyme concentration will mask, to a 
certain degree, small changes in binding affinity between 
chimeric substrates and the ribozyme. To evaluate effects on 

10 binding, the rS, -ldC and -2dU substrates were cleaved at 
varying Dl concentrations. The k^ and values for 
cleavage of each substrate were calculated from the 
dependence of k^, on Dl concentration, as described 
previously [Michels and Pyle, 1995] . The results are shown 

15 in Fig. 19 and Table 4. The kc, t and K,, values for rS are in 
close agreement with the previous published results (Fig. 
19A) [Michels and Pyle, 1995]. 



Table 4. Kinetic parameters for D1/D5 ribozyme cleavage of 
20 the -ldC and -2dU chimeric oligonucleotide substrates. 



25 





k b 






Substrate* 


U0* 2 rain' 1 ) 


(nM) 


(10* M" 1 min 1 ) 










rS(all- 
RNA) 


2.5±0.070 


4.9±0.66 


5.1±0.83 


-ldC 


0.19±0. 0060 


14±1.7 


0.14±0.021 


-2dU 


0.51±0.029 


10±2.4 


0.51±0.15 



rate are shown. 

30 kinetic parameters and standard errors are derived from fits 
of the data shown in Figure 19. Kinetic data on the other 
chimeric substrates is presented in Figures 1,6 and Table 3. 
c Kat/^n errors were calculated by propagation of the # and 



errors . 



35 
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As observed in the previous experiments, the -IdC substrate 
has a 13- fold decrease on the chemical step (k c «) , but also 
has a 3-fold decrease on the ground state binding (K J , 
resulting in a 36 fold effect on k Clt / K , (Fig. MB. Table 
5 4) The -2dU substrate has a 5-fold decrease in the 
chemical step and a 2- fold decrease in the ground state 
binding (Fig. 19C, Table 4), which may be attributed to 
duplex destabilization caused by deoxynucleotide 
substitution. The k eat / K . for the -2dU substrate is 
10 therefore 10-fold lower than that of rS. For both the 
chimeric substrates, the small decreases in binding affinity 
(K J observed here confirm that previous rate comparisons 
were performed under conditions in which the ribozyme 
components were indeed saturating. This is also supported by 
15 the fact that, within experimental error, the decrease in 
between chimeric and rS substrates (13 fold for -IdC, 5 
fold for -2dTJ, respectively) agrees with the decrease (16 
fold for -IdC. 5 fold for -2dU) in the apparent rates 
obtained during the previous kinetic comparisons (Fig 16) . 



20 



aaaiggffl aa safes sffssfca for a faaf <ttoup II ribozyme . 
One can modify the minimal Dl ribozyme construct by adding 
back intron domains 2 and 3 in-cis (D123) , to generate a 
ribozyme that cleaves RNA with high efficiency. Under 
25 standard reaction conditions in which D123 (100 nM) and D5 
(3 M M) are in excess, this D123 ribozyme cleaves the rS 
substrate with an apparent k^, of 0.88 min" 1 (Fig. 20). 
This is 33 -fold faster than the Dl ribozyme under the same 
conditions, and comparable to the rate of other ribozyme 
30 systems. D123 cleaves the -IdC substrate with a k ob . of 
0.17 min" 1 (Fig. 20), which is only five-fold slower than 
cleavage of the ribose linkage and 90 -fold faster than 
deoxynucleotide cleavage by Dl . In both cases, cleavage 
occurs at the proper site (Fig. 17C) . Therefore, faster 
35 group II intron ribozymes cleave both DNA and RNA linkages 
with higher efficiency. 
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There is a smaller experimental window for detecting Dl 
reactivity, which is radically perturbed upon introduction 
of a hybrid EBS1/IBS1 duplex immediately adjacent to a 
ribose cleavage site (as in dIBSl, Table 2) . Therefore , it 
5 was not possible to detect Dl ribozyme cleavage of an all- 
DNA oligonucleotide (dS) . By contrast, the D123 ribozyme 
cleaves dS with a second-order rate constant of 2.6 ± 1.4 x 
10 3 M 1 min" 1 under standard reaction conditions using a 
constant [D5] of 3 /aM. Together with effects from single 
10 deoxynucleotide substitutions, this result indicates that 
any barrier to cleavage of entirely DNA molecules is likely 
to be due to global differences between RNA/RNA and DNA/RNA 
duplexes and not to a specific role for 2' -hydroxy 1 groups 
in the ribozyme mechanism. 

15 

The apparent five- fold preference of D123 for ribose over 
deoxynucleotide linkages compared to the sixteen-fold 
preference observed with Dl does not actually mean that D123 
is less capable of differentiating a ribose from a deoxy 
20 linkage. The ratio between the two cleavage rates by D123 
is smaller because they have different rate-limiting steps. 

The cleavage rate of -ldC by D123 (0.17 min" 1 ) is limited by 
the rate of chemistry, as demonstrated through pH/rate and 
other mechanistic analyses. However, the rate of rS 
25 cleavage by D123 (0.88 min") is not limited by the chemical 
rate of catalysis. Unlike all of the other processes 
reported in this paper, this rate is not pH dependent and it 
has a similar magnitude to rates of folding processes 
observed in stopped-flow fluorescence experiments performed 
30 under identical reaction conditions. Preliminary 
estimates for the chemical rate of rS cleavage by D123 are 
approximately 3 min* 1 , which would indicate that D123 has a 
similar -15-fold preference for ribose observed for the Dl 
ribozyme. Taken together, -these results suggest that 
35 significant transition- state stabilization can be achieved 
by fast group II intron ribozymes without engaging 2'- 
hydroxyl groups for high levels of transition-state 
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stabilization. It is therefore clear that ribozymes in 
mechanistic class B do not necessarily exploit the 
stabilizing effects of the 2'-hydroxyl group at the cleavage 
site. 

5 

Discussion 

At this time, we have only a modest understanding of the 
means by which RNA behaves as an enzyme or potentiates 
chemical events. There are very few natural ribozymes in 

10 existence, and only three in the mechanistic class B to 
which the group II intron belongs [Pyle, 1993; and Cech, 
1993] . For this reason, each ribozyme is an important 
resource, providing a distinctive window on mechanisms for 
transition-state stabilization by RNA. Group II intron 

15 active sites are particularly fascinating because they 
catalyze a variety of sophisticated reactions, including 
forward branching [Peebles et al., 19B6; Schmelzer and 
Schweyen, 1986; and van der Veen et al., 1986], reverse 
branching [M6rl and Schmelzer, 1990; Augustin et al., 1990; 

20 and Chin and Pyle, 1995], site-specific hydrolysis [Michels 
and Pyle, 1995? Jarrell et al., 1988A; Jarrell et al., 
1988B; Pyle and Green, 1994; and Podar et al . , 1995], 
polymerization [MOller et al., 1993], editing [Muller et 
al., 1993] and other reactions on RNA and DNA [M6rl et al. 

25 1992] . To understand more about the basis for catalysis by 
group II intron active-sites, we constructed a series of 
group II intron ribozymes that cleave RNA and DNA 
oligonucleotides in a manner that reflects factors essential 
for RNA folding and transition- state stabilization. In the 

3 0 most minimal system, a group II intron ribozyme was 
constructed of three pieces in which Dl behaves as the 
ribozyme, D5 as a catalytic cofactor, and a short 

oligonucleotide complementary to the EBS sequences acts as 
the substrate (Scheme 2) . This ribozyme catalyzes a single- 
35 step reaction, cleaving the oligonucleotide substrate into 
two pieces using water or activated hydroxide as the 
nucleophile. Previous studies have delineated a kinetic 
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framework for activity by this ribozyme and identified a 
rate-limiting chemical step of reaction [Michels and Pyle, 
1995] . This accessibility of its chemical rate constant 
renders this ribozyme useful for studies of transition-state 
5 stabilization. Because the D1/D5 ribozyme cleaves short 
oligonucleotides, these substrates can be generated on an 
automated DNA/RNA synthesizer, making it possible to 
incorporate unnatural nucleotides and single-atom changes 
for probing the chemical mechanism. There are two major 

10 results from these structure/function studies: First of 
all, the data indicate that a 2'-hydroxyl group at the 
cleavage site has only a small role in lowering the 
activation barrier for reaction. This effect is observed 
also for fast variants of the group II intron ribozyme. 

15 Second, the 2'-hydroxyl groups along the substrate backbone 
do not form an extensive network of sequence -independent 
backbone tertiary interactions upon binding to the ribozyme. 

The chemical role of 2' hydroxyl groups on the substrate: 

20 Previous studies have examined the role of ribose sugars at 
the site of cleavage by other ribozymes in mechanistic class 
B [Herschlag et al., 1993; Smith and Pace, 1993; Herschlag 
et al., 1993; Perreault and Altman, 1992; Robertson and 
Joyce, 1990; Herschlag and Cech f 1990A; Herschlag and Cech, 

25 1990B; and Kleineidam et al., 1993]. These studies have 
shown that a 2 '-hydroxyl group adjacent to the scissile bond 
is important in promoting catalysis. Single deoxynucleotide 
substitution at the cleavage site of the re trahyznena 
ribozyme (derived from a group I intron) resulted in a 590- 

30 fold decrease in the chemical rate of reaction [Herschlag et 
al., 1993A; and Herschlag et al., 1993B] . Analagous 
studies on the RNA component of RNase P resulted in a 
relative chemical rate difference of 3400 fold [Smith and 
Pace, 1993] . It was therefore surprising that, in a family 

35 of group II intron ribozymes, the absence of a 2' hydroxyl 
group at the cleavage site caused only a 16 fold reduction 
in chemical rate (Fig. 21) . This 16-fold differential 
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between ribose and deoxynucleotide substitution is 
consistent with the magnitude of an "induced effect" 
(Herschlag et al., 1993A; Herschlag et al., 1993B; and 
Exner, 197B] . in which the electron withdrawing property of 
5 the 2' hydroxyl group plays a role in stabilizing the 
adjacent 3' oxyanion leaving group that develops in the 
transition state. This effect stems from the fact that the 
-OH is a better electron withdrawing group than the -H. In 
a reaction in which negative charge develops on a nearby 
10 atom (3' -oxygen), the -OH group will better stabilize the 
transition state than the -H group, thus speeding up the 
reaction. 

Although unexpected, there is no fundamental principle 

15 against a limited role of a 2' hydroxyl group at the 
cleavage site. Since the 2' hydroxyl group does not 
participate as nucleophile or leaving group, the chemical 
mechanism of ribozymes in mechanistic class B does not 
implicitly require its presence (Scheme 1) IPyle. 19931. 

20 The closest example of this type of behavior was observed in 
studies of pre-mRNA splicing, a fact that is significant 
because group II introns are widely viewed as possible 
evolutionary predecessors of the eukaryotic spliceosome. In 
that work, it was reported that deoxynucleotide substitution 

25 at the 5' splice site had no effect on either step of 
splicing [Moore and Sharp, 1992] . While the chemical step 
of catalysis is unlikely to be rate-limiting in that complex 
system, it may be significant that there was no apparent 
block to reactivity, especially considering that base 

30 mutations at that site have strong effects on the apparent 
efficiency of splicing (Newman et al., 1985; and Parker and 
Siliciano, 1993). There is also precedent for group II 
intron-catalyzed attack on DNA linkages. Although it was 
reported to be very inefficient, the free bll intron (a 
35 different group II intron from yeast mitochondria) was 
observed to promote cleavage of a DNA substrate [M6rl et 
al., 1992] . 
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Small rate effects were observed for the 2' hydroxyl group 
located one nucleotide 5' -from the cleavage site, at 
position -2U. This 2' -hydroxyl group does not appear to 
participate in specific ground state contacts. However, 
5 deoxynucleotide substitution at this position leads to a 5- 
fold reduction in k chem . This is interesting because it 
implies that 2' hydroxyl groups at nucleotides adjacent to 
the cleavage site can be directly involved in the transition 
state. Perhaps the 2' hydroxyl group at position -2U is 
10 involved in a tertiary contact that generates tension, or 
strain on the adjacent scissile phosphodiester linkage. 
This 2' -hydroxyl group may be close enough to act as an 
outer-sphere ligand for a catalytic metal ion. 

15 Rapid deoxynucleotide cleavage by a group II intron 
ribozyme : 

Dl cleavage of an all-ribose substrate is slower than that 
of many known ribozymes. The degree to which the Dl 
ribozyme is considered "slow" must be put in relative terms 

2 0 since it still provides a 10 7 -fold rate -enhancement over the 
uncatalyzed cleavage reaction [Michels and Pyle, 1995] . 
However, compared to ribozymes derived from group I introns 
and RNase P, its chemical rate for oligonucleotide cleavage 
is inefficient [Herschlag and Cech, 1990; and Beebe and 

25 Fierke, 1994] . One explanation for this inefficiency is 
that the D1/D5 ribozyme does not contain all of the active - 
site functionalities necessary for efficient cleavage, 
including contacts to the 2' -hydroxyl group at the cleavage 
site. This effect might quantitatively account for the -5 

30 kcal/mol by which this ribozyme is chemically slower than 
others [Herschlag et al . , 1993A; Smith and Pace, 1993; 
Herschlag and Cech, 1990; and Beebe and Fierke, 1994] . It 
was therefore important to determine if the 2' hydroxyl 
group at the cleavage site contributes to catalysis by a 

35 fast group II intron ribozyme. 



The D1/D5 ribozyme can be altered such that it catalyzes 
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highly efficient oligonucleotide cleavage. In particular, 
when Domain 3 (-80 nucleotides in aiSg) is covalently 
attached to Domain 1. one observes an apparent cleavage rate 
approaching 1 min" 1 for rS. For the fast variants, such as 
5 D123 described here, deoxynucleotide cleavage is also 
efficient, proceeding with a rate of 0.17 min" 1 (Fig. 20). 
This is 90 -fold faster than deoxynucleotide cleavage by the 
D1/D5 construct. The D123 ribozyme can also cleave an all- 
OA oligonucleotide, apparently overcoming difficulties 
10 associated with a totally hybrid IBS/EBS duplex. The ease 
of DNA cleavage and the similar rates for rS and -IdC 
cleavage by D123 suggest that this ribozyme, like its slower 
relative Dl, does not utilize the 2'-hydroxyl at the 
cleavage site for extensive transition-state stabilization. 



15 



The k c . t value of 0.17 min 1 for -IdC cleavage by D123 is 
orders"of magnitude faster than deoxyribonucleotide cleavage 
by any known natural ribozyme [Lesser et al., 1990; and 
Walsh, 1979] . It is similar to the k c . t for cleavage of DNA 
20 by EcoRI (0.72 min" 1 ) and only an order of magnitude slower 
than the individual rate constant for the bond-breaking step 
[Kunkel et al., 1991; and Scaringe et al., 1990]. The k clt 
for -IdC cleavage by D123 is only two-fold slower than the 
rate observed for the best mutant group I intron ribozyme 
25 that was evolved through 27 generations of mutation and 
selection to cleave a DNA linkage [Tsang and Joyce, 1994] . 
However, the value for the D123 ribozyme is ' 3.2 x 

10 7 M^min" 1 (estimating a K, of 5 nM for S binding to D123) , 
a value which is 10- fold faster than W*. °f the best DNA- 
30 cleaving ribozyme derived through directed molecular 
evolution [Tsang and Joyce, 1994] . This result is primarily 
significant because it was obtained with components of a 
group II intron active site that have not been selectively 
optimized in-vitro for increased catalytic activity or for 
35 activity specifically directed against a DNA linkage. 
Although detailed kinetic characterization of D123 and D13 
ribozymes is currently under investigation, the results 
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presented here suggest that oligomicleotide-cleaving 
ribozymes of the Dl or D123 class, together with catalytic 
D5, exploit a catalytic strategy in which the 2' -hydroxy 1 
group at the cleavage site is not central to stimulation of 
5 catalysis. 

Efficient deoxyribose cleavage by a group II intron ribozyme 
may have implications for the apparent mobility, or 
migration of group II introns into new genomic positions 

10 [Lambowitz, 1989; and Lambowitz and Belfort, 1993] . It is 
widely speculated that free group II introns can bind and 
then reverse -splice into new RNA molecules, provided that 
these share sequence similarity to the parent site [Van Dyck 
et al., 1995; and Schmidt et al . , 1994]. It may also be 

15 important to consider that group II introns could attack and 
integrate into DNA sites, thus becoming directly 
incorporated into a new genome [M6rl et al., 1992] . Despite 
the temptation to speculate on a biological role for DNA 
cleavage by these ribozymes, it is important to note that 

20 the catalytic behavior of group II intron ribozymes that 
cleave oligonucleotide substrates is often different than 
that observed for cis-splicing and reactions in which the 
cleavage site is covalently attached to Dl [Podar et al., 
1995] . Although they are more difficult to analyze 

25 kinetically, cis-cleaving constructs will need to be 
substituted with deoxynucleotides before strong statements 
can be made about the function of 2' -hydroxy 1 groups in 
specific cis-splicing or reverse-splicing reactions. 

30 Ground-state effects of 2'-hydroxyls on the substrate: 

For the D1/D5 ribozyme, the effect of deoxynucleotide 
substitutions on ground state binding was measured by 
changes in the K„ value. This is considered a reasonable 
measure of Dl-S association energy because it was previously 

35 shown to equal K<, for binding of rS [Michels and Pyle, 1995] . 
Because k cat = k chftra for each substrate investigated here, it 
is likely that K^K^ for the chimeric substrates as well. 1^ 
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values for the -ldC and -2dU substrates are only 3 fold and 
2- fold larger than for rS (Table 4) . These effects are 
similar to the magnitude of duplex weakening effect observed 
previously upon single deoxynucleotide incorporation into an 
5 RNA/RNA helix [Bevilacqua and Turner, 1991; and Hall and 
McLaughlin, 1991], suggesting that these 2' hydroxyl groups 
do not form tertiary contacts in the ground state. It is 
also noteworthy that in the group I Tetrahyroena and RNase P 
ribozymes, small ground-state effects of -3 fold were also 
10 observed for substitution at the -1 position [Herschlag et 
al., 1993A; Smith and Pace, 1993; and Herschlag et al., 
1993B] . 

Single deoxynucleotide substitutions within the IBS1 region 

15 at positions -3 to -7 resulted in no observable effect on 
binding or chemistry (Fig. 16, Fig. 21). Being far away 
from the splice site, it is not surprising that these 2' 
hydroxyl groups contribute little, if any, energetic 
stabilization in the transition state. The lack of any 

20 major effects on the rate at 100 nM Dl and 3 |iM D5, which 
represent the lower limit for saturation of the rS 
substrate, also suggests that the ground state has not been 
substantially disturbed and that there are no ground-state 
tertiary contacts between any of these 2' hydroxyl groups 

25 and functionalities on the ribozyme. A lack of backbone 
tertiary interactions in oligonucleotide binding by group II 
intron ribozymes is supported by the fact that ribozyme- 
substrate binding energies parallel the amount of duplex 
stabilization one can calculate for base-pairing of the two 

30 EBS-IBS interactions [Michels and Pyle, 1995; and M6rl et 
al., 1992]. This contrasts with the behavior of group I 
intron ribozymes, in which 2' hydroxyl groups were found to 
make specific tertiary contacts that provide extra binding 
energy that stabilizes both the ground and transition states 

35 [Pyle et al., 1992; Herschlag et al., 1993; Strobel and 
Cech, 1993; Pyle and Cech, 1991; and Bevilacqua and Turner, 
1991] . 
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Iraplications for target sequence- specificity: 
Because 2' hydroxyl groups appear to make minimal 
contributions to stabilizing S-Dl association, group II 
5 introns may be able to obtain a higher level of specificity 
in substrate recognition. Ribozymes that rely on extensive 
tertiary interactions with the backbone of their substrates 
can suffer from a loss of sequence specificity. This is 
because, regardless of base identity, every residue 

10 possesses 2 f hydroxyl and phosphate groups. If these 
functionalities are important in substrate binding, they may 
undermine the energetic differential between matched and 
mismatched base -pairings, allowing a mismatched substrate to 
be "tolerated" long enough in the binding site to undergo 

15 cleavage [Herschlag, 1991] . By excluding 2' hydroxyl 
contacts from the total ribo2yme- substrate stabilization 
energy group II introns may experience an enhanced "signal 
to noise", in which nonspecific background energy is lowered 
and the contrast between canonical and non-canonical base 

20 pairing stands out. 

An exceptionally high level of specificity in substrate 
recognition by group II intron ribozymes may also result 
from the presence of two separate substrate binding domains: 

25 the EBSl and EBS2 regions, which bind to noncontiguous 
portions of Domain 1. By breaking the substrate pairing 
sequences into pieces and putting linkers between them, it 
may be possible to target a highly distinctive sequence <in 
this case, 16 nucleotides including the linker region) 

30 without adding so much binding energy that mismatches within 
the stretch will be tolerated. It is worth noting that 
targeting a unique site within the human genome requires 
base-pairing recognition of a site that is at least 16 
nucleotides in length [Sambrook et al., 1989]. The actual 

35 sequence of the EBS-IBS interactions is apparently 
unimportant, since the EBS-IBS interactions covary to retain 
base-pairing in group II introns of different sequence 
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is not the result of lariat degradation. Instead, it is a 
linear species resulting from hydrolytic, rather than 
branch-point attack at the 5' -splice site. Given that the 
5'-exon product will have the same 3'-hydroxyl terminus, 
5 regardless of whether it was produced through branch-point 
or hydrolytic attack, it should be possible to follow a 
hydrolytic first step with a second step of splicing. In 
the end, this reaction would lead to spliced exons and the 
production of linear intron, much like that observed for 
10 group I intron splicing (Cech, 1993) . 

A hydrolytic pathway, occuring during the first step, was 
previously proposed to account for the abundant linear 
intron produced at high concentrations of KCl (Peebles et 

IS al., 1987; Jarrell et al., 1988b). We reasoned that, if the 
pathway could occur in the presence of KCl (which also 
produces significant amounts of lariat intron) , there might* 
be some propensity for this pathway to occur under all 
conditions. Specific ions would be expected to increase the 

20 relative rate of one of the two possible splicing pathways, 
resulting in apparent dominance of a particular reaction 
product. To visualize this effect qualitatively, we 

compared the relative proportions of linear and lariat 
intron, produced under KCl and NH 4 C1 conditions (Fig. 23B and 

25 23C) . Like (NH«) 2 S0 4 , splicing in the presence of CI 
results predominantly in lariat evolution (Band A, Fig. 
23C) . However, large amounts of linear intron are produced 
as well (Band C, Fig. 23C) . In KCl, which was previously 
reported to promote formation of linear intron (Peebles et 

30 al., 1987; Jarrell et al . , 1988b), we observe the production 
of both linear and lariat species (Fig. 23B) . The amount by 
which lariat dominates over linear products in NH«C1 is 
reversed under KCl conditions, and neither product dominates 

by more than several- fold. This reasoning suggested that 
35 the same set of reactions were going on under all of the 
different reaction conditions. In each case, two pathways 
for splicing were evident: branch-point attack on the 5'- 
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splice site competes with hydrolytic attack on the 5' -splice 
site and, in each case, the free 5'-exon can go on to 
participate in the second step of splicing. 
Kinetic analysis of two competing splicing pathways: 
5 if linear and lariat molecules arise from competing 
mechanisms of cleavage at the 5' -splice site, then self- 
splicing will follow -parallel -first order reaction 
kinetics" as defined previously (Fersht, 1985; Jencks, 
1987). A first-order treatment is appropriate, in spite of 

10 the two-step nature of self -splicing, because intermediate 
products (such as lariat 3'-exon) do not accumulate in 
abundance and the first step of self -splicing is rate- 
limiting (Jacquier and Jacquesson-Breuleux, 1991) . 
Parallel kinetics of self -splicing are illustrated in in 

15 Scheme 1: 

K hy linear intron 

Group II * 
Intron Precursor 
20 lariat intron Scheme 1 

Using this form of kinetic analysis, we sought to define 
reaction efficiency as the actual rate for individual 
pathways of splicing, i.e. branching and hydrolysis (k^ and 

25 k hy ) under a variety of reaction conditions. In order to 
solve for k br and k hy , one must note three features of 
parallel first -order reactions (Jencks, 1987) (1) All 
products will evolve at the same total rate (k^ , which is 
equal to the sum of the individual rates (k hy and k te ) . (2) 

30 The ratio of the different products (f ^ / f br , defined in 
Materials and Methods) will be a constant through time. (3) 
If parallel kinetics apply in this case, then the sum of 
and f br will equal 1 even if equations 1-3 (Materials and 
Methods) are fit independently (as performed here) . Given 

35 these characteristics of a parallel reaction, it should not 
be possible to obtain k hy or k br directly by monitoring 
evolution of a single reaction product. Instead, one must 
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solve for k hy and k ^ by determining the relative fractional 
amounts of each product and using these values to correct k 
t (equations 4 and 5). 

5 Kinetics of self -splicing in 100 mM MgCl 2 : The first data 
sets to be analyzed were those reaction conditions that 
included Mg 2 *, but lacked monovalent ion. Commonly referred 
to as "low-salt" conditions, these are known to stimulate a 
slow self-splicing reaction (Peebles et al., 1986; Peebles 

10 et al., 1987). The most common "low salt" condition 
currently in use is 100 mM MgCl 2 (Chanfreau & Jacquier, 1994; 
Boulanger et al . , 1995). It has been reported that 
mutational sensitivities of splicing under this condition 
show the strongest parallels to effects observed in-vivo 

15 (Boulanger et al., 1995). Extended timecourses under this 
condition revealed a slow reaction following simple kinetic 
behavior that best fit equations 1-3. 

Important features of the time-course include the following: 
2 0 First of all, the overall k T values determined in three 
separate ways - reaction of precursor, evolution of lariat 
or evolution of linear intron - are all in close 
correpondence, differing by less than a factor of two 
(legend, Fig. 24A) . This correspondence implies that the 
25 evolution of lariat and linear intron follow parallel first- 
order kinetics and thereby evolve simultaneously and 
independently from precursor. The evolution of spliced 
and free exons also resulted in a rate equivalent to kj. (data 
not shown) . However, because the exons were smaller than 
30 intronic products, they were less radioactive and thus less 
accurate for tracking the course of the splicing reaction. 

In all cases, the evolution of intronic products has been 
correllated with the same extent of exon splicing. 

35 As in other reactions that did not contain molar quantities 
of ammonium salts, evolution of linear intron dominated over 
lariat intron (Fig. 24A) . The latter observation, together 
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with poor fit to the sequential reaction model (Tinoco et 
al., 1978), makes it unlikely that a sequential pathway for 
linear formation from lariat is the operative mechanism. As 
will be seen by comparing these data with reaction 

5 conditions involving high amounts of monovalent cation 
(below), the entire population of molecules in Mg 1 ' salts 
alone reacts slowly, with rates similar to a subset of 
molecules that react slowly in the presence of monovalent 
ions (Table 5) . For this reason, all of the rate data for 
10 the homogeneous population of molecules reacting in 100 mM 
M g'- or 10 mM Mg" , is located in section B of Table 1. The 
overall rate of splicing observed under these conditions was 
0.0049 min- 1 . From this value, and from the asymptotic 
coefficients for linear and lariat accumulation (f hy and f te , 

15 Table 5), we calculated k*- 0.0040 min'* and k* - 0.0009 min" 
1 (Table 5, legend to Fig.24A) . 
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Table 5 Rate Constants* for Self Splicing 
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rates, and not the rate of a specific chemical or 
conformational step. Although there is evidence to 
suggest that the hydrolytic rate constants shown here 
represent a chemical rate and that the burst rate for 
branching is a lower limit for the chemical step in 
that reaction, it was not within the scope of this 
study to dissect specific rate constants along the 
splicing pathway. 



WO 96/22689 



PCTAJS96/01337 



-lie- 
fa. Fast and slow populations of molecules refer to the two 
conformational states identified during kinetic and 
chromatographic analysis of the precursor RNA. All 
rate constants for reaction in the absence of 
5 monovalent ion (the first two rows) are put in the 

category of Population B because their rates matched 
those of Pop. B under all other reactions conditions. 
However, given the slow rates and limited reaction 
extent in 10 mM MgCl 2 and 100 mM MgCl a , one cannot rule 
10 out the existence of a second, less active, 

conformation under these conditions. 

c. In addition to the salts listed, each reaction 
contained 40 nM MOPS, pH 7.5. 

d. The variables k x , fc r , f , ^ anc^ f for a given 
15 population (A or B) are the total rate of splicing, the 

rate of branching, the fraction of molecules that 
splice through branching, the rate of hydrolytic 
splicing, and the fraction of molecules that splice 
through hydrolysis, respectively. 
20 e. The rate of reaction for the 13% of the total 
population that undergoes a burst of reaction. 

f . The remainder of population A that does not burst, but 
still reacts through branching at a fast rate. 

g. The top value (13%, determined independently, Figure 5) 
25 represents the fraction of molecules that proceed 

through a burst. The bottom value represents the 
remainder of the fast population. 

Kinetics of self-splicing in 10 mM MgCl 2 : The first set of 
3 0 ionic conditions ever used to detect group II intron self- 
spl icing were 10 mM Mg 2 * with 2 mM spermidine (Peebles et 
al., 1986). Under these conditions, lariat and linear 
products form in approximately eguimolar ratios. However, 
in our hands, RNA degradation was so severe that it was 
35 difficult to perceive the accumulation of either spliced or 
free exons (lanes 15-17, Fig. 28) . The data was of 
sufficient quality to quantitate an overall rate of 
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precursor reaction (0.003 min -1 , data not shown), but 
degradation made it difficult to accurately partition this 
rate between relative contributions from linear and lariat 
pathways. Because of this problem, we found an alternative 
5 polyamine for analyzing activity at 10 mM Mg 2 \ We observed 
that stoichiometric amounts of protamine, a small arginine- 
rich polypeptide, promoted significant amounts of reaction 
at 10 mM Mg 2 *, without concomitant RNA degradation. This 
could be attributable to the formation of a tightly folded 

10 RNA structure in the presence of protamine, or an inability 
to remove traces of nuclease from the spermidine. In the 
presence of 5 nM protamine, the accumulation of lariat RNA 
is similar to that observed using 10 mM Mg 2 * with 2 mM 
spermidine, although all the linear products of self- 

15 splicing (such as linear intron, spliced exons and free 
exons), which tend to be more susceptible to degradation, 
can be quantitated because they are readily detected above 
the noise of random RNA degradation. As in the case of 100 
mM Mg 2 \ the reaction obeys parallel kinetics in which the 

20 linear intron product dominates over lariat intron product 
(Fig. 24B) . That both pathways lead almost exclusively to 
the production of ligated exons (rather than free exons) is 
apparent from product analysis (lanes 18-20, Fig. 27) . The 
overall rate of splicing observed under these conditions was 

25 0.0029 min' 1 (as seen for 10 mM Mg2+ and 2 mM spermidine). 
From this value, and from the asymptotic coefficients for 
linear and lariat accumulation (f^ and fb r ), we calculated 
k hy = 0.0025 min" 1 and « 0.00044 miri . These rates are 
both approximately two-fold slower than in the presence of 

30 100 mM Mg 2 *. 

Kinetic analysis of self -splicing with high concentrations 
of monovalent ion: 

Previous studies reported that monovalent cations had an 
35 inhibitory effect on the rate of self -splicing in the 
presence of 10 mM Mg 2 * (Peebles et al., 1986), although high 
concentrations of certain salts had a marked stimulatory 
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effect on rate in the presence of 100 raM Mg 2 * (Peebles et 
al., 1987). Presumably this is because monovalent ions 
facilitate folding at high concentration. However, at high 
concentration, monovalent ions may be able to compete with 
5 low concentrations of catalytically essential divalent ions 
in the active site. Thus, one might be required to 
correspondingly increase the concentration of divalent ion 
in order to exploit the stimulatory effects of monovalent 
ions. Several studies defined three salts that are 

10 particularly stimulatory for splicing in the presence of 100 
mM MgCl 3 at concentrations from 0.2M to 2M: KC1, NH«Cl and 
(NH,) 2 S0 4 (Peebles et al . , 1987; Schmelzer & Muller, 1987; 
van der Veen et al., 1987; Jarrell et al., 1988b). While 
these are the major salts employed in group II self -splicing 

15 studies today, we also examined several others, including 
potassium acetate . 

Biphasic kinetic behavior: Although they appear to favor 
different pathways during the first step of splicing (Fig. 

20 24), plots of the overall self -splicing reaction in 
<NH 4 ) 2 S0 4 , KC1 and NH^Cl are remarkably similar. Each of them 
shows pronounced biphasic behavior (Fig. 25) with a 
transition at approximately 50% reaction extent. That is, 
the data for each timecourse fall on two straight lines, 

25 which intersect at - 50% reaction completion. Because this 
apparent switch from fast to slow reaction occurs at 
approximately the same reaction extent under such radically 
different reaction conditions, it suggests that there are 
two major conformations of the precursor RNA, resulting in 

30 two populations of molecules: a fast population and a slow 
population. The slope of the line during the latter 50% of 
reaction is almost identical under all three reaction 
conditions shown. Thus, the slow population reacts 
similarly, regardless of ionic environment. Given the fact 

35 that the rate transition occurs at the same extent of 
reaction but at different times under the different ionic 
conditions (-10 minutes for NH 4 S0 4 ; -30 minutes for N^Cl, 
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-60 minutes for KC1) , the two conformations of molecules 
appear to be in slow exchange, if they exchange at all. 
This finding indicates that a kinetic treatment of group II 
intron self -splicing in monovalent ions must accomplish two 
5 objectives: It must account for the two chemically distinct 
pathways for splicing (branching and hydrolysis) . At the 
same time, it must also mathematically account for two non- 
interconverting conformations in the precursor RNA. 

10 Kinetic treatment of two RNA conformations: The preceding 
data suggest that biphasic reaction behavior is due to the 
presence of two populations of molecules, irrespective of 
the type of reaction pathway that dominates during splicing. 
In addition, the data suggest that biphasic behavior should 

15 not be attributed to a sequential pathway in which lariat 
serves as an intermediate for production of linear intron. 

This is because self -splicing proceeds almost exclusively 
through a hydrolytic pathway in the presence of KC1, 
although the same biphasic kinetics are observed under that 

20 reaction condition. The most simple kinetic model for the 
results is presented in Scheme 2, in which the precursor 
consists of two non- interconvertible populations (A and B) - 
each of which is catalytically active. Each of these 
populations may be able to splice through the hydrolysis or 

25 branching pathways. Each population may therefore react 
with the same parallel- first -order kinetic behavior outlined 
in previous sections. 

A -* linear intron 
30 Group II + I h k\ r 
Intron X 
Precursor J **k B hy 

B -» lariat intron 

k 9 ^ Scheme 2 

35 The presence of two reactive non-interconverting populations 
requires that the raw data in Fig. 26A-C fit two sets of 
exponential terms instead of one, as shown previously. 
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For this purpose, equations 1-5 describing parallel first- 
order behavior were modified such that product evolution is 
described by an equation containing two sets of terms: one 
describing population A and one describing population B 
5 (equations 6-10). As before, if parallel kinetics is 
followed, k T values for a single population (k, for example) 
will be the same, regardless of whether linear, lariat or 
precursor data is fit. Fractional values for a given 
population will remain a constant throughout reaction, and 
10 the four fractional coefficients describing extent of 
formation of intronic products will add up to one. 

Self -splicing in 100 mM Mg 2 * and 500 mM KC1: The data 
obtained under this condition was best fit using 

15 biexponential equations 6-8, as predicted from the first 
order plots (Fig. 25C) . Both the linear and lariat curves 
fit biexponential equations 7 and 8, rather than a 
sequential model (Tinoco et al., 1978), suggesting that the 
pathway outlined in scheme 2 is correct. Both populations 

2 0 of molecules partitioned between linear and lariat reactions 
(Fig. 26A) . As expected for a parallel first order 
reactions from each population i, all kvalues matched and 
all four fractional values describing each population added 
to one (legend. Fig. 26A) . It is important to note that 

25 even the fast population of molecules (totaling 41% of the 
population in this experiment) reacted both through the 
linear and lariat pathways in competition with one another. 
The apparent rate of the reaction by population A was 0.028 
min' 1 (note that error limits for all rates are provided in 

30 the caption to Fig. 26) . By solving equations 9 and 10 for 
population A, we observe that hydrolysis is the fastest 
process, with an individual rate k= 0.022 min" 1 , compared to 
k= 0.0063 min" 1 for branching (Table 5, legend, Fig. 26A) . 
For the slow conformation (population B) , hydrolysis is also 

35 faster than branching. The overall rate kj. * 0.0042 min' 1 
and, by solving equations 9 and 10, we observe that the k= 
0.0032 min' 1 and k= 0.0011 min 1 . These latter rate constants 
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for B are remarkably similar those observed for splicing of 
the entire precursor population in 100 mM MgCl 2 alone (Table 
5) . 

5 Kinetics of self -splicing in NHiCl: Another monovalent salt 
that is a standard condition for studies of in-vitro 
splicing is NH 4 C1 (Peebles et al., 1987; Muller et al., 1988; 
Jarrell et al., 1988b; Suchy & Schmelzer, 1991; Winkler & 
Kuck, 1991; Peebles et al., 1993). Time-courses of reaction 

10 in 100 mM MgCl 2 and 500 mM NI^Cl are shown in Pig. 26B. 
There are several trends that are apparent from inspection 
of the raw data. First of all, the lariat molecules are 
formed during the first 50% of reaction, after which there 
is a plateau in lariat production. Large amounts of linear 

15 intron are observed only during the second half of reaction. 
A simplistic interpretation of these results might be that 
the linear band is due to lariat decomposition. But, after 
lariat begins to plateau (-50 minutes) , all of the evolution 
of linear intron is due to decrease in precursor (Fig. 27B) . 

2 0 Once lariat evolution is a flat line, the curves for 

precursor and linear intron change inversely with respect to 
one another. It is clear that most of the lariat that is 
produced over time is stable because lariat molecules remain 
the dominant product after 15 hours of reaction (Fig. 27, 
25 lane 5) . Therefore, even by this qualitative analysis, it 
appears that linear intron is produced independent of 
lariat, through a mechanism in which a hydrolytic pathway 
competes with branching. 

30 Decrease in precursor and evolution of lariat are best fit 
by biexponential equations 6 and 8, suggesting that lariat 
is produced from two conformations of precursor molecules. 
However, the evolution of linear intron fits only a single 
exponential (equation 2) . The lack of a biexponential fit 

3 5 for the linear intron suggests that linear intron is 

produced only from the slow population (B) . That B produces 
both linear and lariat intron through parallel first -order 
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kinetics is supported by the fact that k T obtained from 
equation 2 is the same as k, obtained from the fit of the 
lariat data to equation 8. Taken together, these data show 
that fast population (A) of precursor reacts solely through 
5 a lariat pathway, with a rate of 0.038 min' 1 . In ammonium 
salts, the hydrolytic pathway from conformation A is 
completely shut down. Under these conditions, this 
conformer may be able to exclude water and promote an 
exclusively branch-mediated splicing pathway. Conversely, 

10 the slow population (B) obeys parallel first order kinetics, 
partitioning through competing linear and lariat pathways. 
Under these reaction conditions, Scheme 2 must be slightly 
modified to a form in which there is no arrow connecting 
population A with the linear intron (the kpathway) . Taking 

15 the precursor depletion as the most accurate representation 
of overall reaction for k, and using the relative amounts of 
linear and lariat produced by that (f- 0.51, f- 0.10), as 
coefficients for correcting the data, the rates of splicing 
from B were calculated to be k- 0.0032 min- and k= 0.00063 

20 min 1 . These rates of reaction by population B are 
remarkably similar to those observed under the other 
reaction conditions (Table 5) . 

Kinetics of self -splicing in (NH,),S0 4 : Li*e NH.C1. lariat 
formation is radically stimulated by (NH,) a SO« (Fig. 26C) . 

25 It appears that population A splices only through a 
branching pathway, resulting only in lariat intron. One 
salient difference in the kinetic behavior in (NH,),S0 4 salt 
is the existence of a considerable burst early in the 
reaction. When plotted in a semi -log format, in addition 

30 to the sharp demarcation observed at 50% reaction (between 
fast and slow populations of molecules) , there is a large 
burst consuming 13% of population A during the first 60 
seconds of reaction (inset, Fig. 25A) . We observe that, as 
in the case of NH«Cl, the linear species evolves very slowly 

35 and fits a single exponential. After a plateau in lariat 
formation, linear evolves only at the expense of precursor 
(Fig. 26C) . 
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The timecouraes in (NH 4 ) 2 S0« were fit using methods similar 
to that described for NH 4 C1, with one exception designed to 
account for the burst of reaction that immediately consumes 
5 13% of the precursor (Fig. 25A, inset). This fractional 
amount was used to add an additional rate term (for which 
the f = 13%) to equations 6 and 8 so that the remainder of 
the population could be fit precisely. The equations were 
extended in the same way as equations 6 and 8 were adapted 
10 from equations 1 and 3, although the new f value was held 
constant. This served to reduce the uncertainty of adding 
an additional set of variables to the fit. After fitting 
the precursor and lariat data to the modified equations 6-8 , 
it was possible to solve for individual rates constants for 
15 branching by population A. The burst population 

(contributing a fraction of 0.13) reacted with a rate of 
0.22 min" 1 . The remainder of population A (f) reacted solely 
through a branching pathway, with a rate of 0.047 min" 1 
(Table 5, Fig. 26C) . The slow population B (k T = 0.0028 
20 min" 1 ) reacted through both hydrolytic <f= 0.39) and 
branching (f= 0.22) pathways at respective rates of k« 
0.0018 and k» 0.0010 min* 1 (Table 5). Population B, 
therefore shows competitive behavior and rate constants 
similar to the other reaction conditions. Population A, 
25 however, is quite different in that it undergoes a rapid 
burst of lariat formation tnat is among the fastest rates 
seen for reactions by group II intron. This confirms 
earlier observations indicating that, if it weren't for 
conformational states that obscure it, the chemical rate of 
30 branching would appear to be extremely rapid (Chin and Pyle, 
1995). The value of the burst provided here (0*22 min" x ) , 
may be viewed as a lower limit for the chemical rate of 
branching. 

35 The time-courses of splicing under different reaction 
conditions display mechanistic features that argue against 
a sequential pathway: (a) There is no "hump" in lariat 
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formation, which would appear as an initial burst followed 
by a steady decrease in the lariat population (Fersht, 
1985) . Such a rise and fall has been observed before under 
conditions where the branch is destabilized (Chin and Pyle, 
5 1995) . (b) The data fit equations for parallel rather than 
sequential kinetic models. (c) Even when linear intron 
dominates, lariat molecules continue to evolve throughout 
reaction, arguing that their conversion to linear molecules 
is not so rapid that lariats cannot accumulate. (d) Linear 
10 appears to evolve quantitatively at the expense of 
precursor, rather than at the expense of lariat. 

Ribonuclease analysis of the linear intron species: 
Large RNA molecules, whether linear or lariat are inevitably 
15 susceptible to some degree of degradation at long times, 
one cut in a lariat molecule will linearize it, creating a 
species that migrates with hydrolytically generated intron. 
The branched "tail" protruding from the linearized lariats 
is very small, so it is expected to have only minor effects 
20 on mobility (Peebles et al . , 1986). Although the data 
presented here are consistent with a hydrolytic pathway for 
production of linear intron, a fraction of the linear 
species is expected to result from lariat decomposition. 
This is supported by the kinetic data, which show that all 
25 k T values determined from lariat evolution are slightly 
larger than the k* values determined from linear intron 
evolution. This effect stems from an artificially early 
plateau in the evolution of lariat, giving rise to an 
artificially late plateau in linear intron production (due 
30 to slow conversion of lariat to linear) . This results in 
a shorter t l/2 for lariat and a longer t ul for linear intron 
evolution. Because we did not expect the linear intron to 
be purely the result of a hydrolytic pathway, we were 
interested in determining the relative stoichoimetry of 
35 branched to mononucleotides in the linear band. We reasoned 
that linear intron isolated after long times of reaction 
might contain more broken lariat than linear species 
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isolated after short times of reaction. 

In order to obtain direct physical information on the 
branched trinucleotide content of the linear intron, we 
5 subjected the linear intron species to chemical 
characterization, through digestion with PI nuclease 
Lariat intron contains a single branched trinucleotide, of 
composition pA (2' -pG) 3 9 -pU. PI nuclease is unable to cleave 
the phosphodiester linkages within a branched trinucleotide 

10 (Reilly et al., 1989). Therefore, complete digestion of the 
lariat RNA with PI nuclease results in abundant 
mononucleotides and a single branched trinucleotide. These 
two species are readily separable by PAGE and can be 
quant itated by radioanalytic scanning if one of the 

15 components of the branch is labelled with 32 P (Moore & Sharp, 
1992; Chin and Pyle, 1995). When linear intron results from 
lariat breakage, it will contain a branched, 2' -5' linkage 
and PI digestion will result in a particular stoichiometry 
of mono and trinucleotides. 

20 

Of the three nucleotides contained in the branched species 
(pA(2'-pG)3'-pU> , guanosine is the least abundant in the 
rest of the intron. Therefore, we transcribed precursor RNA 
with af- 32 P-GTP and subjected this RNA to three diverse 

25 conditions that promote efficient splicing: 100 mM MgCl 2 with 
500 mM KC1, NH 4 C1 or (NHJjSQ,. Large scale aliquots were 
removed at 50% and 90% extents of total reaction and 
products were separated by preparative -scale PAGE. Lariat 
and linear intron bands were isolated from the gel, and the 

3 0 RNA was digested to completion with PI nuclease. After 
product separation, the nucleotide ratio in the lariat band 
was approximately 120:1 under all conditions (Table 6). 
There are 119 guanosines in the intron and therefore PI 
digestion of the lariat RNA should result in an ideal mono 

35 to trinucleotide ratio of 118:1, which is very close to the 
observed value. This confirmed the precision of analysis 
and suggested that 2' -5' branched linkage was stable during 
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isolation from the gel. 

Table 6: PI ribomiclease analysis of intronic phosphodiester 



Condition 


Lariat Band 


Linear Band 


Fraction of 

linear lnLruu 


(buffer 

sa J.C « 

minutes 


mono/branch 
nucleotides 


mono/branch 

nUClcOLlUcb 


conuaining a 
branch c ( for whole 
molecules) 










(NH 4 ) 2 S0 4 








30 


121 b 


615 b 


1 Q "7-1. 1 Q !fc 


180 


118 


254 


46.4±3.5% 










NH«Cl 








30 


112 


667 


16.8±1.3% 


180 


114 


237 


48.1±3.6% 










KC1 








40 


105 


523 


20.1±1.5% 


180 


119 


297 


40.1±3.0% 



10 



15 



20 



Table 6 Legend: Splicing reactions were conducted under the 
25 conditions indicated, followed by electrophoresis and 
isolation of individual intronic products (RNA was 
internally labelled by incorporation with 32 P-GTP) . The 
separate lariat and linear products were digested to 
completion with PI nuclease and relative proportions of 
30 mononucleotide to branched trinucleotide in each species 
were compared. There are 118 guanosine residues in the 
intron, one of which becomes incorporated in the branched 
trinucleotide upon adenosine attack at the 5' -splice site. 
A ratio of mono/branched nucleotide that is significantly 
35 greater than 118 indicates formation of linear intron 
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through an event that does not involve branching (i.e. 
hydrolysis during the first step of splicing) . a Salt 
indicated was present at 0.5 M, together with 0.1 M MgCl 2 . 
b Variance in each nucleotide ratio is ±5.8% with 95% 
5 confidence, representing error in lariat values show 
representing six different trials, ^rror is calculated with 
95% confidence for four separate trials, represent error 
between values obtained in separate experiments performed on 
different days. 

10 

Pi analysis of the linear intron band showed that, during 
the first half of reaction, only 20% of the linear species 
contained a branched trinucleotide (Table 6) . Random RNA 
degradation should increase with time and therefore it was 

15 notable that at 180 minutes, after the latter half of 
reaction in each case, the amount of linear species that 
arose from lariat molecules was approximately 45% (Table 6} . 

The late timepoints are important because they are an 
internal control, demonstrating that we would observe 

20 branched trinucleotide if it were there. These timepoints 
also show that, when significant amount of RNA degradation 
is allowed to occur, one will observe that the linear band 
contains a significant fraction of branched trinucleotide. 
Nonetheless, in our system, the actual amount of branched 

25 species in the linear product was never dominant, even at 
180 minutes. Thus, the nucleotide ratio obtained after Pi 
digestion of the linear band indicated that the major 
component was true linear intron. 

30 Spliced-exon reopening is unrelated to hydrolytic splicing: 
During in-vitro self -splicing of the ai5g intron, one 
commonly observes free 5' and 3' exons (Fig. 23, 27). It 
has been proposed that the free exons arise from a secondary 
reaction involving rapid hydrolysis of spliced exons 

35 (spliced-exon reopening, or SER) (Peebles et al. r 1987; 
Jarrell et al., 1988b). This latter model is credible given 
that the ratio of 5' to 3 '-exons is always one (data not 
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shovm) and it was unlikely that hydrolysis at both splice 
sites would take place at identical rates. In addition, 
timecourses of SER reveal that spliced exons predominate 
early in reaction, giving way over time to free exons (Fig. 
5 27) . By comparing the amount of spliced to linearized exons 
at various points during the reaction, it is possible to 
discern the relative efficiency of the SER reaction under 
different reaction conditions (Table 7) . 

10 Table 7: Dependence of the SER reaction on ionic conditions. 



15 



20 



25 



30 



Ionic Content* 
of 

Reaction 
Buffer 


Fraction of 
Ligated Exons b 
at 12% 
Reaction 


Fraction of Licrated 
Exons at 25% Reaction 








0.5 M NH 4 C1 


0.21±0.04 c 


0.11*0.03 


0.1 M MgCl 2 












0.5 M (NH 4 ) 2 S0 4 


0.63±0.13 


0.50±0.11 


0.1 M MgCl 2 












0.5 M KC1 


0.071±0.01 


<0.01 d 


0.1 M MgCl 2 












0.1 M MgCl 2 


0.027*0.005 


0.033±0.007 








5 nM 

protamine 


0.56±0.11 


0.36±0.08 


0.01 M MgCl 2 







Table 7 Legend: 

a. Other than the specific salts indicated, all reactions 
contained 40 mM MOPS, pH 7.5 and were conducted as described 
in the experimental section. 
35 b. This value is the fraction of ligated exons/total 
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released exons. Total released exons - the ligated exons 
plus the free exons produced through the spliced-exon 
(ligated-exon) reopening reaction (SER) . That free exons 
were formerly ligated and cleaved through SER is indicated 
5 by previous work, by the higher fraction of spliced exons 
early in reaction (see above), and by the fact that the 
ration of fee 5'-exon is always equal to 1. 

c. Variances shown are the standard errors calculated using 
a 95% confidence level from the mean squared deviation in 

10 six to eight trials under each condition. 

d. Spliced exons were present at levels that approach the 
limit of detection. 

Because conditions that promote hydrolytic splicing are 

15 sometimes accompanied by SER {Jarrell et al., 1988b), one 
might assume that a hydrolytic first step is accompanied by 
an unproductive side reaction that reopens spliced exons. 
If this were true, then a hydrolytic first step might be 
considered detrimental to productive splicing and branching 

20 would be considered the "proper" pathway for self -splicing . 
However, we observe SER under all reaction conditions, 
regardless of whether branching or hydrolysis dominates 
during the first step. Even reactions containing high 
concentrations of ammonium are accompanied by SER. 

25 Particularly NH 4 Cl, which promotes splicing primarily through 
lariat formation, also promotes rapid SER (Fig. 27, Table 
7) . The amount of SER determined in this case indicates 
that spliced exons are hydrolyzed almost as fast as they 
form. The same is true under conditions of 100 mM MgCl 2 , 

30 without any monovalent ion. Conversely, protamine 
stimulates a predominantly hydrolytic splicing pathway that 
leads primarily to spliced exons (Fig. 27, Table 7). An 
obligate relationship between SER and the hydrolytic first 
step of splicing is not evident when a variety of conditions 

35 are considered together. The prevalence of SER is instead 
correlated with the presence of chloride ion, which is often 
provided in molar amounts in the study of self -splicing. 
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However chloride ion is not singularly responsible for SER 
because, although it is less prevalent, the reaction is 
still observed in (NHJaSO* (Pig, 23A) and in reactions where 
chloride has been completely eliminated 
5 (K (acetate) /Mg (acetate) 2 , data not shown). 

DlSCUSi? XQN: 

The existence of two possible pathways for the first step of 
self -splicing : 

10 The evidence for a competitive interplay of two different 
reactions at the 5 '-splice site is based on a variety of 
results, including the kinetic behavior observed during 
timecourses of group II intron self -splicing (Fig. 24 and 
26) . This is supplemented by physical characterization of 

15 reaction products (Table 6) . The kinetic analysis indicates 
that the decrease in precursor RNA and the evolution of 
linear and lariat products follow a kinetic paradigm known 
as the "parallel" reaction (Jencks, 1987) . As such, the 
individual products arise simultaneously from a common 

20 precursor. More qualitatively, it is observed that the 
lariat intron is too stable to serve as the precursor for 
the majority of linearized intron observed and linear intron 
is observed to evolve at the expense of precursor, rather 
than lariat RNA. Physical analysis on the products of 

2 5 reactions that favor both hydrolytic and branching pathways 
indicates that the majority of linear molecules do not 
contain a branched trinucleotide {Table 6) . Particularly 
early in reaction, before RNA degradation sets in, linear 
intron is almost entirely composed of molecules arising 

30 through the hydrolytic pathway. Therefore this species is 
primarily the product of a different splicing pathway should 
not necessarily be considered a "broken lariat". 

The hydrolysis pathway may be a legitimate mechanism for the 
35 first step of group II intron self -splicing. After all, 
RNase P catalyzes phosphodiester cleavage using water or 
hydroxide as the nucleophile (Guerrier-Takada et al., 1983; 
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Pace & Smith, 1990) . Regardless of any hypothetical 
relationship to RNase P, it is possible that components of 
the group II intron active site have evolved specifically to 
activate a water/hydroxide nucleophile. There is little 
5 reason to assume that hydrolytic attack takes place in a 
different active site than branching, particularly since 
branch -mediated and hydrolytic attack at the 5' -splice site 
both proceed with identical preference for an Sp 
phosphorothioate diastereomer in cis-exon constructs (Podar 

10 et al., 1995a). Rather than invoking an active site that is 
specifically designed to activate a water nucleophile, one 
could alternatively argue that the "normal" nucleophile is 
the 2'-hydroxyl group of the branch-point A in Domain 6. If 
this is the case, hydrolysis may occur because the active 

15 site exerts tremendous strain on, or creates significant 
chemical activation of, the scissile bond at the 5' -cleavage 
site. As a result, other nucleophiles in the vicinity 
(particularly 55M water) may cleave the scissile 
phosphodiester linkage. 

20 

Implications of a hydrolytic pathway: 

A general hydrolytic mechanism for group II intron self- 
splicing has several important implications. The most 
important one is that, if this reaction readily occurs in- 

25 vitro, it may have a function in-vivo. Unlike the facile 
reversibility of branch-point attack (lariat formation), 
hydrolytic scission of the 5' -splice site is effectively 
irreversible (Chin and Pyle, 1995) . This has important 
implications for intron mobility, which can occur through 

30 lariat reinsertion into sites resembling the original 5'- 
exon sequences (Lambowitz and Belfort, 1993) . It is 
generally believed that group II introns are infectious 
elements that came to rest in their present locations as a 
result of reactions that resembled reverse -splicing 

35 (Lambowitz, 1989) . A possible role for hydrolytic splicing 
comes from the fact that, once a group II intron has become 
integrated into a cascade of normal gene expression, the 
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cell may try to prevent the accumulation of highly reactive 
lariat introns. Lariat introns are potentially deleterious 
because they can reverse-splice {Augustin et al., 1990; 
MYller et al., 1991), reintegrating into other RNA 
5 molecules. For certain introns with a high propensity to 
migrate, there may have been selective pressure for the cell 
to develop intron mutations that specifically promote 
hydrolytic splicing. In that way, the cell exploits the 
advantages of splicing , without the risk of further genomic 
10 alteration. 

Another possible role for a hydrolytic splicing pathway is 
related to the fact that many group II introns contain open 
reading frames for the translation of proteins called 

15 "maturases" (Lambowitz and Perlman, 1990; Lambowitz and 
Belfort, 1993) . These proteins aid in splicing and mobility 
of group II introns. If translation necessarily occurs after 
splicing (and it is not clear that this is the case) , then 
there might be serious difficulties in using lariat RNA as the 

20 message. Although translation in mitochondria is not well 
understood, potential problems with lariat translation are 
suggested by recent studies showing that cytoplasmic 
translation initiation of circular RNA molecules can only 
occur if these contain particular sequences called internal 

25 ribosome entry sites, or IRES (Chen & Sarnow, 1995). Even if 
the hydrolytic pathway were very inefficient and rare in-vivo, 
it might serve an important function: Hydrolytic splicing 
might generate a small population of linear intron molecules 
that are readily translated. This would result in the 

30 accumulation of maturase proteins, which would then facilitate 
additional rounds of self -splicing through a predominantly 
lariat pathway. By this model, linear introns may serve to 
"seed" lariat -mediated splicing by generating the proteins 
required for rapid self -splicing in-vivo. 

35 

A hydrolytic splicing pathway in-vivo is consistent with 
reaction conditions that favor hydrolytic splicing in- vitro. 
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Intracellular salt concentrations are approximately 20 mM MgCl 2 
and 140 mM KC1 (Alberts et al., 1989) which, together with 
millimolar levels of polyamines and cationic peptides, are most 
analogous to the in-vitro conditions that favor the hydrolytic 
5 splicing pathway. Indeed, the in-vitro conditions that promote 
lariat formation (molar ammonium ion) are the least 
physiological. Nonetheless, group II introns were first 
identified because of an abundance of stable circular RNA 
molecules observed during electron microscopic analyses of 

10 cellular RNA (Arnberg et al., 1980; Halbreich et al., 1980). 
There is good evidence that splicing in-vivo proceeds through 
a lariat pathway, although there is no evidence suggesting that 
this is the exclusive route or that the linear pathway does not 
occur. Indeed, lariat molecules are expected to appear more 

15 abundant because they do not have free ends and are therefore 
refractory to attack by cellular nucleases. 

Electron microscopy (Arnberg et al., 1980; Halbreich et al., 
1980; Hensgens et al., 1982; Schmidt et al., 1987), combined 
20 with enzymatic analysis of cellular RNA (van der Veen et al . , 

1986) , showed that excised group II introns were lariats, 
containing an internal 2' -5' bond and a free 3' terminus. 
However, abundant linear RNA molecules of the same contour 
length were also identified in the RNA fractions used for 

25 electron microscopy (Arnberg et al . , 1980; Schmidt et al., 

1987) . m group II intron il from the mitochondria of 
filamentous fungus Podospora anaerina, two forms of the intron 
were isolated (Schmidt et al., 1987). Mapping of the intron 
with different DNA probes indicated that these two intron forms 

30 contained the same sequences. RNAse H digestion converted the 
slower migrating form into one that comigrated with the faster 
form, indicating that the RNA transcripts were topological 
isomers of the same molecule. Debranching treatment also 
converted the slower band into one which comigrated with the 

35 faster band, indicating that the RNA in the slower band is a 
lariat form of the RNA (Schmidt et al., 19B7) . It is 
interesting to note that only one band results from the RNase 
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H treatment, and it is not accompanied by a smear that would be 
expected if the linear species contained randomly broken lariat 
molecules. The true test of a hydrolytic pathway in- vivo does 
not, however, rest with physical analysis of cellular RNA 
5 because alternative explanations (such as instability of the 
2' -5' link) can found for the results. Instead, evidence for 
the existence of this pathway lies with genetic studies on 
mutant yeast strains in which the branch-point A has been 
deleted from essential group II introns like ai5g. These 
10 studies, in which the bulged A has been inactivated in a 
variety of ways, would be the most definitive way to unravel 
the function of hydrolysis in- vivo. 

Four group II introns contained within chloroplast tRNA" 1 genes 
15 from four diverse organisms (Wicotiana tobacvm, Marchantia 
polymorpha, Zea mays, and Hordeum vulgare) are known to be 
missing the branchpoint A in Domain 6 (Michel et al., 1989). 
These four introns might be good candidates for a studying a 
class of group II intron that splices entirely through a linear 
20 pathway in vivo. Not all group II introns may be alike in 
their ability to splice hydrolytically in-vivo, so it might be 
important to test branch-point deletion on a variety of 
essential yeast group II introns. It is also possible that the 
different pathways are used under different growth conditions, 
25 or at different times during organellar biogenesis. 

Mechanistic features of the self -splicing reactions: 

One of the most pronounced effects explored in this study is 

the presence of two nonexchangable conformations in the 

30 precursor RNA. Population A (fast) can splice through either 
hydrolytic or branching pathways, although it dramatically 
prefers the branching pathway in the presence of ammonium ions. 
••A" forms only in the presence of high monovalent ion, and it 
appears to be a more reactive form of the intron. Population 

35 B (slow) can also utilize either pathway, although it prefers 
hydrolysis. "B" appears to be some sort of default form of the 
intron, with a similar kinetic profile under all reaction 
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conditions including low salt. Indeed, under low salt 
conditions, this may be the only conformation of the intron. 
Any future analysis of splicing rates will have to contend 
kinetically with these multiple conformations. Our means of 
5 dealing with this was to separate the conformations 
mathematically using multiple exponentials in our rate 
equations. Preliminary gel -filtration results (data not shown) 
suggest that it may be possible to physically separate the 
conformations and study them individually. 

10 

Conformations A and B might be related to the incubation 
procedures used in this study. Unlike any other group II 
intron splicing work published to date, the reactions shown 
here were preceded by a step in which the RNA was heated to 
15 95 ic, in order to melt alternative inactive conformations 
formed during RNA storage (walstrum & Uhlenbeck, 1990; Pyle and 
Green, 1994; Uhlenbeck, 1995). In order to determine if the 
conformational distribution was related to pre-heating, we 
conducted self -splicing reactions on precursor RNA that was not 
20 subjected to 95 |C preincubation (data not shown). We observed 
behavior very similar to that reported here with one important 
difference: Population A reacted with a rate -50% slower than 
that observed under the other conditions reported here. 
However, the same relative proportions of A and B were 
25 observed. This suggests that a short 95 ;C preincubation does 
not adjust the relative proportion of A and B (or allow their 
interconversion) , but it aids in activation or folding of 
population A. The observation of two conformers is consistent 
with a recent report from another laboratory in which, under 
30 different preincubation conditions than reported here, 40% of 
the ai5g intron was found to be less reactive (Podar et al., 
1995b) . 

The rate constants reported in Table 5 are useful benchmarks 
35 for structure/function studies on mutants of group II introns. 
However, more experiments would be required to determine if 
these rate constant represent specific events in the first step 
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of splicing, such as chemistry or conformational change. It 
may be possible to draw parallels with mechanistic analyses 
performed on mult i- component ribozyme constructs of the ai5g 
group II intron. The fast rate constant for hydrolytic 
5 splicing in KCl (k) is virtually identical to the chemical rate 
of Domain 5-catalyzed hydrolytic scission of the 5'-exon from 
Domain 1-3 (Pyle and Green, 1994) under the same ionic 
conditions. In a study that analyzed trans-D56 attack on the 
5' -splice site (Chin and Pyle, 1995), it was observed that 

10 branching was limited by the rate of conformational sampling of 
a chemically active state. This rate was analogous to kunder 
the reaction conditions reported in Table 5. Similarly, a 
burst in reactivity was proposed to represent a lower limit for 
a chemical step of branching, which is necessarily very rapid 

15 (Chin and Pyle, 1995). 

Although the reaction condition involving protamine had not 
been reported in previous studies of group II intron self- 
splicing, other charge-neutralizing proteins have been found to 

20 replace the high monovalent ion requirements of certain 
ribozymes. For example, RNase P can function efficiently at 
lower salt concentrations in the presence of the C5 RNase P 
protein (Reich et al., 1988). While there may be specific 
proteins that associate with particular group II introns in- 

25 vivo (Wiesenberger et al., 1992; Van Dyck et al,, 1995), we 
were interested in examining the effects of a common 
nonspecific RNA and DNA-binding protein. The amino acid 
sequences of protamines (short arginine-rich peptides) resemble 
isolated RNA- binding domains of about 30 amino-acids in length. 

30 Structural studies of protamine binding to RNA show that the 
polypeptide forms an a-helix, which rests in accessible helical 
grooves of nucleic acid duplexes (Warrant & Kim, 1978; Hud et 
al., 1994). Our first experiments with protamine were 
conducted at concentrations of the peptide. After 

35 decreasing the [protamine] such that there were stoichiometric 
quantities of protamine bound to the RNA, we began to see 
stimulation of catalysis at 10 mM MgCl 2 . Protamine had the 
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greatest stimulatory effect at a concentration of 5 nM (the 
precursor was generally present at 1 nM) . In addition to 
stimulating the reaction in the presence of 10 mM MgCl a (Fig. 
24B,C,D), gels of the reaction products showed that the RNA was 
5 not even slightly degraded during these reactions (Fig. 27, 
lanes 18-20) . Despite the fact that most of the reaction 
proceeds through a hydrolytic pathway during the first step of 
splicing, protamine promotes a high level of exon ligation 
(Fig. 27, lanes 18-20, Table 7), without a concommitant SER 
10 side reaction. Thus, protamine seems to stimulate a stable 
conformation of the intron at low salt and, while it is not the 
most reactive ionic condition analyzed, it promotes clean self- 
splicing without side reactions. 

15 Relationship of results to previous studies: 

Self -splicing by group II introns was first evident in 
reactions showing the evolution of two intronic product bands 
(Peebles et al., 1986; Schmelzer and Schweyen, 1986; van der 
Veen et al., 1986). The slower migrating intronic product was 

20 shown to be a lariat molecule. Initial studies of the fast- 
migrating linear species concluded that, except in KC1, it 
consisted primarily of "broken lariat" (Peebles et al., 1986; 
Peebles et al., 1987; Jarrell et al., 1988b). This is a 
reasonable conclusion because even a single break in lariat 

25 intron would result in an RNA with the same mobility as linear 
intron. However, a hydrolytic pathway for production of 
linear intron was never discounted and several pieces of 
evidence indicated that linear intron might be significant 
reaction product: First of all, the "broken lariat" band was 

30 refractory to debranching treatment by HeLa cell S100 extract 
(Peebles et al., 1986). Second, the lariat product, while 
easily degraded during the isolation procedure, was shown to be 
stable to reincubation under splicing conditions (Peebles et 
al., 1986), indicating that any conversion of lariat band to 

35 "broken lariat" was a process too slow to explain the amount 
of linear intron generated under the reaction conditions. This 
result also indicated that linear species arising solely by 
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"self -debranching" is also unlikely, a result that has been 
confirmed by later work on the relative stability of the 2' -5' 
linkage (Chin and Pyle, 1995), 

5 Additional biochemical data supported the notion that the fast- 
migrating intronic species contained significant quantities 
linear intron. The molecular nature of the linear 
intron/broken lariat band was investigated by Tl ribonuclease 
fingerprinting and by branchpoint trinucleotide analysis with 

10 nuclease PI (Peebles et al., 1987). Data from these studies 
indicates that a fraction of the "broken lariat" band consisted 
of linear species. Quantitation of the branchpoint 

trinucleotides revealed that approximately 80% of the 
linear/broken lariat band seen in KC1 reactions was linear 

15 intron, while 20% was broken lariat. In ammonium sulfate, the 
raw data suggest that approximately 50% of this band was linear 
intron after 30 minutes (Peebles et al., 1987). Jarrell et al . 
concluded that hydrolytic cleavage at the 5' splice site is a 
bona fide pathway utilized by the group II intron under KC1 

20 conditions in vitro (Jarrell et al., 1988b). 

That a hydrolytic pathway could mediate self -splicing was 
particularly evident from experiments studying the role of the 
branchpoint A, which showed that a bulged adenosine is 

25 dispensible for self -splicing in vitro (van der Veen et al . , 
1987) . Self splicing can occur under (NH 4 ) 2 S0 4 conditions in 
vitro when the branchpoint A is either deleted or basepaired in 
a helix. In both cases complete splicing proceeds (albeit more 
slowly) via an apparently linear pathway, since only a linear 

30 intron product was seen and no lariat band could be identified. 
The possibility that rapid cleavage converted all lariat 
product into "broken lariat" was discounted by isolating the 
putative linear band, digesting with T2 ribonuclease, and 
isolating the trinucleotide on two-dimensional TLC (van der 

35 Veen et al., 1987). No branchpoint trinucleotide was recovered 
in the basepaired mutant. Splicing by an intron without a 
functional branch-point has also been observed in our 
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laboratory during mutational studies on the bulged adenosine 
When the branch-point A is replaced by a C, lariat molecules 
cannot be detected. However, hydrolytic splicing (in both 0 5 
M <NH 4 ) 2 SO 4 /100 mM Mg>- and 100 mM Mg* conditions) proceeds with 
5 the same rates as reported here for the wild-type intron 
(Qiaolian Liu and A.M. P., manuscript in preparation). This 
provides further evidence for a linear pathway and confirms 
that, using the parallel kinetic mechanism on the wild-type 
intron, hydrolytic rate constants were correctly deconvoluted 
10 from branching rate constants in the present study (Table 5) . 

This study does not contradict previously published data, but 
extends understanding of the first step of splicing in several 
ways. The data indicate that a hydrolytic pathway for splicing 
15 competes directly with branching at all times, under all 
conditions tested. Certain salts or proteins can 

preferentially promote one pathway or the other, and either 
route can be fast. Parallel reaction kinetics are followed, 
and linear intron should not necessarily be considered a 
20 secondary or degradative product. Hydrolytic cleavage in the 
first step of splicing does not prevent the second step 
although rate reductions in the second step cannot be ruled out 
since the first step remains rate-limiting under all 
conditions. We find that hydrolytic splicing is no more likely 
25 than branch-mediated splicing to be accompanied by the 
destructive SER reaction. The breadth of conditions analyzed 
nere indicate that the hydrolytic pathway is not an artifact of 
high salt reaction conditions, because the linear pathway can 
proceed, and even dominate, under low salt conditions similar 
30 to physiological. Given that a group II intron has the 
chemical capability to react through two distinctive pathways, 
and that the reaction environment regulates the choice of 
dominant pathway, it is important to consider the potential 
involvement of either pathway in the many important biological 
35 functions of group II introns. 
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Plasmids and transcripts: The aiSg intron and flanking exons 
from Saccharonyces cerevisiae were cloned downstream from a T7 
RNA polymerase promoter, and are contained in plasmid P JD20 
(jarrell et al., 1988a). Unspliced precursor RNA was 
5 transcribed from plasmid pJD20 after linearization by 
restriction enzyme Hind III. The precursor RNA is 1497 
nucleotides in length, consisting of 887nt of intron flanked by 
292 nucleotides of 5' -exon and 318 nucleotides of 3' -exon. 
Using methods described previously (Pyle and Green, 1994), 
10 transcriptions were performed in the presence of or- P-UTP. 
resulting in precursor RNA that was internally labelled 
throughout the molecule. Although a small amount of RNA 
splicing takes place during transcription (the buffer contains 
10 mM . the major product was precursor RNA. Transcription 
15 was stopped after 1.5 hours by adding an equal volume (usually 
20ul> of gel loading dye (10M urea, 0.1% xylene cyanol and 
bromophenol blue dyes, 40mM tris pH 7.5, 8.3% sucrose, and 
0 83mM EDTA, final concentrations) and then loading directly on 
a 4% denaturing polyacrylamide gel. Following electrophoreses, 
20 the band was visualized by autoradiography for 2 minutes, 
excised from the gel, and then eluted overnight in MOPS storage 
buffer (MOPS SB: lOmM MOPS, imM EDTA, pH 6.5) . RNA was then 
isolated by ethanol precipitation of the eluate (3 volumes of 
100% ethanol. 0.3M sodium chloride) and then washed with one 
25 volume of 70% ethanol. The RNA was dried, resuspended in 20ul 
of MOPS SB, and stored at -20°C. RNA concentrations were 
determined from their specific activities. 

Splicing Reactions: The splicing reactions were performed 
30 using a total volume of 20 to 40 pi, in 0.65 ml Slickseal 
tubes The concentration of precursor RNA for the experiments 
reported here was 1 nM, although all reactions were repeated at 
0 1 nM precursor to observe any concentration effects that 
might result from a sequential, rather than parallel reaction 
35 pathway. No differences in rate were observed. Reactions were 
initiated by preincubating the precursor RNA (in MOPS SB) at 
95i C for one minute. The RNA sample was then allowed to cool 
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to reaction temperature before combining with an equal volume 
of a salt mix, such that the final ionic reaction conditions 
were as indicated (text, Figure 26, Table 5). The various 
ionic conditions studied were 10 mM MgCl 2 (Aldrich, 99,99%) 
5 with 2 mM spermidine or 5 nM protamine (Sigma) ; 100 mM MgCl 2 ; 
and 100 mM MgCl 2 with 500 mM KCl (99.99% pure, Aldrich) , 500 mM 
NH 4 C1 or (NH4) 2 S0 4 . Before use in splicing reactions, protamine 
was purified by gel -filtration chromatography. Regardless of 
ionic condition, each sample was buffered with 4 0 mM MOPS, pH 

10 7.5 and the reactions were incubated at 42 jC. To measure 
extent of reaction at specific timepoints, aliquot s (1 or 2 ul) 
were removed from the reaction mixture at the specified time 
intervals. The aliquots were placed in equal volume of quench 
buffer (1.8% sucrose, IX TBE, 0.018% xylene cyanol dye, 36% 

15 formamide, and 25mM EDTA) and placed on ice to stop the 
reaction. Samples were then loaded onto a 4% denaturing 
polyacrylamide gel and products were resolved by 
electrophoresis. The gels were dried and quantitated on a 
Packard Instant imager , a radioanalytic scanner. 

20 

For the low-salt splicing reactions conducted solely in 100 mM 
Mg 2 *, or in 10 mM Mg 2 * with spermidine or protamine, the type of 
base used in pH adjustment of the MOPS buffer had significant 
effects on the reaction. Under the above conditions, when pH 

25 of the reaction buffer was adjusted with either KOH or NH 4 0H, 
reaction products were indistinguishable and the effects 
described in Fig. 3 and in the text were observed. However, 
when pH was adjusted with NaOH, hydrolytic cleavage at the 5'- 
splice site was inhibited several -fold (provided no other 

30 monovalent salts were added) . Because pH adjustment introduces 
small quantities of particular monovalent cations that may have 
inhibitory or stimulatory effects, it is important to consider 
the effects of additional components in the reaction buffer 
while measuring splicing at low salt. 

35 

Quantitative analysis of splicing reactions: Radioactivity in 
individual bands was quantitated on the radioanalytic scanner. 
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Counts from a nearby blank position on the same lane were 
subtracted (after normalizing for decreasing counts/mobility in 
the t = 0 lane) as an internal measure of background RNA 
degradation at each timepoint . Thi9 was facilitated by running 
5 separation gels {4% acrylamide) extensively, in order to 
cleanly separate all products. The relative fractions of 
precursor, lariat, linear, and exonic products were plotted 
against time and fit to equations 1-3 or 6-8, as shown in the 
text using Kaleidagraph (Abelbeck Software) . The reaction of 

10 unspliced precursor RNA and evolution of both of its two 
products (lariat and linear intron) were each fit individually* 
Under all reaction conditions, self -splicing proceeded to 
completion (Fig. 27A, lane 15), thus eliminating the need for 
endpoint correction of the data. Equations used for fitting 

15 the data are presented as follows: 

For a single population of molecules reacting through parallel 
kinetics, the independent evolution of two products (lariat and 
linear intron) is described by the following equations, which 
20 were applied in analysis of data from the 100 mM Mg a * and 10 mM 
Mg 2 *, 6 nM protamine reaction conditions (Fig. 24A, B) . 



25 



Prequrspr (t) 
Precursor,, 



-kTc 



Lj.ngar (t) 
Precursor 0 



30 



35 



40 



Lariat tt) 
Pre curs or 0 



where : 



(1 - e' kn ) 



= f br {l - e- m ) 



*by 



X k-r 

f*r_ X k T 

fhy + fbr 



(1> 

(2) 
(3) 
(4) 
(5) 



khy + kj 



br stotal rate of reaction 
fraction of linear intron 
fraction of lariat intron 



It is theoretically correct that + f tor - 1, thus simplifying 
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equation 4 to » f hy • k T . However, it is more accurate to 
use the individual values for f in the denominator, as in 
equations 4 and 5. In so doing, one does not presume that the 
reaction goes to completion (a dangerous assumption since RNA 
5 degradation at long times might obscure a fraction of unreacted 
material) . In addition, if the parallel model for splicing is 
correct, then the sum of f^ and £ r will approach 1, thus 
providing an independent check on appropriate application of 
the kinetic model and the quality of the data (Table 5) . 

10 

For reactions that followed biphasic kinetics (reactions in 
high concentration of monovalent ion, Figures 26A-C) , equations 
1-5 were adapted to account for two reactive populations of 
molecules : 

15 Precursor (t) » f A • e'Vt + ** # e 'Vt «> 
Precursor 0 

Un^r(t) = f\ y e-y c ) + f%(l-e"V t ) (7) 

Precursor 0 

20 

La r lat lfc J = f\r e Vt) + f B br (l-e'Vt) (8) 

Precursor 0 

where : 

25 

kr 1 = k h /+ k br * ■ total reaction rate of population 
i ( i = A, B) . 

f 1 « fraction of population 1 <i = A, B) . 

30 f h /= final fraction of linear product from population 

i (i = A, B) . 

f bl >= final fraction of lariat product from population 
i (i = A, B) . 

35 

The experimental data in Fig. 26A-C were fit independently to 
equation 6-8 in order to obtain the variables k A T , f 1 ^ , and 
f* br , where i = A or B. As in the simplest case presented with 
equations 1-3, the evolution of intronic products and the 
4 0 decrease in precursor were fit independently. That all four 
fractional values f A hy , f A te , t\, and f 8 * add up to one and all 
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k* T values for a given i are generally the same indicates that 
a parallel reaction model is sufficient to explain the data 
(Table 5) . After this treatment, one can calculate the 
individual rates of branching or hydrolysis from each 
5 population (k 1 ^ and k 1 ^, where i « A or B) , using equations 9 
and 10 

kV (f\ y / f \y + (9) 

k\>r= (f\r/f X hy + f\r>* (10) 

10 

Quantitation of branched nucleotide content using PI nuclease: 
Plasmid pJD20 was transcribed as described previously, except 
that the amount of labeled nucleotide incorporated was 
increased in order to obtain a higher signal for any branched 

15 trinucleotide. Specifically, 250 uCi of a 1:1 (ethanol : water) 
mix of or- 32 P-GTP (Amersham) was evaporated to dryness and the 
components for a standard body-labeling transcription reaction 
were added directly into this tube. The components were: 1 jxg 
linearized pJD20 template DNA, 2 jrl of 10X transcription 

20 buffer, 0.5 mM each ATP, CTP and OTP; 0.125 mM GTP, 1 jil T7 RNA 
polymerase (2000 u/jil) and H 2 0 to 20 /il . After vortexing this 
solution to dissolve the dried a- 32 P-GTP in the tube, the 
reaction was allowed to incubate for 2 hours at 37©c before 
loading directly on a 4% denaturing polyacrylamide gel. After 

25 electrophoresis, the full-length precursor band was visualized 
by autoradiography, isolated, crushed, and the RNA was eluted 
overnight at 4©c. The eluted RNA was recovered by ethanol 
precipitation. 

30 In order to collect splicing products, 5 x 10 s cpm of the full 
length transcript was incubated under three of the standard 
high-salt splicing conditions described previously (either 0.5M 
(NH 4 ) 2 S0 4/ 0.5M NHCl, or 0.5M KCl, all at lOOmM MgQl , 40mM 
MOPS pH 7,5) . Reactions were incubated at 42°C and timepoints 

35 were taken at 30 minutes (0.5M (NH 4 ) 2 SO,, and 0.5M NHCl), 40 
minutes (0.5M KCl), and 2 hours (all conditions). Precursor, 
lariat and putative linear products were eluted after 
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electrophoresis of the reactions on a 4% denaturing 
polyacrylamide gel. For each sample of product or precursor, 
the content of branched trinucleotide was analyzed by digestion 
with PI nuclease (Reilly et al. ff 1989) in a buffer containing 
5 25mM sodium acetate pH 5.3, 0.025 pg/ul tRNA, and 0.005 
units/ul PI nuclease (Boehringer) . Reactions were incubated 
for 10 minutes at 37 jC before adding equal volume of quench 
buffer (36% formamide, 0.1 X TBE, 1.8% sucrose, 0.02% xylene 
cyanol dye, 0.02% bromophenol blue dye, 0.02% orange G dye and 
10 25 mM EDTA final concentrations) . Samples were loaded on a 
denaturing 20% polyacrylamide gel. Mononucleotide pG was 
separated from the branched trinucleotide pA(2'-pG)-pU by 
electrophoresis (Chin and Pyle, 1995); gels were then dried and 
quant itated on a Packard Instantimager . 

15 
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WHAT IS CLAIMED IS: 

1. A composition for catalyzed oligonucleotide cleavage 
comprising a synthetic non-naturally occurring 
oligonucleotide compound which compound comprises 
nucleotides whose sequence defines a conserved group II 
intron catalytic region and nucleotides whose sequence is 
capable of hybridizing with a predetermined oligonucleotide 
target sequence to be cleaved, such target sequence not 
being present within the compound, and an appropriate 
oligonucleotide co-factor. 

2. The composition of claim 1, wherein the conserved group II 
intron catalytic region is a group II intron domain I 
catalytic region. 

3. The composition of claim 2, wherein the conserved group II 
intron domain I catalytic region further comprises a 
conserved portion of a group II intron domain II, a group 
li intron domain III, a group II intron domain IV, a group 
II intron domain V, or a group II intron domain VI in cis 
or trans. 

4. The composition of claim 1, wherein the nucleotides whose 
sequence is capable of hybridizing with a predetermined 
oligonucleotide target sequence to be cleaved comprises two 
hybridizing regions each region having 2 to 12 nucleotides 
capable of hybridizing with the oligonucleotide target to 
be cleaved. 

5. The composition of claim 4, wherein each hybridizing 
regions has 6 or 7 nucleotides. 

6. The composition of claim 1, wherein the cleavage is due to 
transesterif ication . 

7. The composition of claim 1, wherein the synthetic non- 
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naturally occurring oligonucleotide compound comprises 
between about 80 and about 1000 nucleotides. 

8. The composition of claim 1 wherein the appropriate co- 
factor is an oligonucleotide compound from a conserved 
portion of a group II intron domain V having the formula: 

3' (X) n (X>„. 5' 
X * X 
X t A 
X * G 
G + C 
(X) c * (X) c , 
G ♦ U (X) d . X 

(X) a * A 
(X) d A X 

wherein each X represents a ribonucleotide which may be the 
same or different; wherein each of (X) n , (X) n . , {X\ , , 
(X) d ., (X) a and (X) d represents an oligonucleotide; wherein n, 
n' , c, c' , d, d' , and a each represents an integer which 
defines the number of nucleotides in the oligonucleotide with 
the provisos that n and n' are greater than or equal to 1; a 
represents an integer which is greater than or equal to 1; d 
and d' represent an integer which is greater than or equal to 
5; c and C represents an integer which is greater than or 
equal to 4; wherein each * represents base pairing between the 
nucleotide located on either side thereof; wherein each t may 
or may not represent base pairing between the nucleotide 
located on either side thereof; and wherein each solid line 
represents a chemical linkage providing covalent bonds between 
the nucleotide located on either side thereof. 
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9. The composition of claim 1 wherein the appropriate co- 
factor is an oligonucleotide compound from a conserved 
portion of a group II intron domain V having the formula: 

3' (X) n <X) n . 5' 
X * X 
X t X 
X * G 
G * C 
(X) c * <X) C . 
G * U (X) d , X 

(X), * A 
<X) d A X 

wherein each X represents a ribonucleotide which may be the 
same or different; wherein each of (X) n , (X) n , , (X) c , (X) c . , 
(X) d ,, (X) a and (X) d represents an oligonucleotide; wherein 
n, n', c, c' , d, d' , and a each represents an integer which 
defines the number of nucleotides in the oligonucleotide 
with the provisos that n and n' are greater than or equal 
to 1; a represents an integer which is greater than or 
equal to 1; d and d' represent an integer which is greater 
than or equal to 5; c and c' represents an integer which is 
greater than or equal to 4; wherein each * represents base 
pairing between the nucleotide located on either side 
thereof; wherein each t may or may not represent base 
pairing between the nucleotide located on either side 
thereof; and wherein each solid line represents a chemical 
linkage providing covalent bonds between the nucleotide 
located on either side thereof. 

10. The composition of claim 1 further comprising a divalent 
cation. 

11. The divalent cation of claim 10, wherein the divalent 
cation is Mg 2 * or Mn 1+ . 

12. The composition of claim 1, wherein the predetermined 
oligonucleotide target sequence to be cleaved is mRNA. 

13. The composition of claim 1, wherein the predetermined 
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oligonucleotide target sequence to be cleaved is DNA. 

14. The composition of claim 13, wherein the mRNA encodes a 
growth factor. 

15. The composition of claim 14, wherein the growth factor is 
an angiogenic factor, a basic fibroblast growth factor, a 
colony-stimulating factor 1, cystic fibrosis transmembrane 
conductance regulator, an epidermal growth factor, an 
erythropoietin, a fibroblast growth factor, a G-protein, 
a granulocyte -macrophage colony stimulating factor, a 
growth hormone, IL-l, IL-1R, IL-2, IL-2R, il-4, il-6, an 
insulin- like growth factor, an insulin- like growth factor 
1, an interferon, an interleukin, a keratinocyte growth 
factor, luteinizing hormone receptor, MDR-1, a nerve 
growth factor, a platelet derived growth factor, a scatter 
factor, a transforming growth factor a, a transforming 
growth factor £, a transforming growth factor, or a tumor 
necrosis factor. 

L6. The composition of claim 12, wherein the mRNA encodes an 
oncogene or a tumor suppressor gene. 

7. The mRNA of claim 16, wherein the oncogene or tumor 
suppressor gene is Jbcl-2, bcr-abl t bek, BPV, c-abl, c-fes, 
c-fms, c-fos, c-H-ras, c-Jcit, c-myb, c-myc, c-mos, c-sea, 
cerbB, DCC, erbA, erbB-2, ets, fig, FSFV gp55, Ha-ras, 
HIV tat, HTLV-1 tat, JCV early, jun, L-/nyc, lc/c, LPV 
early, met, N-myc, NF-i, N-ras, neu, p53, Py mTag, pira-l, 
ras, RB, rel, retinoblastoma- 1 , SV-40 Tag, TGF-a, TGF-0, 
trk, trkB, v-abl, v-H-ras, v-jun, or WT-1. 

8. The composition of claim 17, wherein the mRNA is an mRNA 
associated with a chromosomal translocation. 

9. The composition of claim 17, wherein the oncogene or tumor 
suppressor gene has one or more point mutations. 
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20. The composition of claim 12 , wherein the mRNA is an mRNA 
whose overproduction is associated with a disease or 
condition. 

21. The composition of claim 12, wherein the mRNA is a 
mammalian mRNA. 

22. The composition of claim 12, wherein the mRNA is a yeast 
mRNA. 

23. The composition of claim 12, wherein the mRNA is a 
bacterial mRNA. 

24. The composition of claim 12, wherein the mRNA is a viral 
mRNA. 

25. The composition of claim 24, wherein the viral mRNA is 
associated with a cytomegalovirus, an Epstein-Barr virus, 
a hepatitis B virus, a hepatitis C virus, a herpes simplex 
virus, a herpesvirus, a HIV-1 virus, an immunodeficiency 
virus, an influenza virus, a papillomavirus, a 
picornavirus, a polio virus or a T-cell leukemia virus. 

26. The composition of claim 12, wherein the mRNA is a plant 
mRNA. 

27. The composition of claim 26, wherein the plant mRNA is an 
mRNA associated with alfalfa, apples, asparagus, bananas, 
broccoli, carrots, celery, chicory, coffee, cabbage, 
mustard, com, cottonseed, squash, cucumber, cantaloupe, 
grapes, lettuce, palm, potato, rapeseed, raspberry, 
soybean, sunflower, strawberry, tomato, or wheat. 

28. An oligonucleotide transfer vector containing a nucleotide 
sequence which on transcription gives rise to the 
synthetic non-naturally occurring compound of claim 1. 
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29. The transfer vector of claim 28, wherein the transfer 
vector is a bacterial plasmid, a bacteriophage DNA, a 
cosraid, or an eukaryotic viral DNA. 

30. A host cell transformed by the transfer vector of claim 

28. 

31. The host cell of claim 30, wherein the host cell is a 
prokaryotic host cell or an eukaryotic host cell. 

32. The prokaryotic host cell of claim 31, wherein the 
prokaryotic host cell is an E. coli host cell. 

33. The eukaryotic host cell of claim 31, wherein the 
eukaryotic host cell is a monkey COS host cell, a Chinese 
hamster ovary host cell, a mammalian host cell or a plant 
host cell. 

34 . A method for cleaving a target oligonucleotide comprising 
the steps of: 

a) composition for catalyzed oligonucleotide cleavage 
comprising a synthetic non-naturally occurring 
oligonucleotide compound which compound comprises 
nucleotides whose sequence defines a conserved group 
II intron catalytic region and nucleotides whose 
sequence is capable of hybridizing with a 
predetermined oligonucleotide target sequence to be 
cleaved, such target sequence not being present 
within the compound, and an appropriate 
oligonucleotide co- factor and 

b) contacting the synthetic non-naturally occurring 
oligonucleotide compound with a target 
oligonucleotide to cause cleavage of the target 
oligonucleotide . 
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35. The method for cleaving an oligonucleotide of claim 34, 
wherein the conserved group II intron catalytic region 
further comprises a group II intron domain II, a group II 
intron domain III, a group II intron domain IV, or a group 
II intron domain VI. 

36. The method for cleaving an oligonucleotide of claim 34, 
wherein the nucleotides whose sequence is capable of 
hybridizing with a predetermined oligonucleotide target 
sequence to be cleaved comprises two hybridizing regions 
each region having 2 to 10 nucleotides, 

37. The two hybridizing regions of claim 36 having 6 and 7 
nucleotides respectively. 

38. The method for cleaving an oligonucleotide of claim 34 
wherein the composition further comprises a divalent 
cation. 

39. The divalent cation of claim 38, wherein the divalent 
cation is Mg 24 or Mn 2 *. 

40. The method for cleaving an oligonucleotide of claim 34, 
wherein the predetermined oligonucleotide target sequence 
to be cleaved is mRNA. 

41. A transgenic, non-human vertebrate animal, having one or 
more cells bearing a DNA sequence encoding the composition 
of claim 1 . 

42. The transgenic animal of claim 41, wherein the animal is 
a mammal . 

43. The transgenic mammal of claim 42, wherein the mammal is 
selected from the group consisting of cows, goats, sheep, 
pigs, horses, dogs, cats, and rodents. 
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44. The transgenic animal of claim 41, wherein the animal is 
a fish, a chicken, or a turkey. 

45. The transgenic animal of claim 41, wherein the target 
sequence is a viral target sequence. 

46. The viral target sequence of claim 45, wherein the viral 
target sequence is a target sequence from a picornavirus, 
an immunodeficiency virus, a hepatitis B virus, a 
papillomavirus, an Epstein-Barr virus, a T-cell leukemia 
virus, a hepatitus C virus, a cytomegalovirus, an 
influenza virus, a herpes simplex virus, or a herpesvirus. 

47. The transgenic animal of claim 41, wherein the target 
sequence is a target sequence from a mRNA encoding a 
growth factor. 

48. The transgenic animal of claim 47, wherein the growth 
factor is an insulin-like growth factor, insulin-like 
growth factor 1, a transforming growth factor, 
transforming growth factor 0, transforming growth factor 
a, an interferon, a tumor necrosis factor, a platelet 
derived growth factor, a nerve growth factor, an epidermal 
growth factor, a fibroblast growth factor, a basic 
fibroblast growth factor, colony-stimulating factor 1, an 
interleukin, a granulocyte-macrophage colony stimulating 
factor, erythropoietin, a HGF, a keratinocyte growth 
factor, a scatter factor, or an angiogenic factor. 

49. The transgenic mammal of claim 41, wherein the target 
sequence is a target sequence from a mRNA encoding an 
oncogene or a tumor suppressor gene. 

50. The mRNA of claim 49, wherein the oncogene or tumor 
suppressor gene is bci-2, bcr-ab2, bek, BFV, c-ab2, c-fes, 
c-f/ns, c-fos, c-H-ras, c-kit, c-myb, c-jnyc, c-mos, c-sea, 
cerbB, DCC, erbA, erbB-2, ets, fig, FSFV gp55, Ha -ras, 
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FIGURE 2B 
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FIGURE 4B 
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FIGURE 11 A 
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FIGURE 12A 
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FIGURE 12B 
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FIGURE 12C 
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FIGURE 13A 
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FIGURE 14B 
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FIGURE 14C 
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